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Abstract. The top quark is the heaviest elementary particle observed to date. Its large mass makes 
the top quark an ideal laboratory to test predictions of perturbation theory concerning heavy quark 
production at hadron colliders. The top quark is also a powerful probe for new phenomena beyond 
the Standard Model of particle physics. In addition, the top quark mass is a crucial parameter 
for scrutinizing the Standard Model in electroweak precision tests and for predicting the mass of 
the yet unobserved Higgs boson. Ten years after the discovery of the top quark at the Fermilab 
Tevatron top quark physics has entered an era where detailed measurements of top quark properties 
are undertaken. In this review article an introduction to the phenomenology of top quark production 
in hadron collisions is given, the lessons learned in Tevatron Run I are summarized, and first Run II 
results are discussed. A brief outlook to the possibilities of top quark research a the Large Hadron 
Collider, currently under construction at CERN, is included. 
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1. Introduction 

The top quark is the by far heaviest of the six fundamental fermions in the Standard Model 
(SM) of particle physics. Its large mass made the search for the top quark a long and tedious 
process, since accelerators with high centre-of-mass energies are needed. In 1977 the discovery of 
the bottom quark indicated the existence of a third quark generation, and shortly thereafter the 
quest for the top quark began. Searches were conducted in electron-positron (e+e~) and proton- 
antiproton {pp) collisions during the 1980s and early 1990s. Finally, in 1995 the top quark was 
discovered at the Fermilab Tevatron pp collider. Subsequently, its mass was precisely measured to 
be Mtop = (178.0 ± 4.3) GeV/c^ [T]. The relative precision of this measurement (2.4%) is better 
than our knowledge of any other quark mass. The top quark is about 40 times heavier than the 
second-heaviest quark, the bottom quark. Its huge mass makes the top quark an ideal probe for new 
physics beyond the SM. It remains an open question to particle physics research whether the observed 
mass hierarchy is a result of unknown fundamental particle dynamics. It has been argued that the top 
quark could be the key to understand the dynamical origin of how particle masses are generated by the 
mechanism of electroweak symmetry breaking, since its mass is close to the energy scale at which the 
break down occurs (vacuum expectation value of the Higgs field = 246 GeV) [2] . The most favoured 
framework to describe electroweak symmetry breaking is the Higgs mechanism. The masses of the 
Higgs boson, the W boson and the top quark are closely related through higher order corrections to 
various physics processes. A precise knowledge of the top quark mass together with other electroweak 
precision measurements can therefore be used to predict the Higgs boson mass. 

At present, top quarks can only be directly produced at the Tevatron. The physics results of 
the Tevatron experiments CDF and D0 will therefore be the focus of this article. We review the 
experimental status of top quark physics at the beginning of the Tevatron Run II which will yield 
considerably improved measurements in the top sector. We discuss the first Run II analysis results 
and summarise the lessons learned from Run I data taken between 1990 - 1995. 

The outline of this article is as follows: In chapter [5] we give a brief introduction to the SM of 
particle physics and stress the importance of the top quark for higher order corrections to electroweak 
perturbation theory. In particular, wc discuss electroweak precision measurements used to predict 
the top quark and Higgs boson mass within the SM. Chapter [3] covers the theoretical description of 
SM top quark production at hadron colliders and the top quark decay. In chapter |4] we elaborate 
on analysis techniques used to detect top quarks in particle detectors. A short description of the 
Tevatron experiments CDF and D0 is provided. In chapter [5] we recall the early searches for the top 
quark in the 1980s and the discovery at the Tevatron in 1994/95. In chapter [6] we present various 
techniques to measure the top-antitop pair production cross section. Chapter [7] covers the top quark 
mass measurements and chapter [8] summarises various advanced tests of top quark properties. In 
chapter [5] we discuss the search for anomalous (non SM) top quark production. Previous reviews of 
top quark physics can be found in references [31 [H [U [HI El [HI [S] ■ 

2. The top quark in the Standard Model 

2.1. The standard model of particle physics 

The Standard Model (SM) of particle physics postulates that all matter be composed of a few basic, 
point-like and structureless constituents: elementary particles. One distinguishes two groups: quarks 
and leptons. Both of them are fermions and carry spin 1/2. The quarks come in six different flavours: 
up, down, charm, strange, top and bottom; formally described by assigning flavour quantum numbers. 
The SM incorporates six leptons: the electron (e~) and the electron- neutrino (i^e), the muon (/x^) and 
the muon-neutrino (t'p), the tau (t~) and the tau- neutrino {vt). They carry electron, muon and tau 
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quantum numbers. Quarks and leptons can be grouped into three generations (or families) as shown 
in table [1] which also contains the charges and masses of the particles. The three generations exhibit 
a striking mass hierarchy, the top quark having by far the highest mass. Understanding the deeper 
reason behind the hierarchy and generation structure is one of the open questions of particle physics. 
Each quark and each lepton has an associated antiparticle with the same mass but opposite charge. 
The antiquarks are denoted u, d, etc. The antiparticle of the electron is the positron (e+). 

The forces of nature acting between quarks and leptons are described by quantized fields. The 
interactions between elementary particles are due to the exchange of field quanta which are said 
to mediate the forces. The SM incorporates the electromagnetic force, responsible for example for 
the emission of light from excited atoms, the weak force, which for instance causes nuclear beta 
decay, and the strong force which keeps nuclei stable. Gravitation is not included in the framework 
of the SM but rather described by the theory of general relativity. All particles with mass or 
energy feel the gravitational force. However, due to the weakness of gravitation with respect to 
the other forces acting in elementary particle reactions it is not further considered in this article. The 
electromagnetic, weak and strong forces are described by so called quantum gauge field theories (see 
explanation below). The quanta of these fields carry spin 1 and are therefore called gauge bosons. 
The electromagnetic force is mediated by the massless photon (7), the weak force by the massive VF^, 
Mw = (80.425 ± 0.038) GeV/c^ [TD], and the Z", Mz = (91.1876 ± 0.0021) GeV/c^ [TU], and the 
strong force by eight massless gluons (g). Quarks participate in electromagnetic, weak and strong 
interactions. All leptons experience the weak force, the charged ones also feel the electromagnetic 
force. But leptons do not take part in strong interactions. A thorough introduction to the SM can be 
found in various text books of particle physics [HI [12l [131 [HI [H] • 

2.1.1. Electroweak interactions. In quantum field theory quarks and leptons are represented by 
spinor fields ^' which are functions of the continuous space-time coordinates . To take into account 
that the weak interaction only couples to the left-handed particles, left- and right-handed fields 
'S'l = 5(1 ~ 75) Siud 'i'R = i(l -I- 75) * are introduced. The left-handed states of one generation 
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are grouped into weak-isospin doublets, the right-handed states form singlets: 
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The weak-isospin assignment for the doublet is: up-type quarks (u,c,t) and neutrinos carry T3 = -|-^; 
down- type quarks (d,s,b), electron, muon and tau lepton have T3 = — i. In the original SM the right- 
handed neutrino states are omitted, since neutrinos are assumed to be massless. Recent experimental 
evidence [ini HZl [TH], however, strongly indicates that neutrinos have mass and the SM needs to be 
extended in this respect. 

The dynamics of the electromagnetic and weak forces follow from the free particle Lagrangian 
density 



£0 = i * 7^9,, * 

by demanding the invariance of £0 under local phase transformations: 

*L > QW<y-(xyT+ig' li{x)Y — ^ Q^g'l3{x)Y 



(1) 



(2) 



For historical reasons these transformations are also referred to as gauge transformations. In ([2]) the 
parameter ol{x) is an arbitrary three-component vector and T = (Ti, T2, T3)* is the weak-isospin 
operator whose components Ti are the generators of SU{2)l symmetry transformations. The index L 
indicates that the phase transformations act only on left-handed states. The matrix representations 
are given by = ^ where the Ti are the Pauli matrices. The Ti do not commute: [Ti, Tj] = i eijk Tk- 
That is why the SIJ (2) l gauge group is said to be non-Abelian. (3{x) is a one-dimensional function of x. 
Y is the weak hypercharge which satisfies the relation Q = T3 + Y/2, where Q is the electromagnetic 
charge. Y is the generator of the symmetry group U{1)y- Demanding the Lagrangian Co to be 
invariant under the combined gauge transformations of SU{2)l x see ([2]), requires the addition 

of terms to the free Lagrangian which involve four additional vector (spin 1) fields: the isotriplet 
{Wi ,W2 ^,f for SU{2)l and the singlet for C/(l)y. This is technically done by 
replacing the derivative df^ in Cq by the covariant derivative 



^a'\B^Y 



iW 

9 ■ W,, 



X W 



(3) 



The field tensors 



and = - d^B^. Since 



and adding the kinetic energy terms of the gauge fields. ^ 
W^^ and Bf,^ are given by W^^ = d^W^ - 9,.W^ 
the vector fields and B^ are introduced via gauge transformations they are called gauge fields 
and the quanta of these fields are named gauge bosons. For an electron-neutrino pair, for example, 
the resulting Lagrangian is: 
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(4) 



This model developed by Glashow \19l, Weinberg and Salam [5D1 [H] in the 1960s allows to describe 
electromagnetic and weak interactions in one framework. One therefore refers to it as unified 
electroweak theory. 
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2.1.2. The Higgs mechanism. One has to note, however, that Ci describes only massless gauge 
bosons and massless fermions. Mass-terms such as ^M^B^B^^ or —m^^ are not gauge invariant and 
therefore cannot be added. To include massive particles into the model in a gauge invariant way the 
Higgs mechanism is used. Four scalar fields are added to the theory in form of the isospin doublet 
$ — (fp^ where and (fP are complex fields. This is the minimal choice. The term Ch = 

— is added to Ci. The scalar potential takes the form — <^^<^ + \ ($^$)^. 

In most cases particle reactions cannot be calculated from first principles. One rather has to use 
perturbation theory and expand a solution starting from the ground state of the system which is in 
particle physics called the vacuum expectation value. The parameters /i and A can be chosen such 

that the vacuum expectation value of the Higgs potential V is different from zero: |$vac| ~ \J ~^ 1^1^ 
and thus does not share the symmetry of V . The scalar Higgs fields inside $ arc redefined such that 
the new fields, ^{x) — i£,i{x),£,2{x),^3{x)Y and H{x), have zero vacuum expectation value. When 
the new parameterization of $ is inserted into the Lagrangian, the symmetry of the Lagrangian is 
broken, that is, the Lagrangian is not an even function of the Higgs fields anymore. This mechanism 
where the ground states do not share the symmetry of the Lagrangian is called spontaneous symmetry 
breaking. As a result, one of the Higgs fields, the H{x) field, has acquired mass, while the other three 
fields, ^, remain massless [22l [23]. 

Applying spontaneous symmetry breaking as described above to the combined Lagrangian 
C2 = Ci + Ch and enforcing local gauge invariance of £2, makes three electroweak gauge bosons 
acquire mass. This is the aim of the whole formalism. One field, the photon field, remains massless. 
The resulting boson fields after spontaneous symmetry breaking are, however, not the original fields 

and B^ but rather mixtures of those: the — (W^ =F i the Z° and the photon field 

A,: 

Af^ \ ^ / cos Ow sm0w \ { 
Z^. J ^\ -sinOw cosBw j V 

The mixing angle 9w is the Weinberg angle defined by the coupling constants g'/g = tanOw While 
the original massless vector fields have two degrees of freedom (transverse polarizations), the new 
massive fields acquired a third degree of freedom, their longitudinal polarization. The longitudinal 
modes of the and the Z° are due to the disappearance of the ^ states from the theory. The total 
number of degrees of freedom is thus conserved. 

Spontaneous symmetry breaking also generates lepton masses if Yukawa interaction terms of the 
lepton and Higgs fields are added to the Lagrangian: 



(5) 



^lepton ^ < 

■^Yukawa — 



(6) 

Here the Yukawa terms for the electron-neutrino doublet are given as an example. Ge is a further 
coupling constant describing the coupling of the electron and electron-neutrino to the Higgs field. In 
this formalism neutrinos are assumed to be massless. 

Quark masses are also generated by adding Yukawa terms to the Lagrangian. However, for the 
quarks both, the upper and the lower member of the weak-isospin doublet, need to acquire mass. 
For this to happen an additional conjugate Higgs multiplet has to be constructed: $c = iT2$* = 
((/)° , —(/<")*. The Yukawa terms for the quarks are given by: 

i — 1 j — 1 

The Uj and dj are the weak eigenstates of the up-type (u, c, t) and the down-type {d, s, b) 
quarks, respectively. Couplings between quarks of different generations are allowed by this ansatz. 
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After spontaneous symmetry breaking the Yukawa terms produce mass terms for the quarks 
which can be described by mass matrices in generation space: (wi, M2, ws)/? Al" (wi, U2, ^3)^ and 
(di, c?2, rfa)^ M'^ {di,d2,d^)\ with M^j = |*vac| • Gij and Mfj = |*vac| ■ Gij. The mass matrices are 
non-diagonal but can be diagonahzed by unitary transformations, which essentiaUy means to change 
basis from weak eigenstates to mass eigenstates, which are identical to the flavour eigenstates u, c, 
t and d, s, b. In charged-current interactions {W^ exchange) this leads to transitions between mass 
eigenstates of different generations referred to as generation mixing. It is possible to set weak and 
mass eigenstates equal for the up- type quarks and ascribe the mixing entirely to the down-type quarks: 

d' \ / d\ / Vud Vus Vub\ / d\ 

s' = V s ^ I Vcd Vcs Vcb ] { s ] (8) 

b' J, V Wl V v,b J \ b 

where d' , s' and b' are the weak eigenstates. The mixing matrix V is called the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix [24] . 

2.1.3. Strong interactions. The theory of strong interactions is called quantum chromodynamics 
(QCD) since it attributes a colour charge to the quarks. There are three different types of strong 
charges (colours): "red", "green" and "blue". Strong interactions conserve the flavour of quarks. 
Leptons do not carry colour at all, they are inert with respect to strong interactions. QCD is a 
quantum field theory based on the non-Abelian gauge group SU{2>)c of phase transformations on the 
quark colour fields. Invoking local gauge invariance of the Lagrangian yields eight massless gauge 
bosons: the giuons. The gauge symmetry is exact and not broken as in the case of weak interactions. 
Each gluon carries one unit of colour and one unit of anticolour. The strong force binds quarks 
together to form bound-states called hadrons. There are two groups of hadrons: mesons consisting of 
a quark and an antiquark and baryons built of either three quarks or three antiquarks. All hadrons 
are colour-singlet states. Quarks cannot exist as free particles. This experimental fact is summarised 
in the notion of quark confinement: quarks are confined to exist in hadrons. 

2. 2. Model predictions of top quark properties 

In 1975 the tau lepton was the first particle of the third generation to be discovered [25]. Only two 
years later, in 1977, a new heavy meson, the T, was discovered |26| . It was quickly recognised to 
contain a new fifth quark, the b quark (T = bb). After these discoveries the doublet structure of 
the SM strongly suggested the existence of a third neutrino associated with the tau lepton and the 
existence of a sixth quark, called the top quark. The properties of the 6 quark and b hadrons were 
subject to extensive experimental research in the past 25 years. As a result, charge and weak isospin 
of the b quark are well established quantities: Qb ~ and T3 = — i. The charge Qb was first 
deduced from measurements of the cross section CThad for hadron production on the T resonance at 
the DORIS e+e^ storage ring [27l [28l [29]. The integral over crhad is related to the leptonic width 
Fee of the T: J (Thad dAf = 6 7r^ree/M|. (using the approximation that Fhad ~ Ttotai)- Nonrelativistic 
potential models of the T relate Tee to the charge of the constituent b quarks. 

The weak isospin of the b quark was obtained from the measurement of the forward-backward 
asymmetry ^fb in e+e^ — > bb reactions. The asymmetry is defined as the difference in rate of b 
quarks produced in the forward hemisphere (polar angle 9 > 90°) minus the rate of b quarks produced 
in the backward hemisphere {9 < 90°) over the sum of the two rates. In 1984 the JADE experiment 
measured ApB = (—22.8 ± 6.0 ± 2.5)% [3D] in excellent agreement with the value of -25.2% predicted 
for a b quark that is the lower member of a SM weak isospin doublet. In case of a weak isospin singlet 
the asymmetry would be zero. The assignment of quantum numbers for the b quark allows to predict 
the charge and the weak isospin of the top quark to be: Qt = +| and = +i. 
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Figure 1. (a) An example of a fermion triangle diagram. (b,c) Two example Feynman diagrams 
of radiative corrections involving the top quark, (b) is a radiative correction to the propagator, 
(c) shows a vertex correction to the decay — > bh. 



Another argument which supported the existence of a complete third quark generation comes from 
perturbation theory. In particle physics the terms of a perturbation series are depicted in Feynman 
diagrams. The first order terms are pictured as tree level diagrams. Higher order terms correspond 
to loop diagrams. Certain loop diagrams are divergent. These divergencies can only be overcome by 
summing up over several divergent terms in a consistent manner and have the divergencies cancel each 
other. This formalism is called renormalization. For electroweak theory to be renormalizable, the sum 
over fermion triangle diagrams, such as the one shown in figure [T^, has to vanish [31, J2l[33]. With 
a third lepton doublet in place the cancellation only occurs if also a complete third quark doublet is 
there. While the SM predicts the charge and the weak isospin of the top quark, its mass remains a 
free parameter. 

2.3. Top quark mass and electroweak precision measurements 

Even though the top quark mass is not predicted by the SM it enters as a parameter in the calculation 
of radiative corrections to electroweak processes. With highly precise measurements at hand it 
is therefore possible to indirectly determine the top quark mass from those processes [Ml [3S]. In 
this context radiative corrections denote higher order contributions to a perturbation series, e.g. to 
calculate cross sections of electroweak processes. 

The most precise electroweak measurements are available from e+e~ colliders operating at the Z° 
pole, where -^/i = Mz. Here ^/s is the centre-of-mass energy of the colliding electrons and positrons. 
Mz is the mass of the Z° boson. In the 1990s there were two particle colliders operating on the pole: 
the Large Electron Positron Collider (LEP) [36] at CERN with four experiments (ALEPH [371I38]. 
DELPHI [391 HQ], L3 [mug and OPAL [13]) and the Stanford Linear Collider (SLC) f44^ with the 
SLD [45l|46] experiment. The LEP experiments collected about 4 million Z'^ events each, SLD on the 
order of 0.5 million events. Although the SLD sample is considerably smaller, the experiment reached 
a competitive sensitivity for many measurements, since at SLC it was possible to polarize the colliding 
electrons and positrons. 

Two examples of radiative corrections to the Z*^ propagator and to the vertex involving the 
top quark are given in figure [Th and[TJ:. The top mass plays a particular large role in the corrections 
represented by these diagrams due to the large mass difference between the top quark and its weak- 
isospin partner, the b quark. The correction terms introduce a quadratic dependence on the top mass, 
whereas the dependence on the Higgs mass is only logarithmic. The LEP electroweak working group 
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(LEPEWWG) has combined the measurements of the four LEP experiments to obtain a common data 
set. To test SM predictions and check the overaU consistency of the data with the model the SLD 
data and Tevatron measurements in pp colhsions are included |37] . Among the results of the fits is an 
indirect determination of the top mass and indirect limits on the Higgs mass. The quantities used for 
a constraint fit to the SM are discussed below: 

(i) The mass Mz and the total width Tz of the . The definition of Mz and Tz is based on the 
Breit-Wigner denominator of the electroweak cross section for fermion pair (//) production due 
to Z'^ exchange: cTff{s) — a° ■ sr|/((s — M|)^ + {sTz/Mz)"^) where a'j is the pole cross section 
at ^/s — Mz- 

(ii) The hadronic pole cross section: 

Q _ 12 TT Fee Fhad 

where Fee and Fhad are the partial widths for the Z° decaying into electrons or hadrons, 
respectively. 

(iii) The ratio of the hadronic and leptonic widths at the Z'^ pole: = Fhad/r^^- r« is obtained 
from Z° decays into e+ e~, /x"*" IjL~ or r+ t~ assuming lepton universality (Fee = = Trr = r«) 
which is only correct for massless leptons. To account for the tau mass Trr is corrected, in average 
by 0.2%. 

(iv) The pole forward-backward asymmetry Ap^ {£ e {e,/x,r}) for Z^ decays into charged leptons. 
The asymmetry is defined as the total cross section for a negative lepton to be scattered into 
the forward direction (90° < 9 < 180°) minus the total cross section for a negative lepton to be 
scattered into the backward direction (0° < < 90°) divided by the sum. 

^0,1 _ o-f(^-,90° <9< 180°) -crB(^~,0° <e'<90°) 
~ fTF(^-,90° <e< 180°) +ctb(^-,0° <e< 90°) 

Here 6 is the angle between the produced lepton and the incoming electron. The forward-backward 
asymmetry originates from parity violating terms in the total cross section. Apg can be written 
as a combination of effective couplings Ap^ = j-AeAe with Ai — (2 gv.e 9A,e) / (gy 9a i) where 
gv,e and gA,e are the vector and axial couplings, respectively. 

(v) The effective leptonic coupling Ai derived from the tau polarization measurement. Parity 
violation in the weak interaction leads to longitudinally polarized final state fermions from 
the Z° decay. In case of the decay — > r+r" this polarization can be measured via the 
subsequent parity violating decays of the tau as polarimeter. The polarization is defined in 
terms of the cross section for the production of right-handed and left-handed t^, respectively: 
Pt = {ctr — aL)/{<JR + cl). Averaging Pr over all production angles 9 gives a measure of the 
effective coupling: (P^) = — At- 

(vi) The effective leptonic weak mixing angle sin^ 9^^^^ from the inclusive hadronic charge asymmetry 
(Qfb) which is measured from Z° qq decays. An estimator for the quark charge is derived 
from the sum of momentum weighted track charges in the quark and antiquark hemispheres, 
respectively. 

(vii) The effective leptonic coupling Ae derived from left-right asymmetries measured with 
longitudinally polarized electron and positron beams at SLD. The left-right asymmetry ^lr 
is formed as the number of Z° bosons Nj^ produced by left-handed electrons minus the number 
of Z^ bosons Nn produced by right-handed electrons: 

Nl - Nr 1 



Nl + Nr Pe 
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Physical 
quantity 


Measurement 
with total error 


biVl nt 
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Measurement 
with total error 
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LEP measurements 




LEP and SLD heavy flavour 




Mz [GeV/c^] 


91.1875 ± 0.0021 


91.1875 


Rl 


0.21630 ± 0.00065 


0.21562 


Tz [GeV/c2] 


2.4952 ± 0.0023 


2.4966 


Rl 


0.1723 ± 0.0031 


0.1723 


al [nb] 


41.540 ± 0.037 


41.481 


^FB 


0.0998 ± 0.0017 


0.1040 




20.767 ± 0.025 


20.739 


40,c 
^FB 


0.0706 ± 0.0035 


0.0744 




0.0171 ± 0.0010 


0.0165 


At, 


0.923 ± 0.020 


0.935 


A ((Pr)) 


0.1465 ± 0.0033 


0.1483 


Ac 


0.670 ± 0.026 


0.668 


sin^^^J^ ((Qfb 


,)) 0.2324 ± 0.0012 


0.2314 








SLD only measurements 




UA2, Tevatron and LEP 2 measurements 


At 


0.1513 ± 0.0021 


0.1483 


Mw [GeV/c2] 


80.425 ± 0.034 


80.394 


Tevatron only measurements 




Tw [GeV/c2] 


2.133 ± 0.069 


2.093 


Mtop [GeV/c2] 


178.0 ± 4.3 


178.1 









Table 2. Experimental input to overall fits to the SM by the LEP electro- weak working group |47| . 
The experimental values of the physical quantities, which are described in detail in the text, are 
given with the total error. The results of the constraint fit to all data is given in the columns named 
"SM fit". In the fit the Higgs mass is treated as a free parameter. 



Additionally, one divides by the luminosity weighted polarization of the electron beam, ^lr 
is measured for all three charged lepton final states. The pole values A^pj are equivalent to the 
effective couplings Ae, A^ and Ar, respectively. Assuming lepton universality the results are 
combined to form a single value Ag. 

(viii) The measurement of the top mass by CDF and D0 [T]. 

(ix) The ratios of b and c quark partial widths of the with respect to its total hadronic partial 
width: Rl = Fbt/rhad and i?° = Vcs/^had- 

(x) The forward-backward asymmetries for the decays bb and Z^^ — s- cc: and Apg, defined 
analogously to the leptonic asymmetries. 

(xi) The effective couplings for b and c quarks: Ab and Ac, defined analogously to Ag. 

(xii) The mass and width of the W boson from the combination of measurements in pp collisions by 
UA2, CDF and D0 and measurements at LEP 2. 

The quantities describing the decay bb Ap^ and Ab) exhibit a particular strong dependence 

on Mtop because they are sensitive to weak vertex corrections as the one shown in figure [TJ;. In many 
cases vertex corrections involving a W boson are suppressed due to small CKM matrix elements. The 
vertices in figure[l] however, contain the matrix element Vtb which is close to 1. Therefore, the graphs 
lead to significant corrections depending on the top quark mass. 

The measured values of the above listed quantities are given in table. [51 Other input quantities 
not shown in the table are the Fermi constant Gp, the electromagnetic coupling constant at the Z° 
mass scale, a(Mf ), and the fermion masses. The Higgs mass and the strong coupling constant at the 
Z'^ mass scale, as(Mf ), are treated as free parameters in the fits. Detailed reviews of electroweak 
physics at LEP are, for instance, given in references [481149]. 
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Figure 2. Results of fits to clcctroweak data |47) . Direct Higgs boson searches at LEP2 and the 
top mass measurement at the Tevatron are summarised in these plots, (a) shows the plane "top 
mass versus Higgs mass" . The light shaded area is excluded by direct searches for the Higgs boson 
at LEP2. The dark shaded area is given by the central value of the direct top mass measurement at 
the Tevatron and its errors. The dashed line encloses the region preferred by the electroweak data 
at 68% confidence level, (b) shows the "Mw versus Mtop" plane. The full line encloses the area 
preferred by the SM fit to data from LEPl and SLD. The dashed contour indicates the result of the 
LEP2, UA2 and Tevatron W mass measurements and the direct top quark mass measurement. The 
plot also shows the SM relationship of the masses as function of the Higgs boson mass. 



If the top mass is left floating in the fit it is possible to obtain a prediction for Aftop from the 
model based on all other measurements. Comparing this prediction with the direct measurement from 
CDF and D0 is an important cross-check of the SM. The fit yields a value of A/top = 179+g^ GeV/c^ 
which is in very good agreement with the measured value of Aftop = (178.0 ± 4.3) GeV/c^. The result 
of the fit is depicted in figure [2^ where Aftop is plotted versus Mh- The plot shows that the overlap 
region of the direct measurement with the indirect determination prefers a low value for the Higgs 
mass. It is obvious from the diagram that there is a large positive correlation, about 70%, between 
the top quark and the Higgs boson mass. 

The constraint on the Higgs mass becomes significantly stronger if the measured value for Aftop 
is included in the SM fit. Only Mh and as{M'^) are free parameters in this case. The resulting value 
for the Higgs mass is: Mh = 114^45 GeV/c^, a value close to the exclusion limit obtained by direct 
searches at LEP2 which yield Mh > 114.4 GeV/c^ :50j. The SM fit (not using the direct search 
result) yields an upper limit of Mh < 260 GeV/c^ at the 95% confidence level. The SM prediction 
for all observables resulting from the constrained fit is given in table [51 

The last fit discussed here is the one made with Mw and Aftop left as free parameters. The result 
is shown in figure [2)3 which also includes SM predictions for Higgs masses from 114 to 1000 GeV/c^. 
The direct measurement of Mw and Aftop are in fair agreement with the indirect determination. A 
low Higgs mass is preferred. The dependence of the W mass on the top and Higgs mass is introduced 
via loop diagrams like the ones shown in figure [3] As becomes clear from figure [2] increasing the 
precision of the direct measurements of Af^y and Aftop is a precondition to gain higher leverage for 
SM predictions of the Higgs mass. 

Summarising this section on electroweak precision measurements we point out two major issues 
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Figure 3. Examples of loop diagrams of the W propagator introducing a dependence of the W 
mass on the top and Higgs mass. 



with respect to the top quark. First, the very good agreement of the predicted value for Mtop based 
on higher order electroweak corrections with the direct measurement is a significant success of the 
SM. Second, more precise measurements of Mtop and Mw are needed to obtain tighter constraints on 
the Higgs boson mass. This is one of the major physics motivations for the Run II of the Tevatron. 

2.4- The top quark in flavour physics 

Flavour physics describes the transitions between quarks of different flavour. These transitions always 
involve the exchange of a W boson. They are charged current interactions. Higher order transitions of 
this kind involve loops in which particles can occur that are much heavier than the hadrons involved 
in the interaction. Diagrams with down-type quarks (d, s and b) in the loops cancel each other via the 
GIM mechanism very effectively, since they are nearly degenerate in their masses. The large degree of 
mass splitting in the up-quark sector prevents the cancellation and leads to sizeable loop contributions, 
e.g. in radiative B meson decays. One other example where the top quark plays a prominent role is 
the mixing of neutral B mesons. 

The phenomenon of particle-antiparticle mixing has been experimentally established in the neutral 
kaon system {K^-K'^) and the B^~B^ system [HI [52]. Mixing is also expected to take place in the D°- 
DO and B°-'B2 system, but was not observed yet. Only limits were set p l [54 l l55 l [56 l [57 t l58 l [59 l [60 l [6T]. 

Mixing changes the flavour quantum number (strange, charm, bottom) of the mesons by two units 
{AS = 2, AC = 2, AB = 2). a neutral meson is produced in a well defined flavour eigenstate |P°) or 
|P°) with P E {K, D, Bd, Bg}. This initial state evolves due to second order weak interactions into a 
time-dependent quantum superposition of the two flavour eigenstates: \P{t)) = a{t) |P°) -|- b{t) |P°). 
The time-evolution of the coefficients a(t) and b{t) is given by the effective Hamiltonian matrix 
H = M-i/2r : 



M and T are 2x2 Hermitian matrices denoted as mass and decay matrix, respectively. CPT invariance 
requires Mu = A/22 = A/ and Fn = F22 — F, where M and F are the mass and the decay width of 

the flavour eigenstates. 

In case of the B meson systems {Bj^-B^ and B^-Bf) the transition amplitude for the mixing 
matrix is dominated by box diagrams involving the top quark, see figure SI All up-type quarks 
{q = u, c, t) can be exchanged in the box and contribute in general to the mixing amplitude. In the 
case of B mesons it is, however, found that the mixing is strongly dominated by the diagrams with 
two top quarks in the loop as shown in figure ID All other diagrams involving at least one u- or c-quark 
can be neglected with respect to the top-top diagram [62l [63l [64] . As a side remark: In the neutral 
kaon system the case is different. Charm and top quark box diagrams as well as intermediate virtual 
pion states contribute significantly to the mixing. For B mixing the off-diagonal elements of M and 




(9) 



Top quark physics in hadron collisions 



14 




q = d,s 



q = d, 



q = d,s 



t^' tit 



q = d,s 



Figure 4. Lowest order Feynman diagrams for B meson mixing. A (bd) or a {bs) oscillates 
into a (bd) or a (bs), respectively. 
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'i/2 Ti2 {q G {d, s}), are in very good approximation given by [62 t [64 l [10] : 



M- 



12,9 



12 7r2 



'"T'bJb.Bb^Vb ■ (Kg VtbY ■ So{xt) 



^12,q — 



q2 



+ v;^VtbV:^v,t o 



(10) 



(11) 



where Gp is the Fermi constant and Mw the W mass. rUb and rric are the h quark and c quark masses, 
respectively, tob^ represents the masses of the and B^ mesons. Jb^ is the decay constant and Bb^ 
the bag parameter introduced as a correction factor to hadronic matrix elements. The Vij are the CKM 
matrix elements, xt = {Mtop/MwY is the squared ratio of top quark mass over W mass. The function 
So{xt) is an Inami-Lim function [BS] and can be well approximated by So{xt) = 0.784 a;^ ''^ 
The parameters r]B and rj'g represent QCD corrections to the box diagram. The mass eigenstates 
Bh (h for heavy) and Bi (1 for light) diagonalize the effective Hamiltonian H and are given by: 
\Bh) =p\B°) -q\W) and \Bi) = p\B°) +q\W) with q"^ = M^^-^/^^h and = M12 -1/2^2- 
Mixing experiments do not determine the masses m/i and m; but rather the mass difference between the 
two Am = nih — mi which is in good approximation (about 1% accuracy) theoretically predicted to be 
Am = 2 IM12I [64] • This quantity depends on the top quark mass via r]BSo{xt)- To give a flavour of the 
dependency: r]BSo{xt) changes from roughly 1.1 for Aftop = 150 GeV to 1.7 at Mtop — 200 GeV [57] . 
Historically, the ARGUS measurement of B^-B^ mixing yielded a mass difference Am^ found to 
be surprisingly large [52j . Using the dependence of Am^ on the top quark mass several authors 
interpreted this measurement as a hint of a large top quark mass ^68, £9, ,70, ,71 , 72 . 
The dependence on Aftop drops out if the ratio 



Am^ 



Anid ruBJi^BB^ 



Vts 



Vtd 



(12) 



is taken. Once Am^ is measured in Bs mixing this relation will allow to extract the absolute value of 
the ratio Vts/Vtd with good precision. 



3. Top quark production at hadron colliders 

In this chapter we present the phenomenology of top quark production at hadron colliders. We limit 
the discussion to SM processes. Anomalous top quark production and non-SM decays will be covered 
in chapter [5] Specific theoretical cross section predictions refer to the Fermilab Tevatron, running at 
= 1.8 TeV (Run I) or Vs = 1-96 TeV (Run II), or to the future Large Hadron Collider (LHC) 
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at CERN {-s/s — 14 TeV). In the intermediate future the Tevatron and the LHC are the only two 
coUiders where SM top quark production can be observed. 

The two basic production modes of top quarks at hadron coUiders are the production of tt 
pairs, which is dominated by the strong interaction, and the production of single top quarks due to 
electroweak interactions. 

3.1. tt production 

We discuss only top quark pair production via the strong interaction, tt pairs can also be produced 
by electroweak interactions if a Z*^ or a photon are exchanged between the in- and outgoing quarks. 
However, at a hadron collider the cross sections for theses processes are completely negligible compared 
to the QCD cross section. The cross section for the pair production of heavy quarks has been calculated 
in perturbative QCD, i.e. as a perturbation series in the QCD running coupling constant as{y?). The 
results are applicable to the bottom and to the top quark. In the following we will refer only to tt 
production, the scope of this article. 

3.1.1. The factorization ansatz The underlying theoretical framework of the calculation is the parton 
model which regards a high-energy hadron A, in our case a proton or antiproton, as a composition 
of quasi-free partons (quarks and gluons) which share the longitudinal hadron momentum pA- The 
parton i has the longitudinal momentum pi, i.e. it carries the momentum fraction Xi = Pi/pA- The 
cross section calculation is based on the factorization theorem stating that the cross section is given 
by the convolution of parton distribution functions (PDF) fi{x,n'^) for the colliding hadrons {A, B) 
and the hard parton-parton cross section &ij : 

a{AB tt) J '^^i^^J fiM^i^l^^)fj.Bixj,ii^) ■ <Ty (ij tt] s,y?) (13) 

The hadrons AB are either pp (at the Tevatron) or pp (at the LHC). The parton distribution function 
fi,A(xi, ^'^) describes the probability density for finding a parton i inside the hadron A carrying a 
longitudinal momentum fraction Xi. 

The parton distribution functions and the parton-parton cross section dij depend on the 
factorization and renormalization scale ji. For calculating heavy quark production the scale is usually 
set to be of the order of the heavy quark mass, here specifically /i = Mtop. Strictly speaking one 
has to distinguish between the factorization scale introduced by the factorization ansatz and the 
renormalization scale due to the renormalization procedure invoked to regulate divergent terms 
in the perturbation series when calculating the parton-parton cross section Uij . Since both scales are 
to some extent arbitrary parameters most authors have adopted the practice to use only one scale 
ji — fiF — l^R in their calculations. If the complete perturbation series could be calculated, the 
result for the cross section would be independent of /x. However, since calculations are performed at 
finite order in perturbation theory, cross section predictions do in general depend on the choice of fj,. 
The /z-dependence is usually tested by varying the scale between /i = Mtop/2 and = 2Mtop- The 
variations in the cross section are quoted as an indicative theoretical uncertainty of the prediction, 
but should not be mistaken for Gaussian errors. 

In the variable s denotes the square of the centre-of-mass energy of the colliding partons: 
s = XiXjipA +PbY- In symmetric colliders, where pa = Pb = P, we have s = AxiXjp^ — XiXjs. 
The sum in (fTS]) runs over all pairs of light partons {i,j) contributing to the process. The 
factorization scheme serves as a method to systematically eliminate collinear divergencies from the 
parton cross section aij and absorb them into the parton distribution functions. A detailed theoretical 
justification for the applicability of the factorization ansatz to heavy quark production can be found 
in reference |73] . According to this analysis effects such as an intrinsic heavy quark component in the 
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Figure 5. Parton distribution functions (PDFs) of u quarks, u quarks, d quarks, d quarks and 
gluons inside the proton. The parametrization is CTEQ3M |75| . The scale at which the PDFs are 
evaluated was chosen to be = 175 GeV {p? = 30625 GeV^). 



hadron wave function and flavour excitation processes like qt — > qt do not lead to a breakdown of the 
conventional factorization formula. 

3.1.2. Parametrizations of parton distribution functions The PDFs are extracted from measurements 
in deep inelastic scattering experiments where either electrons, positrons or neutrinos collide with 
nucleons. A variety of experiments of this kind have been conducted since the late 1960s and provided 
hard evidence for the parton (quark) model in the first place. During the 1990s the experiments 
ZEUS and HI [71] at the ep storage ring HERA at DESY have reached an outstanding precision in 
measurements of the proton structure, which meant a big leap forward for the prediction of cross 
sections at hadron colliders. 

Several parametrizations of proton PDFs have been extracted from the experimental data 
by different groups of physicists. As an example figure \5\ shows PDFs of the CTEQ3M 
parametrization [7S]. Plotted are the distributions most relevant for pp and pp collisions, the PDFs 
for u, u, d, d and the gluons. For antiprotons the distributions in figure[5]have to be reversed between 
quark and antiquark. The scale fj, was chosen to be /x = 175 GeV. One sees that the gluons start to 
dominate in the x region below 0.15. The tt production cross section at the Tevatron is dominated by 
the large x region, since the top quark mass is relatively large compared to the Tevatron beam energy 
{Mtop/^/s = 0.0875). At the LHC {Mtop/^/s = 0.0125) the lower x region becomes more important. 

There are differences between the various sets of PDF parametrizations. To give the reader a quan- 
titative impression we compare the CTEQ3M, MRSR2 and MRSD parametrizations in figure [S] As 
a representative example we plot the ratios /«(a;)MRSR2//«(a::)cTEQ3M and /«(a::)MRSD//«(a::)cTEQ3M 
for up-quarks as weU as /g(a;)MRSR2//g(a;)cTEQ3M and /g(a;)MRSD//g(a;)cTEQ3M for gluons. These 
three PDF parametrizations were chosen because they are used for predictions of the tt cross section 
at the Tevatron ^761 Ell [78| . Typical differences are on the order of 5 to 10%. 

3.1.3. The parton cross section The cross section a of the hard, i.e. short distance, parton-parton 
process ij — > tt can be calculated in perturbative QCD. The leading order processes, contributing 
with to the perturbation series, are quark-antiquark annihilation, qq — s- tt, and gluon fusion, 
gg — > ti. The corresponding Feynman diagrams for these processes are depicted in figure [71 
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Figure 6. Comparison of the CTEQ3M, MRSR2 and MRSD parametrizations for the PDF of the 
up quark and the gluon. The left plot shows the ratios /u(a:^)MR,SR2//u{2;)cTEQ3M (Ml line) and 
/u{a;)MRSD//i/(a;)cTEQ3M (dashed line). On the right hand side we plot /9(a;)MRSR2//s(a:)cTEQ3M 
(full line) and /9(a;)MRSD//9(2;)cTEQ3M (dashed line). The factorization scale is set to fj? = 
30625 GeV^. 
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Figure 7. Feynman diagrams of the leading order processes for tt production: quark-antiquark 
annihilation (qq tt) and gluon fusion {gg — > tt). 



The leading order (Born) cross sections for heavy quark production were calculated in the late 
1970s [751 HOI mi [Ml HSl IHlj) niost of them having charm production as a concrete application in 
mind. The differential cross section for quark-antiquark annihilation is given by 



do- Anal 



[(m^ - if + (m^ - uf + 2to^ 



(14) 



where s, t and u are the Lorentz-invariant Mandelstam variables of the process. They are defined by 
s = ipq +Pqf, t = {pq — Ptf and u = {pq — p if with Pi being the corresponding momentum 4-vector 
of the quark i. m denotes the top quark mass. The differential cross section for the gluon-gluon fusion 
process is given by: 

Tial r6(m^ - £)(m^ — tt) m^(s — 4m^) 



da , 

-TJ [9192 tt) 

at 



+ 



3(m2 -t)(m2 -m) 

4 (m^ - t)(m^ - u) - 2m^(m^ +f) 4 (m^ - t)(m^ - it) - 2m^(m^ + m) 



(m2 - £)2 



+ 



g {m^ — i){m^ — u) — rn^{u — i) ^ 



3 (m? — u)^ 

{m? - t) (m^ ~u) — m^{t — u) 



ty 



(15) 

The invariant variables are in this case s = {pgi + Pg2f, t = [Pgi — Ptf and u = {pgi — pif ■ The 
cross sections in (|14p and (jlSp are quoted in the form given in reference |11] . The invariants t and u 
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Figure 8. Example Feynman diagrams of next-to-leading order (NLO) corrections to quark- 
antiquark annihilation and gluon fusion, a) and c) show virtual corrections, b) and d) gluon 
bremsstrahlung graphs. 

may be expressed in terms of the cosine of the scattering angle 9 in the parton-parton centre-of-mass 
system: 



cose=\jl--l^ f=-|(l-cos^) w = -^(1 +cos^) (16) 

where is the common transverse momentum of the outgoing top quarks. 

A full calculation of next-to- leading order (NLO) corrections contributing in order to the 
inclusive parton-parton cross section for heavy quark pair production was performed independently 
by two groups: Nason et al. in 1988 [SF and Beenakker et al. in 1991 [86l [87], yielding consistent 
results. The NLO calculations involve virtual contributions to the leading order processes, gluon 
bremsstrahlung processes {qq tt + g and gg tt + g) as well as processes like g + q{q) ti+ q{q)- 
Examples of Feynman diagrams of NLO processes are given in figure [S] a) and c) display two virtual 
graphs, b) and d) two gluon bremsstrahlung graphs, corrections arise when interfering those graphs 
with the leading order graphs of figure [T] For the NLO calculation of the hadron-hadron cross section 
a{AB tt) to be consistent one has to use next-to-leading order determinations of the coupling 
constant as and the PDFs. All quantities have to be consistently defined in the same renormalization 
scheme because different approaches distribute the radiative corrections differently among the parton- 
parton cross section, the PDFs and as- Most authors use the MS or an extension of the MS scheme. 
First NLO cross section predictions for production at the Tevatron (y^ = 1.8 TeV) yielded values 

of about 4 pb [siiHziiiiEnj 



At energies close to the kinematic threshold, s = 4M^Qp, the quark-antiquark annihilation process 
is the dominant one, if the incoming quarks are valence quarks, as is the case of pp collisions. At 
the Tevatron 80 to 90% of the tt cross section is due to quark-antiquark annihilation [751 EZl [UT| . 
At higher energies the giuon-gluon fusion process dominates for both pp and pp collisions. That is 
why one can built the LHC as a pp machine without compromising the parton-parton cross section. 
Technically, it is of course much easier to operate a pp collider, since one spares the major challenge 
to produce high antiproton currents in a storage ring. For the Tevatron the ratio of NLO over LO 
cross sections for gluon-gluon fusion is predicted to be 1.8 at Mtop = 175 GeV/c^, for quark-antiquark 
annihilation the value is only about 1.2 Since the annihilation process is dominating, the overall 
NLO enhancement is about 1.25. 



3.1.4- Soft gluon resummation Contributions to the total cross section due to radiative corrections 
are large in the region near threshold (s — iM^^p) and at high energies (s > 400 M^^^). Near 
threshold the cross section is enhanced due to initial state gluon bremsstrahlung (ISGB) [92j. This 
effect is important for tt production at the Tevatron, but not for the LHC where gluon splitting and 
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Table 3. Cross section predictions for production at next-to-leading order (NLO) in perturbation 
theory including the resummation of initial state gluon bremsstrahlung (ISGB). The predictions are 
given for a top mass of Mtop = 175 GeV/c^. The cross sections are given for pp collisions at the 
Tevatron = 1.8 TeV and = 1.96 TeV) and pp collisions at the LHC (v^ = 14 TeV). The 
factorization and renormalization scale is set to = Mtop to derive the central values. It has to 
be stressed that the authors use different sets of PDF parametrizations. Part of the differences can 
be attributed to this fact. In reference 1781 Laenen et al. do not provide a direct prediction for 
Mtop = 175 GeV/c^. To compare their prediction with those of other authors I choose to linearly 
interpolate the given values for A/top = 174 GeV/c^ and Mtop = 176 GeV/c^, thereby neglecting the 
functional form of the mass dependence. Within the assigned errors this simplification is acceptable. 
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flavour excitation are increasingly important effects. The calculation at fixed next-to-leading order 
(a^) perturbation theory has been refined to systematically incorporate higher order corrections due 
to soft gluon radiation. Technically, this is done by applying an integral transform (Mellin transform) 
to the cross section: 

aN{tt)= f dp p^-^ <j{p;ti) (17) 
Jo 

where p = AM^^^/s is a dimensionless parameter. In Mellin moment space the corrections due to 
soft gluon radiation are given by a power series of logarithms InA^. For p. = Mtop the corrections 
are positive at all orders. Therefore, the resummation of the soft gluon logarithms yields an increase 
of the tt cross section with respect to the NLO value. Four different groups have presented cross 
section predictions based on the resummation of soft gluon contributions: (1) Laenen et al. [MlllH], (2) 
Berger and Contopanagos [SI [SI [75], (3) Bonciani et al. (BCMN) (SB [SI [771 [SI] and (4) Kidonakis 
et al. [SZllSHllSSI- Their predictions for Mtop = 175 GeV/c^ are summarised in table [31 

In the region very close to the kinematic threshold non-perturbative effects become dominant and 
the perturbative approach breaks down. The resummation can therefore not sensibly be extended into 
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this region. That is why Laenen et al. introduce a new scale /io with Aqcd <C /^q ^ -^^top where 
they stop the resummation. The concrete choice of scale is to some extent arbitrary. The uncertainty 
quoted by Laenen et al. is derived by varying the scale /ip. To calculate the central value of the 
cross section they use = 0.1 Mtop for gg-annihilation and /io — 0.25 Mtop for gg-fusion. The values 
are not required to be the same, since the respective perturbation series have different convergence 
properties. 

Berger and Contopanagos derive the infra-red cut-off within their calculation and thereby 
define a perturbative region where resummation can be applied. Since /io is derived, it is not treated 
as a source of error in this approach. The theoretical uncertainty quoted by Berger and Contopanagos 
is derived by varying the factorization and renormalization scale ^ between 0.5 Mtop and 2 A/top- The 
central value, given in tabled is calculated using the CTEQ3M PDFs. The uncertainty due to the 
choice of the PDF parametrization is about 4%. Resummation effects are of appreciable size: The 
resummed total cross sections (for -^s = 1.8 TeV and -/s — 2.0 TeV, respectively) are about 9% above 
the NLO cross sections. Berger and Contopanagos predict an increase of 37% in cross section, when 
going from ^/s = 1.8 TeV to ^/s = 2.0 TeV in Run II of the Tevatron. The cross section value for the 
LHC, a = 760 pb, merely reflects an estimate and is not accompanied by an uncertainty. Due to the 
much larger centre-of-mass energy at the LHC the near threshold region is much less important for ti 
production, reducing the significance of ISGB and the need for resummation of these contributions. 

While the first two groups have only resummed leading logarithmic terms (LL), BCMN also 
include next-to- leading logarithms (NLL). They used a different resummation prescription which does 
not demand the introduction of an additional infra-red cut-off /io [95]. The LL result of BCMN 
shows only an increase of about 1% compared to the NLO cross section at the Tevatron [96]. When 
taking NLL terms into account the increase is 4% [77] . In table [3| we quote the NLL results as 
updated in reference [91] with the newest set of PDFs. The errors quoted by BCMN include PDF 
uncertainties, which are evaluated by using sets of PDF parametrizations that provide an estimate of 
"l-cr" uncertainties. In the case of CTEQ [lOOl [101] 40 different sets are available, for MRST [l02] 
there are 30 sets. Kidonakis et al. resum leading and subleading logarithms up to order in an 
attempt to reduce the dependence of the cross-section on the renormalization scale /i compared to 
the NLO calculation. The uncertainty quoted by Kidonakis et al. is dominated by the choice of the 
kinematic description of the scattering process, either in one-particle-inclusive or pair- invariant-mass 
kinematics. The Tevatron predictions given in table [3| are the average of the two choices. The LHC 
prediction is based on the one-particle-inclusive value, since this is believed to be more appropriate 
when the cross section is dominated by gluon-gluon fusion. 

The cross section predictions are strongly dependent on the top quark mass, which is illustrated 
in figure [D The plot shows the predictions for the Tevatron at ^/s — 1.8 TeV, they are in good 
agreement within the given errors. For comparison the Run I measurements of CDF [1031 1104] and 
D0 |105[ 1106] are shown. Within the large errors the measurements agree well with the theoretical 
predictions. It is obvious that a precise cross section measurement has to be accompanied by a precise 
measurement of the top quark mass to provide a basis for a stringent test of the theory. 

3.2. Single top quark production 

Top quarks can be produced singly via electroweak interactions involving the Wth vertex. There are 
three production modes which are distinguished by the virtuality of the W boson {Q^ = — g^, 
where q is the four-momentum of the W): 

(i) the t-channel {q'^ — t ): A virtual W strikes a b quark (a sea quark) inside the proton. The 
W boson is spacelike {q^ < 0). This mode is also known as W-gluon fusion, since the b quark 
originates from a gluon splitting into a bb pair. Feynman diagrams representing this process are 
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Figure 9. The tt cross section at y/s = 1.8 TeV as a function of the top quark mass. The thick 
full line represents the central value predicted by Catani et al. |96II77I . The thick dashed line shows 
the prediction by Berger and Contopanagos 1941 . The thick dashed-dotted line is the prediction by 
Kidonakis et al. |99| . The thin lines indicate the upper and lower uncertainties of the predictions. 
The Run I measurements of CDF |103l Il04l and D0 11051 11061 are shown for comparison. 



shown in figure llOb and figure llOb . ly-gluon fusion is the dominant production mode, both at 
the Tevatron and at the LHC, as wiU be shown in the discussion below. 

(ii) the s-channel [q^ — s): This production mode is of Drell-Yan type. A timehke W boson with 

> (Mtop + rrii,)'^ is produced by the fusion of two quarks belonging to an SU(2) isospin doublet. 
See figure [TUb for the Feynman diagram. 

(iii) associated production: The top quark is produced in association with a real (or close to real) 
W boson (g2 = M^). The initial b quark is a sea quark inside the proton. Figure llOd shows 
the Feynman diagram. The cross section is negligible at the Tevatron, but of considerable size at 
LHC energies where associated production even supercedes the s-channel. 

In pp and pp collisions the cross section is dominated by contributions from up and down quarks 
coupling to the W boson on one hand side of the Feynman diagrams. That is why the (u, d ) quark 
doublet is shown in the graphs of figurc[TOl There is of course also a small contribution from the second 
weak isospin quark doublet, (c, s); an effect of about 2% for s- and 6% for i-channel production [107] . 
Furthermore, we will only consider single top quark production via a Wtb vertex. The production 
channels involving a Wtd or a Wts vertex are strongly suppressed due to small CKM matrix elements: 
0.0048 < \Vtd\ < 0.014 and 0.037 < \Vts\ < 0.043 [TD]. Thus, their contribution to the total cross 
section is quite small: ~ 0.1% and ~ 1%, respectively |108j . In the following paragraphs we will review 
the theoretical cross section predictions for the three single top processes and the methods with which 
they are obtained. 

3.2.1. W-gluon fusion The VF-gluon fusion process was already suggested as a potentially interesting 
source of top quarks in the mid 1980s [1091 1110] and early 1990s |111| . If the h quark is taken to be 
massless, a singularity arises when computing the diagram in figure llOb in case the final h quark 
is coUinear with the incoming gluon. In reality the non-zero mass of the b quark regulates this 
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Figure 10. Representative Feynman diagrams for the three single top production modes, a) and 
b) show Vl^-gluon fusion graphs, c) the s-channel process and d) associated production. We chose 
to draw the graphs in a), b) and c) with the (ii, d) weak-isospin doublet coupling to the W . This is 
by far the dominating contribution. In general, also the (c, s) doublet contributes. The graphs show 
single top quark production, the diagrams for single antitop quark production can be obtained by 
interchanging quarks and antiquarks. 

(a) (b) (c) 

u (rf) d [u] u (d) d (m) u (d) d («) 
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Figure 11. a) shows the leading order Feynman diagram for W gluon fusion if the perturbative 
calculation is performed involving a b quark distribution function, b) and c) show two examples of 
as corrections to the leading order diagram, b) is a correction to the b vertex, c) a correction to the 
light quark vertex. 



coUinear divergence. When calculating the total cross section the coUinear singularity manifests itself 
as terms proportional to ln(((5^ + M^^p)/ml) with = —q^ being the virtuality of the W boson. 
These logarithmic terms cause the perturbation series to converge rather slowly. This difficulty can 
be obviated by introducing a parton distribution function (PDF) for the b quark, 6(a;,/z^), which 
effectively resums the logarithms to all orders of perturbation theory and implicitly describes the 
splitting of gluons into bb pairs inside the colliding hadrons |112[ I113j . Once a b quark distribution 
function is introduced into the calculation, the leading order process is qi + b ^ qj + t as shown in 
figurctTTk. In this formalism the process shown in figure [TUk is a higher order correction which is 
already partially included in the b-quark distribution function. The remaining contribution is of order 
l/\n{{Q^ + AftQp)/m^) with respect to the leading order process in figure [TTk. Additionally, there are 
also corrections of order a^: Two examples of those are shown in figure [TTb and figure [Tib. The leading 
order differential cross section calculated from the Feynman graph in figure [TTk for quark-quark or 
antiquark-antiquark collisions is given by |109j : 



ij Tra^ 2 |2 top f-io\ 



For quark- antiquark collisions the result is: 
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Table 4. Predicted total cross sections for single top quark production processes. The cross sections 
are given for pp collisions at the Tevatron {y/s = 1.8 TeV or ^/s = 1.96 TeV) and pp collisions at the 
LHC {y/s = 14 TeV). The cross sections of the t- and s-channel process are taken from reference I116| 
and were evaluated with CTEQ5M1 PDFs. The uncertainties were evaluated in reference I117| . The 
values for associated production are taken from reference 11181 (Tevatron at y/s = 1.96 TeV and 
LHC) and reference. 11071 (Tevatron at ^/s = 1.8 TeV). All cross sections are given for a top quark 
mass of A/top = 175 GeV/c^. 



Process ^/s a{t — channel) a{s — channel) cr{Wt) 



pp^t/t 1.80 TeV 1.45t° j^pb (0.75 ± 0.10) pb 0.141° °^ pb 
pp^t/t 1.96 TeV 1.98tjjJ^pb (0.88 ± 0.11) pb 0.094t°:J^}^ pb 



pp t 14.0 TeV (156 ± 8) pb (6.6 ± 0.6) pb W.Otj.g pb 
pp~yi 14.0 TeV (91 ± 5) pb (4.1 ± 0.4) pb 14.0l^ j pb 



oi-a, — 5^/(4 tt) is the weak fine structure constant. The Mandelstam variable t is given by 

f.-f(l-!^)(l-^).(l-cos^) (20) 

and s + t + u = ml + M'^^^. Another leading order calculations of the cross section was done van der 
Heide et al. [mj. 

The NLO calculation at first order in comprises the square of the Born terms, (fT8)) and p9| . 
plus the interference with the virtual graphs plus the square of the real graphs with one single QCD 
coupling, e.g. figure ITTh . Bordes and van Eijk presented an NLO calculation based on the formalism 
described above in 1995 jll2| . They predict an enhancement of +28% of the NLO cross section over 
the Born cross section for the Tevatron operating at =1.8 TeV. Bordes and van Eijk used small 
masses for gluons and quarks to regularize infrared and coUinear divergencies. Mass factorization 
was performed in the deep inelastic scattering (DIS) scheme. Stelzer et al. |113[ 1115] performed an 
NLO calculation entirely based on the MS factorization scheme. They predict a decrease of the cross 
section when going from leading order to NLO by about -8% to -10%. 

The latest calculation was done by Harris et al. in 2002 116] and contains full differential 
information, such that experimental acceptance cuts and jet definitions can be applied. Their results 
are summarised in table |4l We quote only these latest results for the cross sections, since previous 
calculations used different PDFs, which by itself leads to big differences in the predicted values. Harris 
et al. compare their results with those given by Stelzer et al. using the latest PDFs. The agreement is 
very good, within 1%. The cross sections given for the Tevatron are the sum of top and antitop quark 
production. In pp collisions at the LHC the ly-gluon fusion cross section differs for top and antitop 
quark production, which are therefore treated separately in table [H The increase in the centre-of-mass 
energy from 1.80 TeV to 1.96 TeV in Run II of the Tevatron is predicted to yield a 33% increase in 
the total cross section. The ratio of awg/o'tt is about 30%, for the Tevatron as well as for the LHC. 
The uncertainties quoted in table |4] are evaluated in reference [117j and include the uncertainties due 
to the factorization scale /i, the choice of PDF parameterization, and the uncertainty in the top quark 
mass. The factorization scale uncertainty is ±4% at the Tevatron and ±3% at the LHC. The central 
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value was calculated with fi^ — + M^^^ for the b quark PDF. The scale for the light quark PDFs 
was set to /i^ = Q^. 

Of course the i-channel single top cross section depends on the top quark mass. The current 
uncertainty in the top quark mass (Am = 4.3GeV/c^) corresponds to about 7% uncertainty in the 
cross section at the Tevatron and 3% at the LHC. The dependence is approximately linear in the 
relevant mass range. 

At first glance it is astonishing that the cross section for M^-gluon fusion is of the same order of 
magnitude as tt production although it is a weak interaction process. There are several issues that 
lead to this relative enhancement [109] : 

(i) The parton cross section of the W-gluon fusion mode scales like as opposed to the tt 
cross section which, as a typical strong interaction process, scales like 1/s. At the Tevatron 
the subprocess energies are not much greater than Mw, so one does not gain from the scaling 
behaviour. However, at the LHC the effect is present. 

(ii) Single top production is kinematically enhanced compared to tt production, since only one heavy 
top quark is produced. For single top quark production the parton distribution functions are 
therefore typically evaluated at half of the value of x needed for the strong process. Since the 
PDFs are monotonically decreasing functions, see figure[5l single top quark production is relatively 
enhanced to ti pair production. 

(iii) The H^-gluon fusion process is enhanced by logarithmic terms originating from the coUinear 
singularity discussed above. 

In general, VK-gluon fusion events have three quark jets originating from the hard interaction: 
(1) the b quark jet from the top quark decay, (2) the light quark jet, and (3) the &-quark jet which 
comes from the initial gluon splitting. The transverse momentum distribution of these jets is shown 
in figure [T2k . The b quark is most of the times the hardest jet, the peak of the distribution is around 
60 GeV. The light quark pt distribution peaks around 25 GeV, but has a long tail to high values. The 
b quark pt distribution peaks at low values. A large share of the b jets will therefore not be identified, 
since an experimental analysis will require some lower pt cut off, typically at 15 GeV. 

An interesting feature of single top quark production (in the t-channel and in the s-channel) 
is that in its rest frame the top quark is 100% polarized along the direction of the d quark (d 
quark) |1191 1120|, 11071 1115j . The reason for this is that the W boson couples only to fermions with 
left-handed chirality. Consequently, the ideal basis to study the top quark spin is the one which uses 
the direction of the d quark as the spin axis [120] . In pp collisions at the Tevatron M^-gluon fusion 
proceeds via ug dtb in 77% of the cases. The d quark can then be measured by the light quark 
jet in the event. The top quark spin is best analyzed in the spectator basis for which the spin axis 
is defined along the light quark jet direction. However, 23% of the events proceed via dg — *■ utb, in 
which case the d-quark is moving along one of the beam directions. For these events the spectator 
basis is not ideal, but since the light quark jet occurs typically at high rapidity the dilution is small. 
In total, the top quark has a net spin polarization of 96% along the direction of the light quark jet 
in i-channel single top quark production [120j . In s-channel events the best choice is the antiproton 
beam direction as spin basis. In 98% of the cases the top quark spin is aligned in the antiproton 
direction |120j . 

At the Tevatron top quarks are not produced as ultrarelativistic particles. Therefore, the chirality 
eigenstates are not identical to the helicity eigenstates. The spin asymmetry A = (A^| — A^j)/ (N^+Ni) 
is 0.91 for the i-channel and 0.96 for the s-channel. 

Since the top quark does not hadronizc, its decay products carry information about the top quark 
polarization. A suitable variable to investigate the top quark polarization is the angular distribution 
of electrons and muons originating from the decay chain t + b, ^ £^ + vg. If 6qi is the 

angle between the charged lepton momentum and the light quark jet axis in the top quark rest frame, 




Figure 12. (a) Transverse momentum (pt) distribution of the quarlt jets from Vl^-gluon fusion 
events: (1) the b quarls from the top quarlt decay (solid line), (2) the hght quark jet (dotted line), 
and (3) the b quark from the gluon splitting. The distributions are calculated for the Tevatron 
at y'S = 2TeV. (b) Angular distribution of the charged lepton in W-gluon fusion events at the 
Tevatron (y/s = 2TeV). 9 is the angle between the lepton momentum and the light quark jet axis. 
For comparison the same distribution is shown for the sum of all background processes {W + jj and 
ti). Both plots, a) and b), are taken from reference |115) . 

(a) (b) (c) (d) 
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Figure 13. Examples of -corrections to the s-channel single top production mode. 

the angular distribution is given by 0.5 (1 + COS 6qi). A theoretical prediction for this quantity is shown 
in figure [T^ [115] . Single top quark events show a distinct slope which differs significantly from the 
nearly flat background. 



3.2.2. s-channel production The s-channel production mode of single top quarks probes a 
complementary aspect of the weak charged current interaction of the top quark, since it is mediated 
by a timelike W boson with > (Mtop + nib)^ as opposed to a spacelike W boson in the t-channel 
process. The leading order s-channel process is depicted in figure [TOb . Feynman diagrams yielding 
corrections of order are shown in figure If 31 The first order corrections depicted in figure 113b and 
figure [T3li have the same initial and final states as the W-gluon fusion diagrams shown in figure [TOk 
and figure [TOb. However, these two classes of diagrams do not interfere because they have a different 
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colour structure. The tb pair in the W^-gluon fusion process is in a colour-octet state, since it originates 
from a gluon. In s-channel production the tb pair forms a colour-singlet because it comes from a W . 
The different colour structure implies that both groups of processes must be separately gauge invariant 
and, therefore, they cannot interfere [1091 11111 1107j . 

Several groups have calculated the cross section for the s-channel production mode in leading or 
next-to-leading order, respectively [1211 11221 1107|. IllGj . The leading order result for the partonic 
cross section is given by [121| 



- ^ ■ \V^,? m ■ . ■ I go - ' r: I With qo^ - " . (21) 



,2„,2, - ^^'°P ( Ml^qi\ ^^^^^^ _s + Ml^ 

[s-M^f 3^/l 



In table m we quote the latest NLO results by Harris et al. |116j which are in very good agreement 
with the earlier calculations by Smith/ Willenbrock [122] and Mrenna/Yuan |123| . The later authors 
have resummed soft gluon emission terms in the s-channel cross section. Their resummed cross 
section is about 3% above the NLO value. The predictions for the s-channel cross section have a 
smaller uncertainty from the PDFs than for the i-channel because they do not depend as strongly on 
gluon distribution functions as the t-channel calculation does. The uncertainty in the top quark mass 
leads to an uncertainty in the cross section of about 10% (Am = 4.3GeV/c^, ^/s = 1.96 TeV)) ^117| . 

The ratio of cross sections for the t-channel and s-channel mode is 2.3 at the Tevatron 
{y/s = 1.96 TeV) and 23 at the LHC. In pp collisions the gluon initiated processes, tt production 
and VF-gluon fusion, dominate by far over s-channel single top quark production which is a quark- 
antiquark annihilation process. The s-channel signal will therefore most likely be obscured at the 
LHC. Thus, it will be essential to observe s-channel single top quark production at the Tevatron. 



3.2.3. Associated production The third single top quark production mode is characterised by the 
associated production of a top quark and an on-shell (or close to on-shell) W boson. Studying 
associated production is interesting, since it probes a different kinematic region of the Wtb interaction 
vertex than s-channel or t-channel production and thereby provides complementary information. 
Predictions for the cross sections are given in table H) The calculation by T. Tait [118j includes 
higher order corrections proportional to 1 / \n{M'^^^ / ml) . The prescription is very similar to the one 
described in section 13.2.11 for the t-channel mode. However, it is not a full NLO calculation. Tait 
provides predictions for the Tevatron {^/s = 2.00 TeV) and the LHC [s/s = 14 TeV). It is obvious 
from these numbers that associated production is negligible at the Tevatron, but is quite important at 
the LHC where it even exceeds the s-channel production rate. The errors quoted in tabled] include the 
uncertainty due to the choice of the factorization scale (±15% at the LHC) and the parton distribution 
functions (±8% at the LHC). The uncertainty in the top mass (Am — 5GeV/c^) causes a spread of 
the cross section by ±9% at the LHC. A second calculation by Belyaev and Boos yields a much higher 
cross section for associated production: a{Wt/t) = 62.0l3^g^pb [124] . This result was obtained using 
the CompHEP program [125j . First studies [118] show that associated single top quark production 
will be observable at the LHC with data corresponding to about 1 fb^^ of integrated luminosity. 



3.3. Top quark decay 

In the SM top quarks decay predominantly into a b quark and a W boson. The decays t d + 
and t ^ s -I- are CKM suppressed relatively to t — > 6 -I- by factors of |Vtdp and | Vts p. If we 
assume the CKM matrix to be unitary the values of these matrix elements can be inferred from other 
measured matrix elements: 0.0048 < \Vtd\ < 0.014 and 0.037 < \Vts\ < 0.043 [10]. In the discussion 
of the following paragraphs we will therefore only consider the decay t ^ b + W^. Potential non-SM 
decays which would signal new physics will be discussed in section [9l 
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At Born level the amplitude of the decay f — > 6 + is given by 

M[t^b + W)^'^h e(V7M t ■ (22) 

The decay amplitude is dominated by the contribution from longitudinal W bosons because the decay 
rate of the longitudinal component scales with M^^^. In contrast, the top quark decay rate into 
transverse W bosons increases only linearly with Mtop- In both cases the couples solely to 6 
quarks of left-handed chirality (a general feature of the SM). Since the b quark is effectively massless, 
compared to the mass scale set by Mtop, left-handed chirality translates into left-handed helicity for 
the b quark. If the b quark is emitted anti-parallel to the top quark spin axis, the must be 
longitudinally polarized, = 0, to conserve angular momentum. If the b quark is emitted parallel 
to the top quark spin axis, the W'^ boson has helicity = — 1 and is transversely polarized. Thus, 
elementary angular momentum conservation forbids the production of W bosons with positive helicity, 
= +1, in top quark decays. The ratios of decay rates into the three W helicity states are given 
by i: 



Aih"^ = -1) : Aih"^ = 0) : Aih"^ = +1) = 1 : — ^ : 0. (23) 

w 

Strong next-to-leading order corrections to the decay rate ratios have been calculated, lowering the 
fraction of longitudinal W bosons by 1.1% and increasing the fraction of left-handed W bosons by 
2.2% |126[ 1127] . Electroweak and finite widths effects have even smaller effects on the helicity ratios, 
inducing corrections at the per-mille level [128j . For the decay of antitop quarks negative helicity is 
forbidden. In the SM the top quark decay rate, including first order QCD corrections, is given by 



r, = — — ^ \Vtb\^ (i-^^i (1 



8^^/2 ''^^ MlJ V Ml 



37r 



(24) 



top / \ top ^ 

with y = {Mw/Mtopf and f{y) = liP- ji - 2.5 - 3y 4.5y2 - iy^Xxiy [Uni [1301 ■ Using 
y = (80.45/174.3)^ we find /(y) = 3.85. The QCD corrections of order as lower the Born decay rate 
by -10%. A useful approximation of (EH) is given by Ttop =i 175 MeV/c^ • [Mto^jMwf [HI [HI]- The 
decay width increases from 1.07 GeV/c^ at Mtop = 160 GeV/c^ to 1.53 GeV/c^ at Mtop = 180 GeV/c^. 
Expression ((24| neglects higher order terms proportional to m^/Mjop and o?^. Corrections of order 
were lately calculated, they lower Ftop by about -2% [132[ 1133] . Because the top quark width is small 
compared to its mass, interference between QCD corrections to production and decay amplitudes has 
a small effect of order ©(agFtop/Mtop) .134j. The decay width for events with hard gluon radiation 
[Eg > 20 GeV) in the final state has been estimated to be 5 - 10% of Ftop, depending on the gluon jet 
definition (cone size AR = 0.5 to 1.0) |135| . Electroweak corrections to Ftop have also been calculated 
and increase the decay width by 6ew — +1.7% |1361 1137j . Taking the finite width of the W boson 
into account leads to a negative correction 5t = —1.5% such that Sew and 5t almost cancel each 
other 

The large top decay rate implies a very short lifetime of Ttop = 1/Ftop ~ 4 • 10^^^ s which is 
smaller than the characteristic formation time of hadrons Tform ~ 1/Aqcd ~ 2 • 10"^'' s. In other 
words top quarks decay before they can couple hadronically to light quarks and form hadrons. The 
lifetime of tt bound states, toponium, is too small, F^f ~ 2 Ftop, to allow for a proper definition of a 
bound state with sharp binding energy. This feature of a heavy top quark was already pointed out in 
the early and mid 1980s .139, Hi [HI]. 

Even though top hadrons cannot be formed, there are other long-distance QCD effects associated 
with hadronization which have to be considered. The colour structure of the hard interaction process 
influences the subsequent fragmentation and hadronization process. In the process e'^e^ — > ti the 
top and antitop quark are produced in a colour-singlet state. In hadronic collisions, on the contrary, 
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the production cross section is dominated by configurations where the t or t forms a colour-singlet 
with the proton or antiproton remnant, respectively. The colour field - or in the picture of string 
fragmentation - the string carries the more energy the further the top quark and the remnant are 
apart. If the distance in the top-remnant centre-of-mass system reaches about 1 fm before the top 
quark decays, the colour string carries enough energy to form light hadrons. Whether or whether 
not a significant fraction of top events exhibit the described "early" fragmentation process, depends 
strongly on the centre-of-mass energy of the hadron collider. While at Tevatron energies early top 
quark fragmentation effects are negligible |141j . they may well play a role at the LHC, where top 
quarks are produced with a large Lorentz boost. If there is no early fragmentation, long-range QCD 
effects connect the top quark decay products, the b quarks or the quarks from hadronic W decays. 
Even if early fragmentation happens, the fragmentation of heavy quarks is hard, as seen in c and b 
quark decays, i.e. the fractional energy loss of top quarks as they hadronize is small. Therefore, it 
will be quite challenging to observe top quark fragmentation experimentally, even at the LHC. 

Within the constraints discussed above we can assume that top quarks are produced and decay 
like free quarks. The angular distribution of their decay products follow spin ^ predictions. The 
angular distribution of W bosons from top decays is propagated to its decay products. In case of 
leptonic W decays the polarization is preserved and can be measured [142j . The angular distribution 
of charged leptons from W decays originating from top quarks is given by 

1 dr _ 3 sin' ee + M^{1- cos 9i)^ 

rdcose, ~4 Ml^ + 2M^ ^ ' 

where tt — 6'^ is the angle between the b quark direction and the charged lepton in the W boson rest 
frame [l43j . 

4. Experimental techniques 

Advanced experimental techniques are needed to detect and reconstruct top quark events in hadronic 
collisions. Large scale general-purpose detectors are employed for that task, their overall structure is 
quite similar. Early searches for the top quark were conducted with the UAl and UA2 experiments at 
the CERN SppS. The detectors CDF and D0 are currently in operation at the Fermilab Tevatron and 
we will discuss those as typical examples of collider detectors in more detail in sections 14.21 and 14.31 
respectively. In the future the LHC will also feature two general-purpose experiments, CMS and 
ATLAS, which are currently under construction at CERN. 

Usually, general purpose collider detectors feature rotational symmetry with respect to the 
nominal beam axis and forward-backward symmetry with respect to the nominal interaction point in 
the centre of the detector. Therefore, a right-handed coordinate system is chosen such that the origin 
is at the centre of the detector and the z-axis points along the symmetry axis parallel to the beam 
(in case of the Tevatron the proton beam). The x- and y-axes define the transverse plane, the y-axis 
points vertically upwards. It is often practical to replace the x and y coordinates by the azimuth angle 
(j) = arctan(2//x), measured in the transverse plane with respect to the cc-axis, and the polar angle 9, 
measured with respect to the z-axis (6 — 0). <j) takes values from to 27r, 9 from to tt. Instead of 9 it 
is often handy to use the pseudorapidity 77, which is defined asr] = — ln(tan0/2). For massless particles 
ry is equal to the rapidity y ~ 1/2 \n.[{E + pz) / (E — pz)) , which is an important quantity because the 
rapidity difference Aj/ between two particles is Lorentz-invariant. It is also common to calculate the 
angle between two particles in terms of the distance in the rj-cj) plane: Ai? = At^^ + /S.(p. Table [5] 
gives a summary on the definition of kinematic variables. 
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4-1. The Tevatron Collider 

As most of the knowledge about the top quark was obtained from measurements of pp collisions 
at the Tevatron we briefly discuss this accelerator here. The Tevatron is located at the Fermi 
National Accelerator Laboratory (Fcrmilab) near Chicago. In order to reach energies of 980 GeV 
per beam, a system of several accelerators is needed. In the first stage of acceleration, a Cockcroft- 
Walton pre-accelerator is used to generate negatively charged hydrogen ions out of hydrogen gas and 
then accelerate them via electric fields up to an energy of 750 keV. Afterwards, the ions enter an 
approximately 150 m long linear accelerator, where they are accelerated up to 400 MeV by oscillating 
electric fields. Before leaving this acceleration stage, the ions pass through a carbon foil, which removes 
their negative charges (electrons). As a result one obtains a beam of protons that is subsequently 
bent in a circular path by the magnets of a circular accelerator, called the booster. On its way out of 
the booster, the beam has an energy of 8 GeV. In the next stage, the protons enter the Main Injector, 
a multitask accelerator completed in 1999. This machine accelerates protons up to 150 GeV. Some 
protons are accelerated to 120 GeV and used for antiproton production. They are forced to collide 
with a nickel target, which is installed at the antiproton source facility. The interactions with the 
target produce a variety of particles, among them many antiprotons, which are being collected, focused 
and finally stored in the Accumulator Ring. As soon as a suSicient number of antiprotons has been 
produced, they are sent to the Main Injector, which accelerates them up to an energy of 150 GeV. In 
the final stage, the proton and antiproton beams with 150 GeV energy are injected in the Tevatron, 
a circular accelerator with a circumference of about 6 km, which is the most powerful operational 
hadron accelerator worldwide. Each beam is accelerated to an energy of 0.98 TeV which is equal to 
(anti-)protons reaching velocities of 0.9999995 times the speed of light. These beams arc forced to 
collide with each other at two interaction regions, producing a centre of mass energy of y/s = 1.96 TeV. 
The general purpose experiments CDF and D0 are placed at these collision points. The performance 
of a collider is described with a quantity called luminosity, C. The event rate for a certain process 
with cross section a is given by the product N = a-C. For a certain time interval the number of events 
produced by this process is given hj N = a ■ J Cdt, where J Cdt is the integrated luminosity. The 
peak luminosity is reached at the begin of a store after protons and antiprotons have been injected. 
Over the period of a store the luminosity slowly decreases as collisions and beam gas interactions lead 
to a lowering of beam currents, that is a loss of protons and antiprotons stored in the Tevatron. The 
maximum value of luminosity that has been achieved until May 2005 is 13.0-10^"^cm~^s~^, which is 
about 60% above the design value of 8-10^^cm~^s~^. 
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Figure 14. Elevation view of one half of the CDF II detector. 



4-2. The Collider Detector at Fermilab (CDF) 

CDF is in operation since 1987. The first Tevatron run (Run I) lasted until 1995, when a substantial 
upgrade program of the accelerator complex and the collider detectors began. The installation of 
the renewed experiment finished in 2001 and after a commissioning period of about 1 year CDF II 
is taking physics quality data. Run II is expected to last until 2009 and accumulate an integrated 
luminosity corresponding to 4.4 - 8.5 fb~^ [144j . Although we include physics results obtained in Run 
I, we describe CDF in its upgraded form after 2001 (CDF II), see reference |145| for a brief overview. 
A detailed description of the Run I detector can be found elsewhere [1461 1147| . Figure [T3] shows an 
elevation view of one half of the CDF II detector. 

Tracking System In CDF II particle collisions take place within a luminous region which is 
approximately described by a Gaussian distribution with a width of 30 cm. The inner part of CDF 
II is dedicated to reconstruct the trajectories of charged particles and measure their momenta. The 
entire tracking volume is immersed into a solenoidal magnetic field of \B\ = 1.4 T. The transverse 
momentum px of charged particles is measured by determining the curvature of their trajectories in the 
magnetic field. The tracking system has two major components, the central drift chamber (COT) [148j 
and the silicon tracker that itself is composed of three subsystems: Layer 00, the SVX II |149) and 
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the Intermediate Silicon Layers (ISL) |150| . Layer 00 consists of one layer of radiation-tolerant single- 
sided silicon strip detectors which are directly glued onto the beam pipe. SVX II extends radially from 
ri = 2.4 cm to r2 = 10.7 cm, it is 96 cm long and provides full angular coverage in the pseudorapidity 
region of \ri\ < 2.0. The SVX II is segmented in three cylindrical barrels with beryllium bulkheads 
at each end for mechanical support, cooling the modules and facilitating the readout. The SVX II 
features five layers of double-sided silicon strip detectors which provide measurements in the r-(j> and 
r-z views. The ISL consists of one layer of double-sided strip detectors in the central region, |?7| < 1.0, 
and two layers in the forward and backward regions. The central ISL is located at a radius of 22 cm 
and allows to robustly link tracks between the SVX 11 and the central drift chamber (COT). The COT 
is a cylindrical open-cell drift chamber with inner and outer radii of 40 and 137 cm. It is designed 
to find charged particles in the pseudorapidity region of \r]\ < 1.0 with transverse momenta as low 
as 400 MeV/c. Four axial and four stereo super-layers with 12 sense wires each provide a total of 
96 measurements. The drift gas is a 50:50 Argon-Ethane admixture and the drift field is 1.9 kV/cm, 
yielding a maximum drift time of 180 ns. When combining measurements in the silicon tracker and 
the COT the momentum resolution for charged particles is Spr/Px < 0.17% GeV~^ c. 

Between the COT and the solenoid a Time-of-Flight (TOP) system has been installed to improve 
particle identification |151j . The TOF consists of scintillation panels which provide both timing and 
amplitude information. The main purpose of the TOF system is to distinguish between charged pions 
and kaons with momenta p < 1.6 GeV/c, which is important to reconstruct clean samples of exclusive 
b hadron decays and to tag the flavor of b hadrons in CP violation and oscillation measurements. 

Calorimetry The solenoid of CDF is surrounded by calorimeters, where measurements of the 
electromagnetic and hadronic showers are performed. The calorimeters are segmented in azimuth 
and in pseudorapidity to form a projective tower geometry which points back to the nominal 
interaction point. In addition, the calorimeters are segmented longitudinally to distinguish between 
electromagnetic (EM) and hadronic (HAD) showers. All CDF calorimeter systems are sampling 
calorimeters consisting of a sandwich of absorber (lead or iron) and scintillator material. The 
inner part of the calorimeters is used to measure and identify electromagnetic showers. The central 
electromagnetic calorimeter is called CEM |152j , the forward (plug) part PEM |153| . To improve the 
spacial resolution a layer of proportional wire chambers is located at the shower maximum in the 
electromagnetic calorimeter. Additional proportional chambers located between the solenoid and the 
CEM sample the early development of electromagnetic showers in the magnet. The outer parts of the 
calorimeters measure hadronic showers. There are three hadron calorimeter subsystems: the central 
(CHA) jl54j . the wall (WHA) and the plug (PHA) calorimeter. Table [6] summarises information on 
CDF calorimetry including the coverage in pseudorapidity, the radiation thickness and the energy 
resolution. 

Muon System In general, high-p^ muons can traverse the calorimeter loosing only a small fraction 
of their energy due to ionization. They act as minimum ionizing particles. Most hadrons, on the 
other hand, interact strongly within the calorimeter volume and produce hadronic showers. This effect 
allows to identify muons by placing additional scintillation counters and wire chambers surrounding the 
calorimeter. The CDF muon identification system is composed out of four subdetectors: the Central 
Muon Detection System (CMU) |,155j . the Central Muon Upgrade (CMP), the Central Muon Extension 
(CMX) and the Intermediate Muon System (IMU). The muon system covers the pseudorapidity region 
of 1 77 1 < 1.5. To reach the muon system the muons must have a minimum transverse momentum of 
1.4 GeV/c. 



Top quark physics in hadron collisions 



32 



Table 6. Summary of CDF calorimeter properties in Run II. The energy resolutions for the 
electromagnetic calorimeters are for incident electrons and photons; in case of the hadron calorimeter 
for incident isolated pions. The © signifies that the constant term is added in quadrature. The 
transverse energy Et and the energy E are measured in units of GeV. 
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Trigger and Data Acquisition The CDF trigger and data acquisition system features three distinct 
levels. Level I is implemented in customized hardware, it uses input from the muon system, the 
calorimeters and COT tracks reconstructed by the extra- fast tracker (XFT). The maximum accept 
rate (i.e. output rate) of Level I is 50 kHz. Level 2 is a software trigger implemented on an Alpha 
processor. It allows to identify physics objects like electrons, photons, muons and jets. In addition, it 
is possible to select events with tracks which have large impact parameter; a feature which resembles 
a revolution in enriching samples for bottom- and charm-physics at hadron colliders. The maximum 
output rate of Level 2 is approximately 300 Hz. The highest trigger level. Level 3, runs part of the 
offline reconstruction code on a PC farm and classifies events in different data streams. The output 
rate to permanent storage is about 75 Hz. In most analyses top quark events are selected by requiring 
a high-pT muon or electron. The respective triggers are quite simple and have high trigger efflciencies 
between 95 and 100%. 

Luminosity Measurement An important set of measurements at a hadron collider is to determine the 
cross sections for the production of heavy particles, such as b quarks, top quarks or W and Z bosons. 
To accomplish these measurements it is crucial to know the integrated luminosity of a given data set 
with good precision. At CDF the luminosity measurement is performed by Cherenkov Luminosity 
Counters (CLC) |156j . The detector consists of long conical, gaseous Cherenkov counters that are 
installed at small polar angles in the proton and antiproton directions and point to the collision 
region. The Cherenkov counters measure the number of particles in the CLC acceptance as well as 
their arrival time for each bunch crossing. 

The number of pp interactions per bunch crossing follows a Poisson distribution with mean /x. 
The probability of empty crossings is given by 7^(0) = e^^. By measuring the fraction of empty 
crossings with the CLC the average /x is determined [145j . The instantaneous luminosity is then given 
by £ = / •/J./cTinei, where / is the average bunch-crossing frequency and Uinei the total inelastic pp cross 
section. At the Tevatron in Run II we have / — 1.7 MHz. To avoid a trigger bias the events for the 
luminosity measurement are taken with a random beam-crossing trigger (zero bias trigger) running 
at approximately 1 Hz. The luminosity measurement has a systematic uncertainty of 6% mainly due 
to the uncertainty on the total inelastic cross section and the acceptance of the CLC. 
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Tracking Detectors 

Figure 15. One quadrant of the D0 detector. 
4-3. The D0 Experiment 

The D0 detector is a large general purpose detector primarily designed to study high mass states 
and large transverse momentum phenomena. The detector design was optimised with respect to 
three general goals: (1) excellent identification and measurement of electrons and muons, (2) good 
measurement of parton jets at large pt and (3) a well-controlled measure of missing transverse energy 
(^t) as a sign of the presence of neutrinos. Figure [T5l shows one quadrant of the D0 detector. In the 
following paragraphs we discuss the subsystems of D0 in more details. 

Tracking System The inner tracking system comprises two major subsystems: a high precision silicon 
microstrip tracker (SMT) [157] and a scintillating fibre tracker |158j . The inner tracker is encased in a 
superconducting solenoid, which provides a magnetic field of 2 Tesla. The SMT consists of six barrel 
modules and 16 disk modules. Each barrel module has four radial layers of strip detector ladders 
which are arranged parallel to the beam axis. The layers are evenly spaced between radii of 2.5 and 
10 cm. Single- and double-sided detectors provide measurements in the r-cj) and the r-z views. Layers 
2 and 4 feature double-sided detectors which have axial strips on each side with a small stereo angle of 
2° between the views. There is one disk module with 12 wedge shaped double-sided detectors at one 
end of each barrel. These disks are called F-disks. The silicon wafers on a disk module are oriented 
perpendicular to the beam axis and allow tracking at high pseudo-rapidities. At each end of the 
barrel assembly there are additional three F-disk modules. Thus, there is 12 F-disk modules in total. 
Further along the beam pipe, on each side of the detector, there are two larger disks of single-sided 
strip detectors (H-disks). 

The Central Fibre Tracker (CFT) consists of scintillating fibres mounted on eight concentric 
cylinders which extend radially from 19.5 cm to 51.5 cm. The detector is divided into 80 sectors in 
azimuth. Each sector contains 896 fibres, resulting in a total of 71680 channels. The CFT provides 
full angular coverage of the central region up to < 1.7. Charged particles passing through a fibre 
produce scintillation light that will travel along the fibre in both directions. At one end, an aluminium 
mirror reflects the light back into the fibre. At the other end, the fibre is joined to a wavelength 
shifting waveguide that transmits the light to silicon avalanche photodiodes, which convert the light 



Top quark physics in hadron collisions 



34 



into an electrical signal. The momentum resolution of the D0 tracking system can be parametrized 
by Apx Ipt — ■\/0.015^ + (0.0014^^)^, where pT is measured in GeV. 

Calorimetry D0 has a sampling calorimeter based on depleted uranium, lead and copper as absorber 
materials and liquid argon as sampling medium. The calorimeter consists of three cryostats: the Centre 
Calorimeter (CC), which covers the pseudorapidity region of \r]\ < 1.2 and two endcap calorimeters 
(EC) which extend up to |?7| w 4. The electromagnetic section of the calorimeter has a thickness 
of approximately 20 radiation lengths Xq. Calorimeter cells in the central region have a size of 
Aij X Acj) = 0.1 X 0.1, except for the third layer of the electromagnetic section, where the grid size 
is A?7 X A(j) = 0.05 x 0.05 to provide better spacial resolution in the maximum of electromagnetic 
showers. At 77 = 0, the CC has a total of 7.2 nuclear interaction lengths (A), at the smallest angle of 
the EC the total is 10.3 A. The energy resolution of the calorimeter is AE/E = 15%/ VE 0.4% for 
electrons and photons. For charged pions and jets the energy resolution is 50%/ VE and 80%/-\/E, 
respectively {E is measured in GeV). 

To improve electron identification D0 features a preshower detector which consists of two 
subsystems: the central system covering the region with |?7| < 1.3 (CPS) and the forward systems 
covering 1.5 < \ri\ < 2.5 (EPS) [159] . Both subdetectors are based on a combination of lead radiator 
and scintillator layers. Electrons and photons shower in the radiator, while muons or charged pions 
only deposit energy due to ionization. In the central region the lead is mounted on the outer surface 
of the solenoid amounting to a total of 2Xq. The scintillator is placed between the lead layer and the 
calorimeter. In the forward region, which is not fully covered by the tracking system, the radiator 
is sandwiched between two scintillator layers. Photons do not deposit energy in the first scintillator 
layer and can thereby be distinguished from electrons. 

Muon System The muon detector of D0 consists of three major components: (1) a toroid 
magnet which provides a magnetic field of \B\ ~ 1.8 T, (2) the Wide Angle Muon Spectrometer 
(WAMUS) [160] covering the central region of \r]\ < 1.0, and (3) the Forward Angle Muon Spectrometer 
(EAMUS) covering the area of 1.0 < |?7| < 2.0 |161j . Both spectrometers are based on layers of 
proportional drift tubes and scintillation counters. Being immersed in a magnetic field allows the 
muon system to provide a momentum measurement for muons independent of the inner tracker. This 
feature of D0 is a major difference to CDF where the muon momentum is determined by matching 
a muon stub to a track in the central tracker. The total amount of material in the calorimeter and 
the iron toroids varies between 13 and 19 nuclear interaction lengths, making the background from 
hadronic punch-through negligible. 

Trigger and Data Acquisition D0 has a three-tier system of triggers to select events for recording. 
The Level 1 trigger is a hardware based system with an output rate of 10 kHz |162j . Data coming 
from the preshower detector and the fibre tracker (CFT) are combined to form the Central Track 
Trigger (CTT) |163j . In addition, there are conventional calorimeter and muon triggers. Level 2 of 
the trigger system delivers an output rate to 1 kHz [164j . At the second trigger level information of all 
detector components can be combined at improved resolutions. Level 3, running high level software 
algorithms on a PC farm, further reduces the rate to 50 Hz. 

Luminosity Measurement At D0 the luminosity is measured by two hodoscopes located on the inside 
face of the calorimeters |165j . Each of these hodoscopes is made of 24 wedge shaped scintillators with 
fine mesh photo-multiplier tubes mounted on the face of each wedge. The hodoscopes cover the 
pseudorapidity region of 2.7 < \ri\ < 4.4. The luminosity counters are designed to distinguish between 
single and multiple interaction events. The measurement of the rate of single inelastic interaction 
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events is used to deduce the Poisson mean /i of the number of interactions per bunch-crossing. The 
instantaneous luminosity is then calculated by using C — f ■ [ijainei as discussed in Sec. 14.21 The 
luminosity measurement has an uncertainty of 7%. 

J^.^. Top quark signatures in tt events 

In this section we discuss the experimental signature of top quark events. We constrain the discussion 
to the decay mode which dominates in the SM: t ^ W + b with a branching ratio close to 100%. 
We further concentrate on the signatures of tt events, since pair production is the main source of 
top quarks at the Tevatron and at the LHC. The signature of single top quark events and top quark 
decays via flavour changing neutral currents are presented in chapters 18.61 and [HI 

Once both top quarks have decayed, a ti event contains two W bosons and two b quarks: 
W^W~bb. Experimentally, tt events are classified according to the decay modes of the W bosons. 
There are three leptonic modes (ez^e, fiVp,, tVt) and six decay modes into quarks of different flavour 
(ud, us, ub, cd, cs, cb). We distinguish four event categories: 

(i) Both W bosons decay into light leptons (either ev,, or ^Vfj) which can be directly seen in the 
detector. This category is called dilepton channel. 

(ii) One W boson decays into ei^e or fii/^. The second W decays into quarks. This channel is called 
lepton-plus-jets. 

(iii) Both W bosons decay into quarks. We refer to this mode as the all hadronic channel. 

(iv) At least one W boson decays into a tau lepton (ti't) which itself can decay either leptonically 
(into eVe or /xf^) or hadronically into quarks. 

In good approximation we can neglect all lepton masses with respect to the W mass and write: 

At lowest order in perturbation theory the decay rate into a quark-antiquark pair, qiq2, is given by 
the rate into leptons T^y multiplied by the square of the CKM matrix element iVg^g^P and enhanced 
by a colour factor of 3, which takes into account that quarks come in three different colours: 

T{W^qiq2) = 3\Vg,g,\^Tl, . 

The hadronic decay width Fhad of the W is summed over all six quark-antiquark modes 

Thad = 3r^ • ^ iKjigJ^ with qiq2 e {ud,us,ub,cd,cs,cb}. 

9192 

The unitarity of the CKM matrix demands that the sum over the CKM elements squared is equal to 
2 and one obtains Fhad — 6 T^,. As a result, each leptonic channel has a branching ratio of 1/9, while 
the hadronic decay channel into two quarks has a branching ratio of 6/9. For the tt decay categories 
we get thus the probabilities as listed in Tab.[71 (i) Dilepton mode: 4/81, (ii) lepton-plus-jets: 24/81, 
(iii) all-hadronic channel: 36/81, (iv) tau modes: 17/81. These four different types of events can 
be isolated by their distinct event topologies. 

Dilepton Channel This final state includes (1) two high px leptons, electron or muon, (2) a large 
imbalance in the total transverse momentum (missing transverse energy, ^t) due to two neutrinos, 
and (3) two b quark jets. The dilepton event category has low backgrounds, especially in the e/x 
channel, since Z° mediated events do not contribute. However, the drawback of the dilepton channel 
is its low branching ratio of about 5%. There is a small contribution from tau events to the dilepton 
channel, if the tau decays into e or fi. This cross feed has to be taken into account when calculating 
acceptances for dilepton analyses. Since dilepton events contain two neutrinos which contribute to the 
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Table 7. Categories of tt events and their branching fractions. 
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Figure 16. Event display of a CDF II muon-plus-jets event. The isolated line on the left hand side 
shows the muon trajectory. The arrow on the lower right indicated the direction of the .^-p. The 
event features four hadronic jets. 

l^T, the top quarks cannot be fully reconstructed. This is a drawback if one wants to measure the top 
quark mass. However, this disadvantage is partially compensated by the precisely measured lepton 
momenta, in contrast to the only fair measurement of jet energies in the lepton-plus-jets channel. 

Lepton-Plus-Jets Channel The lepton-plus-jets channel is characterised by (1) exactly one high-pr 
electron or muon, which is also called the primary lepton, (2) missing transverse energy, (3) two b 
quark jets from the top decays, and (4) two additional jets from one W decay. An event display of a 
lepton-plus-jets candidate event measured at CDF II is shown in figure [161 The big advantage of the 
lepton-plus-jets channel is its high branching fraction of about 30%. The backgrounds are considerably 
higher than in the dilepton channel, but still manageable. Several strategies of background suppression 
have been developed and are discussed in section [^?^ Again, as in the dilepton channel, there is some 
cross feed from tau modes which has to be taken into account for acceptances. In lepton-plus-jets events 
the momentum of the leptonically decaying W boson can be reconstructed up to a twofold ambiguity. 
The transverse momentum of the neutrino is assumed to be given by • Two solutions for the z 
component of the neutrino are obtained from the requirement that the reconstructed invariant mass 
of the lepton and the neutrino be equal to the well known W mass: Mti, — Mw- To fully reconstruct 
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the momenta of the top and antitop quark in the event, one has to assign the measured hadronic jets 
to the quarks. Without identification of b quark jets there are 24 possible combinations, including the 
ambiguity of the neutrino reconstruction. If one jet is identified as likely to originate from a b quark 
12 combinations remain. If there are two identified b jets the ambiguity is down to 4 options. This 
illustrates the importance of b quark jet identification for top quark physics. 

All Hadronic Channel The all hadronic channel has the largest branching ratio out of all tt 
event categories, about 44%. However, backgrounds are quite considerable. While there is the 
advantage that all final state partons are measured, one has to deal with numerous combinations 
when reconstructing the top quark momenta. Even if two jets are identified as b quark jets, 12 
possible combinations remain. Another drawback is that jet energies as measured in the calorimeter 
have large uncertainties. The combination of these disadvantages leads to the conclusion that the all 
hadronic channel, even though a clear tt signal has been established here, proves not very useful to 
further investigate top quark properties. 

Tau Modes Top quark events containing the decay W — s- rv-r are difhcult to identify and have not 
been seen to date. The search for this decay mode is briefly described in section [8731 

4-5. Particle detection and identification 

In the previous section we have explained tt signatures in terms of primary partons. In this section 
we will discuss how these partons are detected and how their properties are measured. 

4-5.1. Electrons and muons Charged leptons emerging from a decay oi a W boson have high 
transverse momenta that can be measured by the tracking system with good resolution. Electrons 
interact within the first few segments of the calorimeter and form electromagnetic showers of photons 
and electron-positron pairs. To select electrons from W decays one typically requires a transverse 
energy, Et = E ■ sin0, of at least 20 GeV. The energy resolution (at CDF for example) is 
a{ET)/ET = 13.5%/ y/Er © 2%, where Et is measured in GeV. The shower should be mainly 
contained in the electromagnetic part of the calorimeter. A typical requirement is that 90% of the 
total measured energy is observed in the electromagnetic calorimeter. Moreover, there must be a 
track that, if extrapolated to the calorimeter, matches the location of the electromagnetic shower and 
has a momentum consistent with the shower energy (E/p ^ 1). Additional requirements include the 
transverse and longitudinal shapes of the shower, which are known from test beam data. Energetic 
photons that interact with detector material prior to their entry into the tracking system and produce 
an electron-positron pair can mimic high-p^ primary electrons and pose a serious background. One can 
identify and veto these conversion pairs if both, the electron and the positron track, are reconstructed 
and fulfil a matching criterion with respect to a common vertex. The major background to high-pr 
electrons arises from hadronic showers developing early in the calorimeter and depositing most of their 
energy in the electromagnetic section. Such a case can, e.g., be caused by photons from tt^ decays. The 
hadronic jet background can be considerably suppressed by exploiting the fact that leptons coming 
from a heavy boson decay are isolated from other jet activity in the event. Thus, it is asked that the 
electron shower is isolated from other energy deposits in the calorimeter. A typical requirement is 
that the energy measured in an annulus of radius AR — 0.4 around the electron shower is less than 
10% of the shower energy. 

Muons can be reliably identified as tracks that penetrate the calorimeter as well as additional 
shielding material and reach the outmost layers of the collider experiment, the muon system. 
The track is typically required to have a transverse momentum of at least 20 GeV/c and match, 
when extrapolated through the calorimeter, a short track segment measured in the muon system. 
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Uncertainties due to multiple scattering in the calorimeter and the shielding material have to be 
considered accordingly in the matching criteria. The transverse momentum of high-pj^ muons is 
measured with a typical resolution of <7{pt)/pt = 0.1%^^, where pt is given in GeV/c. The 
energy measured in the calorimeter segment containing the muon is required to be consistent with 
the deposition expected from a minimum ionizing particle. Muons are required to be isolated from 
additional calorimeter activity as described above for electrons. Real muons from cosmic rays have to 
be distinguished from muons originating from hadronic collisions and removed from the data sample. 
One requires that the muon track extrapolates to the primary collision point of the event. Timing 
information of the tracker and the calorimeter are used to ensure that the passage of the muon 
through the primary interaction region falls into a narrow time window around the beam crossing, 
that is known from the accelerator clock signal. Other backgrounds to the muon signal arise from 
pion or kaon decays in flight, or from hadrons that traverse the calorimeter and the absorber material 
without producing a hadronic shower (punch-through). 

High-pT electrons and muons can be speedily reconstructed with good precision and are used in 
the trigger systems of collider experiments to select events containing heavy gauge bosons and other 
high-pT events in real time. The primary data samples for tt dilepton and lepton-plus-jets events are 
defined in this way. 

4-5.2. Neutrinos Neutrinos interact only weakly with matter and therefore cannot be directly 
observed in a collider detector. Instead their presence is inferred from an imbalance in the total 
transverse energy of the event. As described in section r3.1. li the hard scattering in hadronic collisions 
happens between two partons whose momentum fraction is a-priori unknown. The remnants of the 
colliding hadrons, the spectator partons that do not participate in the hard interaction, have little 
transverse momentum and escape undetected down the beam pipe. In contrast to e+e^ collisions one 
can therefore not simply invoke energy and momentum conservation before and after the collisions. 
However, the total transverse momentum is conserved and is known to be zero before the collision. 
Any imbalance in the vector sum of transverse momenta can therefore be attributed to the presence of 
neutrinos that carry away momentum undetected. In practice, one cannot determine the momenta of 
all particles produced in the collision, but rather measures the imbalance of energy in the calorimeter. 
To each calorimeter cell i one assigns a transverse energy vector = {Ei sinOi cos^i, Ei sinfl^ sin^i) 
and calculates the sum over all cells: = ~X]i-^T- Here 6i and (fii are the angular coordinates 
of calorimeter cell i. Since muons loose only a minimal amount of their energy in the calorimeter, 
has to be corrected for identified muons. The limiting factor on the resolution of — I-^t | 
is the uncertainty in the measurement of hadronic jet energies in the calorimeter. Therefore, one 
usually parametrizes the resolution in terms of the scalar sum of all transverse energies in the 
event, '^Et, measured in GeV. For tt dilepton candidate events CDF has measured a resolution of 
(t(-^t) = 0.7 • Et |166j . D0 determined its resolution for minimum bias data and quotes 
ct(^t) 1-08 GeV + 0.019 J^^t [MT]. 

4.5.3. Jets of quarks and gluons The strength of the strong force increases with distance and prevents 
elementary particles carrying colour charge, i.e. quarks and gluons, from existing as free objects. 
Quarks and gluons are confined to exist in hadrons, which are colour singlets. In collider experiments 
quarks and gluons participate in the hard interaction and subsequently form coUimatcd jets of hadrons 
which can be observed in a detector. Jets tend to preserve the direction of motion of the original parton. 
In the calorimeter they are detected as an extended cluster of energy. To compare measurements with 
theoretical predictions a precise mathematical prescription is necessary how to form jets out of energy 
clusters measured in the calorimeter. Different jet algorithms are available, but at hadron colliders 
calorimeter cells are conventionally combined within a cone of fixed radius Ai? in r/-0 space, because 
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jets are approximately circular in the rj-cj) plane. Moreover, in this paranietrization the size of a jet of 
a particular E-^ is independent of the jet rapidity. The cone size Ai? is chosen differently depending 
on the physics analysis. On one hand the cone must be big enough, such that most of the energy 
associated to the original quark is contained in the jet cone. On the other hand the cone size should 
be small enough, such that energy depositions corresponding to different primary partons are resolved 
individually rather than merged to one jet. This is especially a concern for tt events which have many 
jets in their final state. At the Tevatron the optimal choice was found to be Ai? = 0.4 for CDF and 
Ai? = 0.5 for D0. The difference between the two experiments reflects a small dependence on the 
calorimeter design and geometry. 

Both Tevatron experiments feature an energy resolution for jets of a{ET)/ET ~ 100%/\/Et, 
where Et is measured in GeV. Several systematic effects compromise the jet energy measurement: 
(1) intrinsic large fluctuations in the response of calorimeters to hadronic showers, (2) nonlinear effects 
in calorimeter response, (3) calorimeter non-uniformities and energy loss in uninstrumcnted regions, 
such as cracks between modules, (4) increase in energy due to the overlap of multiple hard interactions 
in one bunch crossing, (5) energy of the underlying event feeding into the jet cone (The underlying 
event consists of particles coming from the fragmentation of partons that do not participate in the 
hard scattering of the primary pp or pp interaction.), and (6) energy loss due to the use of a finite cone 
size in jet reconstruction. To understand these effects experimentally turns out to be one of the major 
challenges in collider physics. Appropriate correction methods have to be derived and the detector 
simulation has to be tuned to describe the measurements. 

Calorimeter calibration commonly starts with the electromagnetic section for which an absolute 
energy scale can be derived by comparing electron energies measured in the calorimeter to their 
momenta measured in the tracking system or by reconstructing resonances with well known masses 
such as the Z boson, the 7r° or the J/tp- In a second step the hadronic part of the calorimeter 
can be calibrated against the electromagnetic part by studying photon-plus-one-jet events, where the 
transverse energies should balance. Dijet events are further used to calibrate one hadronic region of 
the detector relative to another better understood region. Contributions from the underlying event are 
investigated by Monte Carlo simulations and comparison to data taken under special trigger conditions 
and luminosities. 

Finally, it has to be noted that one cannot achieve a one-to-one correspondence of primary quarks 
and gluons and the observed jets in all kinematic regions and for all event topologies of a certain 
physics process, and one cannot expect to do so. To define jets a minimum of transverse energy is 
required. One may start counting soft jets with 8 GeV. However, this can complicate discrimination 
between signal and background. In most ti analyses jets are typically defined with Et > 15 GeV. 
This threshold causes, e.g., some ti lepton-plus-jets events to have only three instead of four jets. 
Another source of inefficiency in jet reconstruction is the merging of jets which are close together 
in r]-(j) space and cannot be resolved separately. Conversely, if a parton radiates a gluon with large 
relative transverse momentum, that gluon can be reconstructed as an additional jet. 

4- 5.4- Tagging of b quark jets As pointed out in section 14.41 tt events feature two energetic b 
quark jets, while the heavy flavour content in background events is relatively low. Therefore, the 
identification, also called the tagging, of 6 quark jets is an important tool to suppress background, 
enrich top quark data samples, and to facilitate the full reconstruction of top quark momenta by 
reducing the number of possible combinations of final state objects. In this section we briefly present 
three b tagging methods employed at hadron colliders. The flrst two methods are based on the 
relatively long lifetime of b hadrons, which is about 1.5 ps. Since b hadrons emerging from top quark 
decays have relatively high momenta, their long lifetime allows them to travel several mm before 
decaying. Tracks from a b hadron decay therefore typically originate from a secondary vertex that is 
displaced from the primary interaction point. The first b tagging method reconstructs the secondary 
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vertex of b hadrons within a jet and bases the h tag on the significance of the displacement. The 
second method searches for tracks with large impact parameter with respect to the primary vertex, 
but does not require the reconstruction of the secondary vertex itself. The third method identifies 
leptons at intermediate momenta (typically: 2GeV/c < px < 20GeV/c) from semileptonic h decays 
within jets. Since momenta of these leptons are much smaller than those of leptons coming from W 
or Z boson decays, this method is also called soft lepton tag. 

Secondary vertex tag The secondary vertex method identifies b jets by establishing a displaced 
secondary vertex in the jet from the decay of a long lived b hadron. The following explanation 
goes along the lines of the algorithm employed by CDF, but the D0 version is quite similar. In a 
first step the space point of the hard primary interaction is precisely reconstructed. One starts by 
clustering the z coordinates of all tracks at the point of their closest approach to the origin (perigee) . 
After applying several quality requirements this procedure yields an estimate of the z position of the 
primary vertex, Zpv. The trajectory of the proton and antiproton beams through the collision region, 
also called the beamline, is known with high precision, 0(1 /im), from subsidiary measurements made 
on a run-by-run basis. Combining Zpv with the beamline yields a first estimate of the primary vertex 
position The precision of ajpv is subsequently improved by including tracks with low impact 

parameter significance into a full vertex fit. As a result, the uncertainty on the transverse position of 
the primary vertex ranges from about 10 - 30 /xm, depending on the number of tracks and the event 
topology. 

The subsequent steps of secondary vertex tagging are performed on a per-jet basis. Only tracks 
within the jet cone are considered. To ensure good track quality cuts involving the transverse 
momentum, the significance of their impact parameter d relative to Xp^, d/ad, the number of silicon 
hits attached to a track, the quality of those hits, and the quality of the final track fit are applied. 
Tracks consistent with coming from the decays — > 7r"'"7r~ or A — + tt^p, or photon conversions are 
not accepted as good candidate tracks. For a jet to be taggable at least two good tracks have to fall 
within the jet cone. One iterates over the set of good tracks and tries to fit them to a common vertex 
with a certain minimum fit quality (x^ per degree of freedom). 

Once a secondary vertex is found in a jet, the two-dimensional decay length L2D of the secondary 
vertex is calculated. L2D is the projection of the vector pointing from the primary to the secondary 
vertex onto the jet axis in the xy plane. The sign of L2D is determined by the angle between the jet axis 
and the secondary vertex vector in the xy plane. If the angle is < 90° the sign is positive, if the angle is 
> 90°, the sign is negative. Large positive values of L2D are predominantly attached to vertices of real 
b hadron decays, while displaced vertices with negative L2D are mainly due to track mismeasurements 
or random combinations of tracks. To reduce the background from false secondary vertices, called 
mistags, a good secondary vertex is required to have a minimum decay length significance of typically 
i2D/o'2D > 3, where cr2D is the estimated uncertainty on L2D- The uncertainty cr2D is estimated for 
each vertex individually, a typical value is about 200 /zm. Figure flTb shows the distribution of £20 as 
obtained in a CDF measurement of the tt cross section [168| . The distribution has the form of a falling 
exponential, as expected from a lifetime distribution. The content of the first bin is exceptionally lower 
due to the cut on the decay length significance, i2D/o'2D > 3. More detailed descriptions of secondary 
vertex b tagging can be found in references |166|,I103] (CDF Run I), reference |168j (CDF Run II) and 
references |1691 1170) (D0 Run II). In the latest CDF analysis (Run II) the efficiency of tagging at 
least one b quark jet in a tt event is (53±3)% |168j . This number includes the geometric acceptance 
of tracks from b decays (taggability). The estimated uncertainty is purely systematic. 

Impact parameter tag The impact parameter b tag algorithm compares track impact parameters to 
measured resolution functions in order to calculate for each jet the probability that there are no long 




Figure 17. (a) Two-dimensional decay length L2D of secondary vertices as obtained in a CDF 
measurement of the tt cross section |168l . (b) Distribution of the signed impact parameter significance 
of tracks in a jet data sample, where there is at least one jet with Et > 50GeV per event. The 
data are fit with a resolution function consisting of two Gaussians plus an exponential function, 
separately for the positive and negative sides |103j . 



lived particles in the jet cone. For light quark or gluon jets this probability is uniformly distributed 
between and 1, while it is very small for jets containing tracks from displaced vertices of heavy flavour 
decays. The impact parameter significance 6*0 is defined as the ratio of the impact parameter d and its 
uncertainty ad- The sign of d is defined to be positive if the point of closest approach to the primary 
vertex lies in the same hemisphere as the jet direction, and negative otherwise. Figure fTTb shows 
the 5*0 distribution of tracks in a jet data sample. The distribution is fit with a resolution function. 
The negative side of the resolution function is used to compute the probability P{So) that the 5*0 of 
a particular track is due to detector resolution. The probability that a jet is consistent with a zero 
lifetime hypothesis is essentially given by the product of the individual probabilities P{So) of tracks 
contained in the jet cone including an appropriate normalization. Tracks used in the calculation 
of the jet probability are required to satisfy certain quality criteria involving a cut on the impact 
parameter, the transverse momentum and the number of hits in the silicon detector. The impact 
parameter tagging method has the advantage of providing a continuous variable, the jet probability, 
to distinguish between light and heavy flavour jets, which easily allows it to choose a certain value 
of b tagging efficiency by adjusting the cut value on the jet probability. In addition, the probability 
value can be used in subsequent steps of an analysis, e.g. when applying multivariate techniques. 

Soft Lepton Tag The Soft Lepton Tag (SLT) is based on the identification of leptons originating from 
semileptonic b decays b —^ £ + X, where the lepton is either an electron or a muon. Contributions 
come from the decay mode b — > ciiy or the sequential decay 6 — > c — > siv. The semileptonic branching 
ratio of b and c hadrons is measured to be about 10% per lepton species (e or fj,) |10| . This means 
that about 80% of tt events have at least one lepton from a semileptonic decay. In contrast to leptons 
from W decays the leptons from b decays are not isolated, but are rather contained in the cone of 
the b quark jet. Their px spectrum is much softer, typically well below 20GeV/c. Therefore, the 
detection and identification of these leptons is much more difficult. In Run II, D0 and CDF have 
used only muons for the soft lepton tag so far. D0 selects muons in the transverse momentum range 
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of 4GeV/c < pt < 15 GeV/c and achieves a b jet tagging efficiency of 11% |170j . 
4-6. Generation and simulation of Monte Carlo events 

Collider experiments are intricate devices and predicting the experimental response to a certain physics 
process is non-trivial. To extract sensible information from data one needs first an accurate modelling 
of the event kinematics and topology on parton and hadron level, and second a detailed simulation of 
the detector response to particles interacting with the detector material. 

There are two commonly used general purpose Monte Carlo event generators, Pythia [171j 
and Herwig [172j . which provide an exclusive description of individual events at hadron level. 
They are based on leading order matrix elements for the hard parton scattering convoluted with 
parametrizations of parton distribution functions and include approximate treatments of higher order 
perturbative effects, initial and final state gluon emission, parton shower, hadronization, secondary 
decays and the underlying event. The strength of Pythia and Herwig is the modelling of the parton 
shower, where outgoing partons are converted into a cascade of gluons and qq pairs with energy and 
angular distributions determined by the DGLAP equations, which describe the evolution of quark 
and gluon densities [173[ I174i I175j . The parton shower is terminated when the virtual invariant mass 
of the parton (Q^) falls below a threshold value chosen such that a perturbative treatment is valid 
above the threshold. The evolution of the shower beyond this stage is determined by non-perturbative 
physics. The partons are turned into colourless hadrons according to phenomenological hadronization 
(or fragmentation) models. Particles originating from interactions of the beam remnants are also 
included in the model and are called the underlying event. A set of parameters is tuned to reproduce 
hadron multiplicities and transverse momentum spectra of the underlying event as measured in soft 
pp collisions, so called minimum bias events. To obtain a proper modelling of b and c hadron decays, 
heavy flavor jets are interfaced to the decay algorithm QQ which comprises tables of the most up to 
date branching fractions. 

While Pythia and Herwig provide a good description of tt events |176| . they are insufficient to 
model multijet and W + multijet background events. Particularly critical is the generation of + 
heavy flavour jets. In Run 1 of the Tevatron W + jets events were generated with the Vecbos Monte 
Carlo program |177j , in Run II the Alpgen program is used [178; , which generates high multiplicity 
partonic final states based on exact leading order matrix elements. The parton level events are passed 
to Herwig for showering and hadronization and to QQ for the decay of heavy flavour hadrons. To 
reproduce the entire W + multijets spectrum several Alpgen samples have to be accurately merged 
together. 

Monte Carlo event generators output a list of stable particles which can be fed to a full detector 
simulation that reproduces the interaction of those particles with the detector material. The simulation 
code is based on the Geant package [179] which allows to implement a detailed geometry description 
of the detector, track particles through the given detector volumes and compute their energy loss as 
well as multiple scattering. In the so called digitization step the energy depositions obtained from 
Geant are converted into raw detector data which are in the same format as physics data recorded 
from real collisions. Most detectors are described by parametric models that convert energy deposits 
into detector signals. The parameters of the models are tuned to correctly reproduce the detector 
response. The silicon detector response, for example, can be modelled by a simple geometrical model 
based on the path length of the ionizing particle and a Landau distribution measured from physics 
data. The raw data obtained from Monte Carlo events are then subjected to the same reconstruction 
code as data from collisions. These reconstructed Monte Carlo data can then be analysed in the same 
way as if they would contain physics events. Monte Carlo samples generated in this way are used to 
compute acceptances for certain physics processes and to compare kinematic distributions between 
data and predictions. 
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5. The quest for the top quark 

Immediately after the discovery of the b quark in 1977 the existence of a weak isospin doublet partner, 
the top quark, was hypothesised. The mass of the sixth quark was unknown and a wealth of predictions 
appeared based on many different speculative ideas, see for example references |1801ll811ll82| . Typical 
expectations were in the mass range of about 20 GeV/c^, which became accessible two years later with 
measurements at the PETRA e''"e~ collider. 

5.1. Early searches 

In e'^e~ annihilations tt pairs can be produced by the exchange of a virtual photon or Z boson. The 
mass hierarchy of Standard Model quarks is beautifully reflected in a measurement of the ratio of the 
production rate of hadrons to that for muon pairs, R — (T(e+e^ hadrons)/cr(e+e^ jj.'^jj,^), as a 
function of the centre-of-mass energy y/s. The ratio R is given by i?(s) = 3 J^Qh where the sum 
is over all quark flavours with production thresholds greater than y/s. This corresponds to a step 
function, which is only altered near production thresholds where strong resonance effects occur. If 
one compares R above the tt production threshold to its value below (but above the threshold for bb 
production) one expects to see an increase by AR = 4/3. 

Between 1979 and 1984 the five experiments at the PETRA colhder, CELLO, JADE, MARK J, 
PLUTO, and TASSO, measured R at several centre-of-mass energies ranging from 12 to 46.8 GeV 
in steps of about 20 to 30 MeV [HSl HHH [1851 UlSl [M]. No indication of top quark production 
was found. The final lower limit of the top quark mass derived from PETRA measurements was 
23GeV/c2. 

At the KEK laboratory in Japan a dedicated accelerator, TRISTAN, was built to search for the 
top quark in e+e^ annihilations, with slightly increased centre-of-mass energy than PETRA. The 
maximum energy reached was 61.4 GeV. Between 1987 and 1990 several consecutive searches for top 
quark production were presented by the two experiments operating at TRISTAN (AMY and VENUS). 
Again, no evidence for the top quark was found and the final lower limit on the top quark mass was 
30.2 GeV/c2 p8l[T89] . 

In 1989 the Large Electron Position collider (LEP) at CERN and the Stanford Linear Collider 
(SLC) started e^e" collisions at the Z'^ pole, ^/s « 90 GeV. The experiments at these accelerators - 
ALEPH, DELPHI, OPAL and L3 at LEP, and SLD at SLC - searched for Z° tt events of various 
topologies [1901 11911 11921 1193j . The best lower limit on the top quark mass obtained from these direct 
searches was 45.8GeV/c^. 

In the mid 1980s experiments at hadron colliders joined experiments at e'^e~ machines in the quest 
for the top quark. In contrast to e^e^ annihilations, searches in hadronic collisions have to deal with 
high backgrounds and thus model independent analyses are not feasible. Instead the experiments 
had to concentrate on SM signatures. Initially, one assumed the top quark mass to be below the 
mass of the W boson. Top quark production via the electroweak interaction was expected to be the 
dominating process at the CERN SppS collider operating at ^/s = 546 GeV, later at ^/s = 630 GeV. 
The experiments UAl and UA2 searched for events were an on-shell W boson is produced which 
decays into a top and a bottom quark: pp tb/tb. The top quark was reconstructed in its 

semileptonic decay mode into a b quark, a charged lepton (electron or muon) and a neutrino. 

First results reported by UAl in 1984 seemed to be consistent with the production of top quarks of 
mass (40± 10) GeV/c^ |194j . UAl observed six events with one isolated lepton, 3 electron and 3 muon 
events, and 2 hadronic jets. The effective invariant mass distribution of the lepton, the two jets and the 
missing transverse energy shows a pronounced peak at the W mass, while the mass of the lepton, the 
missing transverse energy and the jet with the lowest transverse energy clustered around 40 GeV/c^. 
The number of background events where the lepton was faked by a hadron or background events from 
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heavy flavour production were believed to be negligible. Thus, these six events were interpreted as 
a first indication of a top quark originating from a W boson decay. However, subsequent analyses 
by UAl were based on a more complete evaluation of the backgrounds and did not support this 
result [195] . Since UAl did not apply a minimum cut on the missing transverse energy in the event, 
backgrounds due the fake leptons, Drell-Yan production of lepton pairs, £'^£~ , as well as backgrounds 
from bb and cc production were found to be more significant than originally estimated. 

The final UAl top quark analysis was published in 1990 and used data corresponding to an 
integrated luminosity of 5.4 pb~^. The observed events were found to be consistent with the non- 
top backgrounds, and a lower limit on the top quark mass of 60GeV/c^ was set |196| . A slightly 
better limit was achieved by the second experiment operating at the CERN SppS, UA2, based on an 
integrated luminosity of 7.1 pb~^ [197j . UA2 used only events containing an isolated central electron 
with p^ > 12 GeV/c and > 15 GeV. After all selection cuts 137 events were retained, in good 
agreement with the expectation of 154.4 ± 14.6 background events. To further enhance the sensitivity 
of the analysis the transverse mass of the electron and the missing transverse energy 



was formed, where A.(j)ei, is the azimuthal angle between the electron and vectors. The limit on 
the top quark mass were obtained by fitting the observed M^'^ distribution to template distributions 
of background events alone or a combination of background and a top quark signal of certain mass. 
As a result, top quark masses below 69 GeV/c^ were excluded. 

It is important to note that experiments at hadron colliders do not provide a direct limit on the 
top quark mass, but rather on the top quark production cross section. The mass limits are derived 
using the theoretically predicted production cross sections, that depend on the top quark mass, and 
the predicted branching ratios. Systematic uncertainties on these predictions have thus to be properly 
accounted for. 

5.2. Searches and discovery at the Tevatron 

In 1988 the CDF experiment at the Tevatron joined the race for discovery of the top quark. Due 
to the higher centre-of-mass energy at the Tevatron of ^/s =1.8 TeV top quarks are predominantly 
produced as ti pairs. The first CDF top quark search uses a data sample with an integrated luminosity 
of 4.4 pb"^ accumulated in Run which lasted from 1988 to 1989 |1981I199| . CDF pursued a similar 
strategy as UA2 and used the M^"^ distribution to discriminate between ly+jets background events, 
where the W boson decays on shell, and top quark decays where the electron and the neutrino originate 
from the exchange of a virtual W (this was still under the assumption that Mtop < Mw + Mb). The 
analysis is based on isolated central electron events with an electromagnetic transverse energy above 
20 GeV and I^t > 20 GcV. In addition, at least two hadronic jets with Et > 10 GeV and |?7| < 2.0 
are required. This analysis pushed the lower limit of the top quark mass to A/top > 71 GeV / c^ using 
predictions of the cross section by Altarelli et al. [89l [85] . 

In addition to the search in the electron-plus-jets channel CDF performed an analysis in the 
electron-plus-muon (efi) dilepton channel [200j . By requiring two leptons from different families, 
backgrounds from Drell-Yan and Z° production as well as M^-|-jets events are strongly suppressed. 
The remaining background is mainly due to Z° t~^t~ events. Subsequently, CDF extended the 
dilepton analysis to include the ee and /i/x channels. Backgrounds from Drell-Yan and events are 
reduced by additional requirements: (1) The dilepton azimuthal opening angle A(j)u is required to be 
below 160°. (2) ee or events with a dilepton invariant mass in the window 75 < Mu < 105 GeV/c^ 
or with < 20 GeV are rejected. The lepton-|-jets analysis was also further improved. High-px 
muons from the decay of a boson were included and a soft muon b tag for at least one of the 
jets was required to reduce the MK-|-multijet background. The dilepton and the lepton-|-jets analysis 
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were combined to derive a common limit on the ti cross section cr^j. Using theoretical expectations 
for att [SO], and assuming SM decays for the top quark, the cross section limit was translated into 
a lower limit on the top quark mass of Mtop > 91GeV/c^ at the 95% confidence level |201[ 1202] . 
It is important to realize that with the top quark mass being above the sum of the W and the h 
quark mass, the transverse mass is no longer a suitable discriminant and a strategy change is 
necessary. Therefore, later CDF analyses were based on counting events only, rather than fitting the 
M^" distribution. 

In 1992, with the start of Tevatron Run I, the D0 experiment joined the hunt for the top quark. 
In April of 1994 D0 published its first top quark analysis setting the last lower limit on the top quark 
mass before its discovery j203j . The data sample was recorded in 1992/93 and corresponds to an 
integrated luminosity of 15 pb~^. In this analysis D0 uses the and the ee mode of the dilepton 
channel and the lepton-plus-jets channels. In e/j, events one electron with i?T > 15 GeV and one muon 
with pt > 15 GeV/c, > 20 GeV and at least one jet with i?T > 15 GeV are required. For ee 
candidates tighter cuts are applied, the two electrons must have i?T > 20 GeV and > 25 GeV. 
To suppress Z° — > e+e~ background the cut on missing transverse energy is raised to 40 GeV within 
a mass window of \Mf.e — Mz\ < 12 GeV/c^, where Mee is the dielectron invariant mass. In the 
electron+jets mode one electron with Et > 20 GeV, at least four jets with Et > 15 GeV, and 
> 30 GeV are required. In muon+jets events the cut on the missing transverse energy is slightly 
relaxed: > 20 GeV. To further reduce the VF+jets background D0 performs a topological analysis 
of the events. The event shape is characterised by the aplanarity A which is a quantity proportional 
to the lowest eigenvalue of the momentum tensor of the observed objects. In e+jets events a cut of 
A > 0.08 is apphed, in ^+jets events A > 0.1 is demanded. After all analysis cuts two dilepton and 
one e+jets event are observed, consistent with background estimates. The intersection of the derived 
upper limit on the ti cross section with the theoretical prediction [92] yields a lower limit on the top 
quark mass of 131 GeV/c^ 

In 1993 and 1994 CDF saw mounting evidence for a top quark signal. The detector upgrade 
for Run I, mainly the addition of a silicon vertex detector, was the keystone for the discovery of the 
top quark at CDF. The new silicon detector allowed for the reconstruction of secondary vertices of 
b hadrons and a measurement of the transverse decay length L^y with a typical precision of 130 /zm. 
Secondary vertex b tagging proved to be a very powerful tool to discriminate the top quark signal 
against the Vl^+jets background and increase the sensitivity of the lepton-plus-jets tt analysis. In July 
1994 CDF published a paper announcing first evidence for tt production at the Tevatron based on 
events in the dilepton and the lepton-plus-jets channel [2041 1166] . The analysis uses a data sample 
with an integrated luminosity of (19.3 ± 0.7) pb~-^. Two e^i dilepton events pass the selection cuts, 
but no ee or fifj, events. The expected background in the dilepton channel is 0.56lg'^3. In the lepton- 
plus-jets channel three or more jets with Et > 15 GeV and jryj < 2.0 are required. Secondary vertex 
and soft lepton b tagging reduce the W-|-jets background. Six events with a secondary vertex tag are 
observed over a background of 2.3 ± 0.3 events. Seven events have a soft lepton tag (muon or electron) 
over a background of 3.1 ± 0.3 events. The soft lepton and the secondary vertex tagged samples have 
an overlap of three events. 

Since the cross section of M^-|-jets production cannot be reliably predicted with sufficiently small 
uncertainties, special techniques had to be developed for estimating backgrounds in the lepton-plus- 
jets search directly from data. This technique assumes that the heavy quark content (6 and c) of jets 
in the Vl^-|-jets sample is the same as in an inclusive jet sample. This assumption is a conservative 
overestimate of the backgrounds, since the inclusive jet sample contains heavy quark contributions 
from direct production (e.g. gg — > 66), gluon splitting (where a final state gluon branches into a 
heavy quark pair), and flavour excitation (where an initial state gluon excites a heavy sea quark in 
the proton or antiproton), while heavy quarks in Ty-fjets events are expected to be produced entirely 
from gluon splitting. Tag rates are measured in an inclusive jet sample and parametrized by the Et 
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Figure 18. Reconstructed top mass distributions as published (a) in the CDF evidence paper of 
1994 [2041 and (b) in the CDF discovery paper [2051 . The solid histogram shows CDF data. The 
dotted line shows the shape of the expected background, the dashed line the sum of background plus 
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and track multiplicity of each jet. The parametrization is used to derive an estimate on the number 
of expected background events based on the Vl^+jets sample before applying the h tag algorithms. 
The probability that the estimated background has fluctuated up to the total number of 12 candidate 
events, taking into account that three events have double tags, is found to be 0.26%. This corresponds 
to a 2.8(7 excess for a Gaussian probability function. 

Assuming that the excess of h tagged lepton-plus-jets events is due to it production a value 
for the top quark mass is estimated using a constrained kinematic fit. In the jets sample 

one additional soft jet with > 8 GeV and \r]\ < 2.4 is required. Seven events of the b tagged 
lepton-plus-jets sample pass this requirement and are fitted individually to the tt hypothesis. For 
each event the 12 possible combinations of the jets, the lepton and the missing transverse energy are 
considered. The combination with the best is chosen. The resulting top mass distribution, shown in 
figure [T5k. is fitted to a sum of the expected distributions from W^-Fjets and production for different 
top quark masses. The fit yields a value of Mtop = (174 ± lOtH) GeV/c^. The corresponding log 
likelihood distribution is depicted in figure [T5b. In November 1994 the D0 collaboration confirmed 
the evidence seen at CDF. An update of the previous D0 analysis, now with an integrated luminosity 
of (13.5 ± 1.6) pb~^, added soft muon b tagging [2061 1167] . In total, D0 observed nine events over a 
background of 3.8 ± 0.9. 

As Run I continued more data were accumulated and finally, in April 1995, CDF and D0 were 
able to claim discovery of the top quark |205il207j . CDF used a data sample corresponding to 67 pb~^ 
and significantly improved its secondary vertex b tagging techniques. The efficiency to identify at least 
one b quark jet in a tt event with more than three measured jets was found to be (42 ± 5)%, almost 
double the previous value of the 1994 analysis. In the new analysis the background estimate is also 
considerably improved. While the mistag rate due to track mismeasurements is again measured with 
samples of inclusive jets, the fractions of VF-|-jets events that are Wbb or Wcc are disentangled from 
Monte Carlo samples, applying measured tagging efficiencies. There are 27 jets with a secondary 
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Figure 19. Fitted mass distribu- 
tion for candidate events (histogram) 
of the D0 discovery paper |207j . Over- 
laid is the expected mass distribution 
for top quark events with Mtop = 
199 GeV/c^ (dotted curve), the back- 
ground (dashed curve), and the sum of 
top and background (solid curve) for 
(a) the standard and (b) the loose event 
selection. 
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Figure 20. History of the quest for 
the top quark. The shaded histogram 
shows experimental lower limits on the 
top quark mass. Theoretical predic- 
tions based on flavour symmetries are 
shown (dots) as well as predictions 
based on electroweak precision mea- 
surements (squares). Direct measure- 
ments of Mtop by CDF and D0 are 
represented by the triangle (triangles 
pointing up: CDF; triangles pointing 
down: D0 ). 



vertex b tag in 21 W+ > 3 jets events. The estimated background is 6.7 ± 2.1 b tags. The probabiUty 
for this observation to be a background fluctuation is 2 • 10~^. The 1995 dilepton and soft lepton b 
tag lepton-plus-jets analyses are only slightly changed compared to those of 1994. Six dilepton events 
are observed over a background of 1.3 ± 0.3. There are 23 soft lepton tags observed in 22 events, with 
15.4 ± 2.0 b tags expected from background sources. Six events contain both a jet with a secondary 
vertex and a soft lepton tag. The probability for all CDF data events to be due to a background 
fluctuation alone is 1 • 10^^, which is equivalent to a 4.8 a deviation in a Gaussian distribution. Again 
the top quark mass is kinematically reconstructed for W+ > 4 jets events as described above. The 
mass distribution is shown in figure [TSb. The best fit is obtained for Aftop — (176 ± 8 ± 10) GeV/c^. 

Simultaneously to CDF the D0 collaboration updated its top quark analyses based on data 
with an integrated luminosity of 50 pb~^ |207| . The updated analysis is very similar to the previous 
searches, involving the dilepton channel, soft muon b tagging and the topological analysis. From all 
channels, D0 observes 17 events with an expected background of 3.8 ± 0.6 events. The probability 
for this measurement to be an upward fluctuation of the background is 2 • 10~^, which corresponds 
to 4.6 standard deviations for a Gaussian probability distribution. To measure the top quark mass, 
lepton+4 jets events are subjected to a constrained kinematic fit. The resulting top quark mass 
distribution is shown in figure [TOl for the standard cuts (a) and looser selection requirements (b). 
A likelihood fit to the observed mass distribution yields a central value for the top quark mass of 
Mtop = 1991:2? (stat.) ± 22 (syst.) GeV/c^. 

In September 1995 CDF completed the series of publications establishing the top quark discovery 
with a complementary kinematic analysis without b tagging, that used the Et of the second and third 
highest Et jets to calculate a relative likelihood for each event to be top quark like or background 
like |2081 1209j . The probability for the observed data to be due to a background fluctuation was found 
to be 0.26%. 

Finally, 17 years after the discovery of the b quark its weak isospin partner, the top quark, 
was firmly established. In figure [201 we summarise the long lasting quest for the top quark. The 
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figure shows the lower experimental limits on the top quark mass (histogram) [1861 12101 12111 11851 
12121 11951 11981 12011 1203j , theoretical predictions mainly based on flavour symmetries within the quark 
mass matrix (dots) |182[ 12131 12141 12151 12161 1217' '218j and predictions based on electroweak precision 
measurements at LEP and SLC (squares) [219, 220* [Ml Si [2H as discussed in section O 
The direct measurements by CDF [20l [2051 ISH [104] and D0 [20Z1 IllSl [2261 HOS] are also included 
(triangles). The good agreement between the measured top quark mass and the prediction obtained 
from electroweak precision measurements constituted a major success of the Standard Model. 

6. tt cross section measurements 

Within the SM the tt cross section is calculated with a precision of about 15% [99l [91], see also 
section [5.1.41 The SM further predicts that the top quark decays to a boson and a h quark with 
a branching ratio close to 100%. Measuring the cross section in all possible channels tests both, 
production and decay mechanisms of the top quark. A significant deviation from the SM prediction 
would indicate either the presence of a new production mechanism, e.g. a heavy resonance decaying 
into tt pairs, or a novel decay mechanism, e.g. into supersymmetric particles. The ti cross section 
depends sensitively on the top quark mass. In the mass interval of 170 < A/top < 190GeV/c^ the 
cross section drops by roughly 0.2 pb for an increase of 1 GeV/c^ in Mtop. This theoretically predicted 
dependence can be exploited to turn a cross section measurement into an indirect determination of 
the top quark mass. A 15% measurement of the cross section is approximately equivalent to a 3% 
measurement of Mtop. One can also turn the argument around and use the measurements of the cross 
section and the top quark mass and test their compatibility with the theoretically predicted cross 
section and its mass dependence as indicated in figure [H in section [5. 1.41 

The tt cross section measurements are very fundamental to top quark physics at the Tevatron, 
since these analyses isolate data samples that are enriched in ti events and lay thereby the foundations 
for further investigations of top quark properties. The ti cross section has been measured in the 
dilepton, the lepton-plus-jets and the all hadronic channel. So far the ti tau modes evaded observation. 
In this section we discuss the experimental methods applied in ti cross section measurements in more 
detail. We highlight representative Tevatron analyses for different channels and analyses methods, 
either from CDF or D0. A comprehensive summary of all measurements and the combined cross 
section results is presented in section WM 

6.1. Dilepton channel 

As mentioned in section 14.41 the dilepton channel comprises ti final states with two high pt charged 
leptons (electrons or muons), and two b quark jets. This clean signature allows to select a ti sample 
with high purity, reaching a signal to background ratio between 1.5 and 3 depending on the analysis. 
In 2004 CDF published its first ti cross section measurement at ^/s = 1.96 TeV in the dilepton channel 
based on a data sample with an integrated luminosity of C-mt = (197 ± 12) pb""'^ [227] . i.e. about twice 
as much data as used in Run I. Two complementary analyses are performed. The first one requires 
both leptons to be specifically identified as either electrons or muons, while the second technique allows 
one of the leptons to be identified only as a high-px isolated track, thereby significantly increasing 
the lepton detection efficiency at cost of a moderate increase in the expected backgrounds. In the 
following we will concentrate the discussion on the second, the lepton-plus-track analysis. 

At trigger level the data samples used in dilepton analyses are selected by finding events that 
either contain a central electron or muon candidate with i?T > 18 GeV, or an electron candidate in 
the forward calorimeter with > 20 GeV and > 15 GeV. In the offline analysis two oppositely 
charged leptons with Et > 20 GeV are required. One lepton, the "tight" lepton has to pass strict 
lepton identification criteria and be isolated. Tight electrons have a well measured track pointing at 
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an energy deposition in the calorimeter. For electrons with \rj\ > 1.2, this track association is made by 
a calorimeter-seeded silicon tracking algorithm. Tight muons must have a well-measured track linked 
to a track segment in the muon chambers and an energy deposition in the calorimeters consistent 
with that expected for muons. The second lepton, the "loose" lepton, is defined as a well-measured, 
isolated track with px > 20 GeV/c and a pseudorapidity of \rj\ < 1.0. 

Candidate events must have a J^t > 25 GeV. The value of is corrected for all loose leptons 
if the associated calorimeter Et is less than 70% of the track px- Events are rejected if the vector 
-^T lies within 5° of the loose lepton axis and < 50 GeV. Jets are counted with a threshold of 
Et > 20 GeV and within the pseudorapidity range of \r]\ < 2.0. At least two jets defined in that 
way are required. Candidates being compatible with cosmic ray muons or photon conversions are 
removed. To remove dilepton pairs due to Z° boson production the cut on the missing transverse 
energy is tightened to > 40 GeV in a window of ±15 GeV/c^ around the Z° mass. 

The dominant backgrounds to ti dilepton events are Drell-Yan {qq Z/Y* ^ £^£^) production, 
misidentified (fake) leptons in — > iiy + iets events where one jet is falsely reconstructed as a lepton 
candidate, and diboson {WW, WZ, and ZZ) production. The Drell-Yan background is estimated by 
a combination of fully simulated Pythia Monte Carlo events and CDF data. A sample of Z boson 
candidates in the mass range of 76 - 106 GeV/c^ is selected. The numbers of events passing the 
nominal tt selection or a Drell-Yan specific selection are obtained, respectively. These two numbers 
provide the normalization for the expected Drell-Yan background contributions to the ti dilepton 
sample as obtained from Monte Carlo. The rate of lepton misidentification is obtained from an 
inclusive jet sample after removing sources of real leptons such as W and Z boson decays. The 
accuracy and robustness of the derived lepton fake rate is checked with several control samples and 
good agreement is found within the statistical uncertainties. The diboson background estimates are 
derived by calculating the acceptances from Monte Carlo and applying the theoretical cross sections 
which have a relatively small uncertainty. The total number of background events expected in the 
data sample is 6.9 ± 1.7 of which 61% are due to Drell-Yan processes, 22% are due to fake leptons, 
and 17% are due to diboson production. CDF observes 18 events in data. The tt cross section is 
calculated according to the formula 

^Gvt ' ■^int 

where iVobs is the number of observed events, A'bkg is the number of expected background events, 
and Ecvt is the event detection efficiency which includes the kinematic acceptance, trigger and 
reconstruction efficiencies, and the branching ratio into dilepton events. The event detection efficiency 
is Ecvt = (0.88 ± 0.12)%, including a branching ratio of BR(W^ £u) = 10.8% for the W boson decay 
into a charged lepton plus neutrino. The kinematic acceptance is evaluated for a top mass of Mtop = 
175GeV/c^. The resulting value of the cross section is 7.0^2 3 (stat.)l|^'3 (syst.) ± 0.4 (lumi.)pb, 
which is in very good agreement with the theoretical prediction of 6.7^Q gpb [91j . The distributions 
of various kinematic variables are found to be consistent with the SM prediction as obtained from 
Pythia Monte Carlo. As an example, the distribution of the scalar sum of all transverse energies in 
the event is shown in figure [2U A Kolmogorov-Smirnov test of this distribution yields a p-value of 
75%. 

6.2. Lepton-plus-jets channel 

The signature of the tt lepton-plus-jets channel comprises a high-pT electron or muon, missing 
transverse energy and four jets, see section 14.41 The branching fraction of this channel is about 
30% which is one of the advantages over the dilepton channel. However, W + multijet backgrounds 
are large and call for dedicated strategies to improve the signal to background ratio. In this section 
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Figure 21. The scalar sum of all 
transverse energies in the event. The 
histograms are stacked. The shaded 
histograms are the backgrounds, the 
hatched histogram shows the tt signal 
distribution scaled to the measured 
cross section. The dots plus error bars 
are data from the first CDF ti dilepton 
analysis in Run II |227| . 
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Figure 22. Scatter plots of aplanarity 
versus for (a) the W + 4 jets 

background sample, (b) QCD multijet 
background events, (c) tt Monte Carlo 
events from Herwig, and (d) D0 
lepton plus jets data. The dashed lines 
represent the threshold values used for 
the selection in the D0 topological 
analysis 105 . 



we discuss four important methods to reduce backgrounds in the lepton-plus-jcts channel and obtain 
a measurement of the ti cross section. 

6.2.1. Topological analyses The production of two very heavy objects in tt events is reflected in 
several kinematic distributions that describe the general event topology. Some analyses take advantage 
of the fact that the scalar sum of transverse energies, iJx, tends to be much higher in ti events than 
in background events. This is due to the fact that jets originating from the decay of a heavy object 
are typically much harder, i.e. higher in i?T, than those from gluon radiation. Other analyses utilise 
the higher spherical symmetry of ti events. More recent analyses combine several kinematic variables 
into one kinematic discriminant using a neural network [228" . 

We discuss here in more detail the final D0 topological analysis used to measure the ti cross 
section in Run I at = 1.8 TeV |105j . This analysis marks the final result of the technique used in 
D0 to discover the top quark in the lepton-plus-jets channel in 1995 [207j . The data were selected 
by a lepton-plus-jets trigger, requiring one loosely identified electron or muon and at least one jet. 
The electron data set corresponds to a luminosity of 119.5 pb^^, the muon data set to 107.7 pb^^. 
The offline selection asks for an isolated, well identified electron with > 20 GeV and \ri\ < 2.0 
or an isolated, well identified muon with > 20GeV/c and \ri\ < 1.7. The missing transverse 
energy as calculated in the calorimeter has to be above 25 GeV for electron-plus-jets events and above 
20 GeV for muon-plus-jets events. The pseudorapidity of the reconstructed W boson is required to 
be < 2.0. Jets are counted with a threshold of Et > 15 GeV and \r]\ < 2.0. Events with four 

or more jets are accepted. To avoid an overlap with the soft muon ti analysis events with jets having 
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a soft muon h tag are excluded. 

Three topological variables are used to further discriminate between signal and background: (1) 
the -STt variable which is defined as the scalar sum of all transverse jet energies using the jet definition 
as mentioned above, (2) the scalar sum of I^t and the lepton which is denoted El^^ and (3) the 
aplanarity A which is an event shape variable quantifying the "flatness" of an event and is defined as 
I times the smallest eigenvalue of the normalized laboratory momentum tensor M of the observed 
physics objects. The tensor M. is defined by: 

where is the three momentum of the object k and the indices i and j are the x, y and z coordinates. 
The sum runs over all objects under consideration, i.e. all jets defined by the cuts > 15 GcV and 
\t]\ < 2.0, and the W boson reconstructed from the lepton and ^t- Large values of A are indicative 
of spherical events, whereas small values correspond to more planar events. The following cuts on 
topological variables are used: Ht > 180 GeV, EJ^ > 60GeV and A > 0.065. These cuts result from 
an optimization procedure such that the smallest error on the measured cross section and the best 
signal to background ratio S/B are achieved. 

The main backgrounds to the topological ti analysis arise from M^+jets events and QCD multijet 
events that contain a misidentified electron or an isolated muon and mismeasured ^t- The background 
estimate proceeds in three major steps. In a first step the QCD multijet background is estimated as 
a function of jet multiplicity from data samples where all selection cuts except the three topological 
cuts are applied. The electron and the muon samples are treated separately because processes that 
give rise to a misidentified electron or an isolated muon are significantly different. 

Jets that produce showers with a large fraction of the energy in the electromagnetic calorimeter 
can sometimes pass the electron selection criteria and be falsely identified as an electron. To determine 
the rate of multijet events containing such misidentified electrons one examines the spectrum of 
events that pass the electron trigger but fail the full offline electron identification. These events 
correctly describe the shape of the distribution of QCD multijet events, while the normalization 
has to be found by matching the number of events at low J^t between this background sample and 
the sample where the full electron identification has been applied. The number of events in the tail of 
the normalized distribution above > 25 GeV provides an estimate on the number of QCD multijet 
background events. 

Muons from semileptonic 6 of c quark decays are normally within a jet, i.e. they are nonisolated. 
However, occasionally the decay kinematics is such, that there is insufficient hadronic energy to 
produce a jet. In this case the muons from heavy quark decays appear to be isolated and constitute 
the major source of QCD multijet background in the muon-plus-jets sample. The rate of isolated 
muons from heavy quark decays is measured using multijet samples with < 20 GeV. 

The multijet background is estimated for each inclusive W + jet multiplicity sample. In the 
e+ > 4 jets sample the background is estimated to be 7.2 ± 2.2 events, in the /i+ > 4 jets sample 
13.9 ± 4.4 events. 

In the second step of the background estimate the rate oi W + multijet events is computed 
using a fit to the jet multiplicity spectrum that remains after the subtraction of the QCD multijet 
background. This method exploits a simple exponential relationship between the number of events 
and the jet multiplicity 

ajW + n jets) ^ 

a{W +{n~l) jets) " ' ^ ' 

where a is a constant that depends on the specific jet definition, i.e. the clustering algorithm as well as 
the Et and rj requirements. Relation pS)) is known as "Berends scaling" |2291 1230] and based on the 
assumption that the observed jets come from gluon radiation where the emission of each additional 
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jet is suppressed by a factor of the strong coupling constant a^. A fit to the observed D0 data yields 
a scaling parameter of a = 0.18 ± 0.01 for electron data and a — 0.19 ± 0.02 for muon data. These 
fit values are used to predict the number of background events in the W-\- > 4 jets sample and are 
found to be 44.8 ± 8.6 events for electron-plus-jets events and 25.8 ± 4.6 events for the muon-plus-jets 
sample. As a side remark it has to be noted that a cannot simply be identified with as, since the 
variation of as in the PDFs has to be taken into account [231] . If this is properly done the sensitivity 
to as is largely lost. 

In the third step the cut efficiencies for the topological cuts are computed and then applied to 
the previously derived background estimate for the W+ > 4 jets sample. The efficiency on the QCD 
multijet background is again derived from data, while the efficiency on the W + jets background is 
taken from events generated with the Vecbos Monte Carlo program |177j . After applying all cuts 
the total background is predicted to be 4.51 ± 0.91 in the e+ > i jets sample and 4.32 ± 1.04 events 
in the iJ,+ > 4 jets sample, while 9 or 10 events are observed in data, respectively. 

The event detection efficiency Cevt for the tt signal is determined from Herwig Monte Carlo events 
that were passed through the D0 detector simulation. For a top quark mass of Aftop = 170GeV/c^ 
a value of Cevt = (1.29 ± 0.23)% is found for the electron sample, and Ccvt = (0.91 ± 0.27)% for the 
muon sample. Using eevt determined for a top mass of Mtop = 172.1 GeV/c^ as measured by D0 the 
cross section is computed to be (2.8 ± 2.1) pb in the electron-plus-jets channel and (5.6 ± 3.7) pb in the 
muon-plus-jets channel. Figure [22] shows the aplanarity versus Ht plane for the background samples, 
tt Monte Carlo and D0 data. The plots illustrate the effectiveness of reducing the background via 
the cuts on the topological variables Ht and A. 

6.2.2. Secondary vertex tag While each tt event features two h quark jets, only about 2% of the 
W -I- jets background contain a h quark. Therefore, the tt signal can be significantly enhanced by 
identifying h quark jets. Three different h jet identification methods are discussed in section [4.5.41 
In CDF the secondary vertex tagging method proved to be the most effective one, yielding the cross 
section measurement with the smallest total uncertainty. 

We present here a summary of the first CDF Run II tt cross section measurement using 
lepton-plus-jets events with secondary vertex b tagging |168| . The b tagging method has already 
been presented in section 14.5.41 That is why we concentrate here on the event selection and the 
background estimate. The analysis uses a data sample triggered by high momentum electrons or 
muons corresponding to an integrated luminosity of 162 pb^^. The electron selection requires an 
isolated cluster in the central calorimeter with > 20 GeV matched to a track with > 10 GeV/c. 
Muon candidates have a track in the drift chamber with pT > 20GeV/c that is matched to a track 
segment in the muon chambers. Events consistent with being photon conversions (electrons) or cosmic 
rays (muons) are rejected. The missing transverse energy is required to be > 20 GeV. By requiring 
one and only one well identified lepton tt dilepton events and Z — > e"'"e~/^^/i~ events are suppressed. 
To improve the removal efficiency for Z bosons, events are also removed if a second, less stringently 
identified lepton is found that forms an invariant mass Mu with the primary lepton within the window 
of 76 < Mu < 106 GeV/c^. Jets are defined as clusters in the hadronic calorimeter with Et > 15 GeV 
and |?7| < 2.0. At least three jets are required for an event to fall into the signal region. One of these 
jets has to be identified as containing a b quark using the secondary vertex tag algorithm. The final 
cut is on the total transverse energy and demands > 200 GeV, which rejects approximately 40% 
of the background while retaining 95% of tt signal events. 

The secondary vertex tag algorithm is described in section [4.5.4l Figure[T7k shows the distribution 
of the two-dimensional decay length of secondary vertices in the CDF data sample compared to the SM 
prediction of tt signal and background. Good agreement is found. The backgrounds to the secondary 
vertex tagged sample are (i) direct QCD production of heavy flavour quarks without an associated 
W boson (non-W^ QCD), (ii) W plus light quark jets events where one jet is falsely identified as 
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heavy flavour (mistags), (iii) W plus heavy flavour jets, (iv) diboson production, single top quark, 
and Z r+r^ production. The estimate of these backgrounds is partially derived from CDF data 
and partially from Monte Carlo. In particular, the number of b tagged + jets background events is 
calibrated with the number of observed W + jets events before b tagging. Therefore, the first step in 
the background calculation is to estimate the number of background events that do not contain a W 
boson in the pretag sample and subtract that background from the observed number of events. 

(i) Non-W QCD background: The non-W background is mainly due to events where a jet is 
misidentified as an electron and the is mismeasured, or due to muons from semileptonic b decays 
which pass the isolation criterion. Since the background sources in the electron and muon sample 
are different, one has to treat these two samples separately. The non-W background estimate uses 
the variable and the isolation variable i?iso , that is defined as the ratio of calorimeter energy i^iso 
contained in an isolation annulus of AR = 0.4 around the lepton (excluding the energy associated to 
the lepton) divided by the lepton energy Eg. The i?iso versus I^t plane is divided into a signal region 
(i?iso < 0.1 and > 20GeV) and three sideband regions. One assumes that the two variables are 
uncorrelated for non-H^ background events and calculates the number of background events in the 
signal region as a simple proportion of events in the sideband regions. The contribution of true W 
and tt events in the sideband regions is subtracted using Monte Carlo predictions normalized to the 
observed number of events in the signal region. In the pretag e+ > 3 jets sample (20 ± 5)% of the 
events are estimated to be non-W events, in the > 3 jets sample it is (7.5 ± 2.3)%. The pretag 
sample is defined as the data sample where all selection cuts have been applied except requiring a 6 
tag. In the final sample, after requiring a secondary vertex tag, about 18% of the total background is 
due to non-VF events. 

(ii) Mistags: The mistag rate of the secondary vertex algorithm is measured using inclusive 
jet samples. Mistags are caused mostly by a random combination of tracks which are displaced 
from the primary vertex due to tracking errors. The main idea is to use the rate of events with 
negative two-dimensional decay length as an estimate of the mistag rate. Corrections due to material 
interactions, long-lived light flavour particles (e.g. and A), and negatively tagged heavy flavour 
jets are determined using fits to the effective lifetime spectra of tagged vertices. The mistag rate is 
parameterized as a function of four jet variables: i?T, the good track multiplicity, ij and (j) of the jet 
as well as one event variable, i.e. the scalar sum of the Et of all jets with Ej- > 8 GeV. To estimate 
the mistag background in the W + jets sample each jet in the pretag sample is weighted with its 
mistag rate. The sum of weights over all jets in the sample is then scaled down by the fraction of 
non-Vl^ events in the pretag sample. Since the mistag rate per jet is sufficiently low, this prediction of 
mistagged jets is a good estimate on the number of events with a mistagged jet. It is found that 34% 
of the background for the final tt selection are due to mistags. 

(iii) W + heavy flavour: The production of a boson in association with heavy fiavour quarks 
is the main background to the tt signal with a secondary vertex tag. Heavy quarks occur in the 
process qiq2 ^ W + g where the gluon splits into a 66 or a cc pair, and in the process gq — > Wc. As 
mentioned in section HTBl the Alpgen Monte Carlo program [178] is used to generate several samples 
of exclusive W + n jets final states. This includes W + bb/cc + n jets, and Wc -I- jets. The Alpgen 
generation is followed by showering with Herwig. To reproduce the entire W + multijets spectrum 
the exclusive samples are merged together taking the appropriate cross sections into account. The 
combination procedure avoids double counting of jets produced at matrix element level with Alpgen 
and hard jets produced by the showering in Herwig. A jet clustering algorithm is run at Monte Carlo 
particle level after showering but before detector simulation. These Monte Carlo jets are matched to 
partons generated at matrix element level. If extra jets occur, that do not match a parton, the event 
is rejected. 

While the shape of the W + n jets spectrum can be reproduced well by the Monte Carlo, the 
absolute normalization has large theoretical uncertainties and is therefore taken from collider data. 




Figure 23. (a) Pseudo-cr distribution for inclusive jet data. The distribution is used to measure 
the heavy flavour fractions in the jet sample. The fitted contributions for the different heavy flavour 
components and secondary interactions in light flavour jets are shown |168l . (b) Jet multiplicity 
spectrum for CDF data (dots), the total background expectation (dashed histogram), and the sum 
of background and the predicted tt signal (cr = 6.7pb, full line). 



The heavy flavour fractions predicted by Alpgen are found to be too small and are calibrated against 
fractions measured in inclusive jet data. The composition of the inclusive jet sample is determined 
by fitting the pseudo-cr distribution for b tagged jets. The pseudo-cr is defined as L2D • M^ik/Pt^, 
where L2D is the two-dimensional decay length (see also section 14.5.41) , Mvtx is the invariant mass 
of all tracks associated to the secondary vertex, and p^^^ is the transverse momentum of the vertex 
4- vector. Figure [23b shows the pseudo-cr distribution for the inclusive jet data as well as the fitted 
heavy flavour components. The template distributions for heavy flavour jets used in the fit are taken 
from Alpgen Monte Carlo samples. The measured heavy flavour fractions in jet data are consistently 
higher by 50% than the Alpgen prediction. Therefore, a correction factor of 1.5 ±0.4 is applied to the 
heavy flavour fractions in the W^-|-multijets sample. The assumption here is that this correction factor 
is universal and can be transferred from the inclusive jets to the W + multijets sample. As a result, 
it is found, that the fraction of -I- 4 jets events coming from the Wbb process is (4.8 ± 1.3)%, the 
fraction with one or two c jets from M^cc is (7.3±2.0)%, and the fraction of VFc events is (6.1±1.3)%. 
It has to be noted that the Wc fraction is directly taken from the Alpgen prediction and is not 
corrected, since it is due to a different physics process. 

The number of W+ heavy flavour background events for the ti analysis is computed by multiplying 
the derived heavy quark fractions by the number of pretag events, after subtracting the non-VF 
background. The heavy flavour contributions to the total background in the tt candidate sample are 
as follows: 25% are due to Wbb, 8% due to Wcc and 6% due to Wc production. 

(iv) Diboson, single top and Z t^t^ : There is a number of smaller backgrounds which can 
be reliably predicted by combining the event detection efficiency for these processes as determined 
from Monte Carlo events with the theoretically predicted cross section. This method is feasible, since 
the cross section predictions for diboson production processes, i.e. WW, WZ and ZZ, and single 
top quark production have relatively small uncertainties. Diboson events can mimic a tt signal if one 
boson decays leptonically and the other one decays into jets, where at least one jet is due to a 6 or c 
quark. In addition, Z -I- jet production can mimic tt events if the Z boson decays into r+r^ and one 
r decays leptonically producing an isolated electron or muon, while the second r decays hadronically. 
The contribution of diboson, Z — + r+r~ and single top events to the total background is found to be 
9%. 



Top quark physics in hadron collisions 



55 



The total background prediction is 13.5 ± 1.8 events, while 48 events are observed in the W+ > 3 
jets sample. The clear excess of events is attributed to tt production. Figure [^5b shows the jet 
multiplicity spectrum for CDF data, the total background expectation, and the sum of background 
and the predicted tt signal. A very good agreement between the data and the SM expectation is 
found. 

To compute the tt cross section according to (f26| , the event detection efficiency ecvt for tt events 
is needed. It is obtained from a sample of Pythia tt Monte Carlo events. The trigger efficiency as 
well as several correction factors that account for differences between data and Monte Carlo events 
are measured with control data samples. Particularly important is the correction factor for the b 
tagging efficiency. A sample of inclusive electrons with pT > 7 GeV is used for this purpose. Many 
electrons in this momentum regime originate from semileptonic b decays. Therefore, the sample is 
enriched in heavy flavour. Using a double tag technique it is found that the average tagging efficiency 
of a 6 quark jet in data is (24 ± 1)%, while the Monte Carlo predicts (29 ± 1)%, thus yielding a 
b tagging correction factor of 0.82 ± 0.06, where the error includes systematic uncertainties. The 
efficiency to tag at least one jet in a ti event, after all other cuts have been applied, is found to be 
(53.4 ± 3.2)%. The uncertainty is almost purely systematic. The overall event detection efficiency is 
Cevt = (3.84 ±0.40)%. This yields a cross section measurement of a{tt) — 5.6li['^ (stat.) 1q g (syst.)pb, 
which is in good agreement with the theoretical expectation of 6.7^Q gpb. The statistical uncertainty 
is larger than the systematic one, but very soon, using more data, the systematic uncertainty will start 
to dominate the measurement. Although, having larger control samples at hand, it will be possible 
to also reduce the systematic uncertainty. 

6.2.3. Impact parameter tag The impact parameter b tagging algorithm is an alternative method 
to identify b quark jets. The method also relies on the long lifetime of b hadrons, but the explicit 
reconstruction of a secondary vertex is not required. The algorithm rather asks for tracks with high 
impact parameter significance and is explained in section [4.5.4l In the CDF Run II analysis presented 
here a jet is b tagged if its probability to be consistent with a zero lifetime hypothesis is below 
0.01 [232j . The b tagging efficiency of the algorithm is measured in inclusive electron data using a 
double tag method. For heavy flavour jets with i?T > 10 GeV the average tagging efficiency is found 
to be (19.7 ± 1.2)%. 

The event selection for the tt cross section measurement is essentially the same as described in 
section (622] for the analysis using a secondary vertex b tag. There is one additional cut which requires 
that the between the and the most energetic jet is within [0.5, 2.5] radians if If!^ < 30 GeV. The 
cut on Ht is omitted. The background is estimated based on the same techniques as the measurement 
using the secondary vertex tag. The mistag rate of light quark jets is measured in an inclusive 
jet data sample and found to be (1.11 ± 0.06)%. In a data sample corresponding to 162 pb~^ 59 
W+ > 3 jets events are observed after all cuts, while the total background is estimated to be 24.7l4'g 
events. The overall event detection efficiency is found to be Covt = (4.09 ± 0.61)% which includes 
the event b tagging efficiency of (57.2 ± 3.9)%. As a result, the tt cross section is calculated to be 
a{tt) = 5.8tli (stat.) ± 1.3 (syst.) pb, which is in good agreement with the result obtained with the 
secondary vertex b tagging algorithm. 

6.2.4. Soft lepton tag The third technique to identify b quark jets searches for electrons or muons 
within jets originating from semileptonic b decays. We present here a CDF analysis using a soft muon 
tag |233j . A data sample with an integrated luminosity of 194 pb~^ is used. The event selection before 
b tagging is the same as for the secondary vertex tag analysis. After all cuts 337 events are retained. 

The soft muon tag algorithm uses a global built from several characteristic distributions that 
separate muon candidates from background. A jet is considered as b tagged if it contains a muon with 
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Pt > 3GeV/c and < 3.5 within Ai? < 0.6 of the jet axis. Events are rejected if the isolated high 
Pt lepton is a muon of opposite charge to the soft muon and the dimuon invariant mass is consistent 
with a J/iJ, an T or a The efficiency to tag a tt event is approximately 15%, with small variations 
between the + 3 jets and W + A jets samples. The rate of background events where a light quark 
jet is misidentified as containing a semileptonic b hadron decay is estimated from a sample of photon- 
plus-jets events. The misidentification probability of a track with pr > 3GeV/c and Ai? < 0.6 is 
found to be about 0.7%. 

In the W+ > 3 jets data sample 20 events with a soft muon tag are found, while the total expected 
background is estimated to be 9.5±1.1 events. The total event detection efficiency is eevt = (1.0±G.1)% 
and the resulting tt cross section is found to be a{tt) — 5.3 ± -1-3.3 (stat.) ti^o (syst.) pb. 

6.3. All hadronic channel 

While the all hadronic channel has the largest branching ratio of all tt event categories, the QCD 
multijet background is overwhelming and it is therefore very challenging to isolate a tt signal. We 
present here a D0 analysis based on a Run II dataset corresponding to 162 pb~^ [234j . To enrich the 
signal content secondary vertex b tagging and a multivariate analysis with artificial neural networks 
is employed. 

The data sample is collected with a dedicated jet trigger at an average trigger efficiency of 74%. 
The event pre-selection asks for six or more jets, that are reconstructed with an algorithm using a 
fixed cone radius of Ai? = 0.5. Jets are counted with Et > 15 GeV and \ri\ < 2.5. Events that contain 
an isolated lepton are rejected to obtain a dataset orthogonal to the lepton-plus-jets analysis. The 
reconstruction of secondary vertices is used to identify b quark jets. A jet is considered to be tagged 
if a vertex with signed decay length significance larger than 7.0 is found. Exactly one b tagged jet is 
required. Double tagged events are rejected to ease the background estimation. The efficiency to tag 
exactly one jet in a hadronic tt event is determined to be 46%. 

The multivariate analysis combines 13 kinematic or event shape variables using two neural 
networks. The background sample for training the networks is obtained from the pretag data sample. 
A parametrization of the tagging rate is used to select events that have high probability to be tagged, 
but do in fact not have a tag. This ensures on one hand that the background sample is similar to the 
tagged events, but on the other hand is disjoint to the candidate sample. Signal Monte Carlo events 
are generated with Alpgen in combination with Pythia for showering and hadronization. After all 
cuts 220 candidate events are observed over a total background of 186 ±5 events. The event detection 
efficiency is Cevt = (2.8 ± 0.8)% including a hadronic branching ratio of 46.19%. The tt production 
cross section is measured to be (7{tt) = 7.7^3 ^ (stat.)l3 g (syst.) ± 0.5 (lumi.) pb. 

6.4- Cross section combination 

In this section we summarise the final Run I tt cross section measurements and give a brief overview on 
the status of the currently available Run II measurements. The measured cross sections are presented 
in table m For Run I the cross sections obtained from combination of all tt event categories for CDF 
and D0 are also shown [103) 1105] . The combined values should be compared to the predicted cross 
section. Within the uncertainties good agreement is found. 

In Run II various new techniques to measure the tt cross section are introduced. While the 
traditional dilepton analysis uses two leptons that meet strict identification criteria, CDF added an 
analysis where the identification of the second lepton is relaxed to increase the acceptance and reduce 
the statistical uncertainty. This analysis is discussed in more detail in section 16.11 Table [8] quotes 
the combined value for both CDF dilepton analyses |227j . A third CDF dilepton analysis uses a 
kinematic fit to the (_^t, Afjct) phase space to disentangle the major SM processes contributing to the 
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Table 8. Measurements of the tt cross section at the Tevatron. The CDF measurements 
assume a top quark mass of Mtop = 175GeV/c'^ 11031 12351 12361 . the D0 measurements use 
Mtop = 172.1 GeV/c^ |105) . In Run II a top mass of Mtop = llbGeV/c^ is assumed for the 
measurements of both experiments. The D0 Run I lepton-plus-jets cross section is a combination of 
the topological analysis and the soft lepton tag analysis. The given uncertainties include statistical 
and systematic contributions. 
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Channel 


CDF D0 


CDF 
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l.Qtll pb 
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5.8 ± 1.8 pb 
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6.6 ± 1.9 pb 




All hadronic 


l.Q+Jd pb 


7.3 ±3.2 pb 


l.b+Jdi pb 


l.l+dl pb 


Combined 


6.5li:Ipb 


5.7 ± 1.6 pb 






Predicted [91] 


5.2l 


pb 


6-7±E!:^ pb 



dilepton sample (global kinematic fit) [237| . D0 has made cross section measurements in a dilepton 
data sample without b tagging |238j and in a e/i sample using a secondary vertex tag |239j . 

The CDF lepton-plus-jets analyses using three different b tagging algorithms are presented in 
section 1121 Based on a data set corresponding to 230 pb~-^ the D0 collaboration has also presented 
measurements using a secondary vertex b tag and an impact parameter tag [240j . CDF has performed 
two analyses that exploit kinematic properties of tt lepton-plus-jets events. One analysis uses the 
sample with secondary vertex tags and adds a likelihood fit to the transverse energy of the leading 
jet [241j . The second analysis is based on the W+ > 3 jets sample before b tagging and uses a neural 
network to construct a powerful discriminant between the backgrounds and the tt signal [228] . 

The D0 measurement of the tt cross section in the all hadronic channel [234] is discussed in 
section 16.31 The CDF Run II analysis in the all hadronic channel is similar to the D0 one, also 
requiring a secondary vertex b tag and exploiting several kinematic variables [2421 . But CDF uses 
simple sequential cuts rather than a multivariate combination via a neural network. At present, 
combined results on the tt cross section in Run II are not yet available. The combination of the 
different tt event categories (dilepton, lepton-plus-jets, all hadronic) is relatively straight forward, 
since the data sets are essentially orthogonal. However, the different analyses within a particular 
category are correlated, sometimes even strongly correlated. The correlation has to be determined 
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and taken into account when a combination of the different analyses is done. 
7. Top quark mass measurements 

As discussed in section I2.3[ the top quark mass is a very important parameter for the description of 
electroweak processes and precision tests of the SM. The main impUcation is that the SM relates the 
masses of the top quark, the W boson and the Higgs boson, such that precise measurements of Mtop 
and Mw imply a prediction for Mh- Therefore, the precise measurement of the top quark mass is one 
of the major goals of Run II at the Tcvatron. However, since Run II results on the top quark mass are 
not yet published, we describe in this chapter the two main methods used by CDF and D0 in Run I 
to determine Mtop: (1) the template method of CDF, which fits template distributions for different 
top quark masses to the distribution observed in data, and (2) the matrix element method of D0, 
which exploits the sensitivity of the leading order matrix element for tt production to Afjop. Both 
methods are also used in Run II to measure Aftop- In section [7.31 we discuss the combination of top 
quark mass measurements in different tt event categories and among the two Tevatron experiments. 

7.1. Template method in tt lepton-plus-jets events 

The best measurement of the top quark mass is achieved in the ti lepton-plus-jets channel, since it 
features a relatively high number of candidate events, moderate background levels and allows for a 
full reconstruction of top quark momenta with reasonable accuracy. That is why we present here the 
final top quark mass measurement in lepton-plus-jets events at CDF in Run I |104| . The data sample 
used in the analysis corresponds to an integrated luminosity of 106 pb~^. The event pre-selection is 
the same as the one in the Run I tt cross section analysis [103) , which is essentially the same as the 
Run II selection described in section 16.2.21 For the full reconstruction of tt candidate events a W 
candidate decaying into ei' or iiv and at least four jets are needed. Therefore, events in the -I- 3 
jets sample are used only if they feature at least one additional jet with Et > SGeV and {rjl < 2.4. 
After pre-selection, but before requiring a b tag the candidate sample consists of 163 events. Four 
subsamples are used to measure the top quark mass: (1) The first subsample contains events with 
two secondary vertex b tagged jets. (2) The second subsample consists of events with exactly one 
secondary vertex tag. (3) The third subsample includes events with one or two soft lepton tags, but 
no secondary vertex tag. (4) The fourth subsample contains events with no b tag and at least four 
jets with Et > 15 GeV and |?7| < 2.0. One very important issue for the top mass measurement is the 
reconstruction of jet energies, which is therefore detailed in the next section. 

7.1.1. Jet energy corrections The reconstruction of the top quark momenta in a ti candidate event 
is based on the assumption that the four leading jets can be identified with the primary partons from 
the top quark decay. To ensure the validity of this assumption several corrections have to be applied 
to the raw jet energies as measured in the calorimeter. While the general concept of jet reconstruction 
is briefly introduced in section [4. 5. 3i we summarise here the jet corrections applied in the CDF top 
mass analysis. 

(i) The relative energy correction accounts for non-uniformities in the calorimeter response as a 
function of 77 and is derived from dijet data, where one jet is measured in the central calorimeter 
(0.2 < |?7| < 0.7) and the second jet is in the forward region (1.1 < \ri\ < 4.2). The forward 
calorimeter is thus calibrated against the response in the central region. The uncertainty on the 
relative correction varies between 0.2% and 4%. 

(ii) Corrections for multiple interactions. A fixed amount of energy AE^^it ~ 0.297 GeV/c is 
subtracted from the jet Et for each additional reconstructed primary vertex in the event. 
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(iii) The absolute energy scale is a multiplicative factor that converts the energy observed in the 
jet cone into the average true jet energy. To derive the absolute correction the calorimeter 
simulation is tuned to agree with testbeam data of single electrons and pions. In a second step 
several parameters controlling the fragmentation process in the ISAJET Monte Carlo generator are 
tuned such, that the multiplicity distribution, the momentum spectrum, and the invariant mass 
distribution of charge particles, as well as the ratio of charged to neutral energy agree between 
ISAJET and dijet data. The absolute correction accounts for the nonlinearity of the calorimeter 
and energy losses near the boundaries of detector modules. The systematic uncertainty on the 
absolute energy scale is about 3%. 

(iv) Underlying event corrections account for extra energy contained in the jet cone due to particles 
coming from the fragmentation of partons that do not participate in the hard scattering of the 
primary pp interaction. A fixed amount of 0.65 GeV is subtracted from the Et of each jet . This 
correction is derived from minimum bias data. 

(v) Out- of- cone corrections account for the energy that is physicswise associated to the jet but falling 

out of the jet cone. The out-of-cone energy is related to low energy gluons emitted from initial 
partons, also referred to as soft gluon radiation. The correction is derived from Monte Carlo 
events. 

(vi) ti specific corrections are applied to the four leading jets and map their momenta to the momenta 
of the quarks from the ti decay. The corrections account for three effects: (a) The difference in 
the jet i?T spectrum of ti induced jets and the flat spectrum assumed in the previous corrections, 
(b) The energy loss in semileptonic h and c hadron decays, where the undetected neutrinos or 
muons that loose only little of their energy in the calorimeter, carry away part of the energy of the 
primary parton. (c) The multijet structure of ti events as compared to dijet events used to derive 
the other corrections. While the jet corrections (i) through (v) arc flavour independent the ti 
specific corrections treat h jets differently than light quark jets. Four jet types are distinguished: 
(1) jets used to reconstruct the hadronic W decay, (2) jets assigned to the h quark from the top 
quark decay without h tag or jets with secondary vertex tag, (3) jets with a soft electron tag, (4) 
jets with a soft muon tag. 

Within the range from 30 to 90 GeV the flavour independent corrections (number i through v) 
amount to an average correction factor of about 1.45. The total systematic uncertainty including 
all corrections varies between 7% for jets with corrected E^ of 20 GeV and 3.5% for jets with 
Et = 150 GeV. 

7.1.2. Event-by- event top mass fitting For top mass fitting the four-momenta of the particles of the 
ti decay chain arc fully reconstructed. The four leading jets are assigned to the four quarks in the final 
state of the hard scattering. The mass of light quark jets is assumed to be 0.5 GeV/c^, the mass of b 
quarks is set to 5.0GeV/c^. The jet-parton assignment bears some ambiguity. If none of the jets is 
tagged as a 6 jet candidate, there are 12 possible jet permutations. Including the two-fold ambiguity 
of the neutrino pz reconstruction, there are 24 combinations. With one b tag that number is reduced 
to 12, with two b tags there are only four possible assignments. For each event all of these possible 
combinations are tested. A kinematic flt based on a criterion is used to find the best combination 
and determine the top quark mass on an event-by-event basis. The expression implements six 
effective kinematic constraints: (1,2) the two transverse momentum components of the ti+ X system 
must be zero, (3) the invariant mass of the lepton and neutrino, Me^, must be equal to Mw, (4) the 
invariant mass of the two light quarks, Mjj, must be equal to Mw, (5,6) the two three-body invariant 
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Table 9. Number of expected and observed number of events in the four top mass subsamples. 



Data sample 


Expected Background 


Expected Signal 


Observation 


SVX double tag 


0.2 


6.1 


5 


SVX single tag 


2.7 


14.4 


15 


Soft Lepton Tag 


5.0 


4.0 


14 


No Tag 


32.4 


11.4 


42 


Total 


40.3 


35.9 


76 



masses, Af^^^- and Mjjj, must be equal to Mtop. The expression is given by 



X 



E 



£,jcts '■^ i'—x-,y 
(Mj-j- - Mw f ^ {M,,j - Mtop f ^ {M,,, ~ Mtop)^ ^ ^29) 

'^Mw ^Mtap ^Mtap 

The first sum in ([29]) is over the primary lepton £ and all jets with raw > 8 GeV and \ri\ < 2.4. 
The second sum is over the transverse components of the unclustered energy Ut, which is defined as 
the vector sum of the energies in the calorimeter towers after excluding the primary lepton and all 
jets with raw Et > 8 GeV and < 2.4. The symbols with a hat represent quantities which are free 
to be altered in the fit procedure. The uncertainty on a quantity X is denoted ax and occurs in the 
denominator of the respective terms. (Tmw is set to 2.1 GeV/c^, crA/t„p is set to 2.5 GeV/c^. For 
each combination of primary physics objects in an event the expression ([29]) is minimized. The 
combination with the lowest is chosen to be the best fit for that event. Events with their lowest 
above 10 are rejected. The fit yields an estimate of Mtop for each event. 
The mass fit performance is tested with ti Monte Carlo events. The best results are obtained for 
events where all four leading jets are correctly assigned to the appropriate quark. A mass resolution 
of 13 GeV/c^ is reached for these events. For events with two secondary vertex tags the correct 
assignment is made in 49% of the events, while for events without b tag the assignment is correct only 
in 23% of the cases. If the four leading partons from the tt decay cannot be uniquely matched to the 
four leading jets within AR < 0.4, the mass resolution deteriorates to an average of 34 GeV/c^. 

7.1.3. Final selection and backgrounds The final cut in the event selection requires that the of the 
best fit for an event is below 10. The background is estimated using the same methods as discussed 
for the tt cross section measurement in section |6^ In addition, the background is renormalized using 
a maximum likelihood fit to the observed rates of events with secondary vertex tag and soft lepton tag 
and their respective expectations. In the fit the sum of tt signal and background events is constrained 
to be equal to the observed number of events. The number of expected and observed events in the 
four top mass subsamples are given in table [5) 

7.1.4- Top mass determination In a last step the best estimate of the top quark mass is determined 
by a maximum likelihood fit to the distribution of reconstructed invariant masses. The shape of 
the reconstructed mass distribution for the tt signal is obtained from Herwig Monte Carlo events. 
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Figure 24. Top quark mass measurement with tlie template metiiod at CDF in Run I 11041 . (a) 
Top mass distributions of signal Monte Carlo for six different values of Mtop (dots). The overlayed 
curves are the fitted template functions for that particular top quark mass, (b) The histogram shows 
the reconstructed top mass distribution for all subsamples combined. The shaded region indicates 
the fitted signal component, the hatched region the fitted background. The inset shows the log 
likelihood as a function of the top quark mass. 



Several samples for different top quark masses are generated. The reconstructed mass distributions are 
referred to as templates, which provides the name for the entire measurement method. The template 
distributions are parametrized using a function that depends on 12 parameters and the top quark 
mass. The parameters are determined by a simultaneous fit to all template histograms. For each 
of the four subsamples a different set of parameters is calculated. Figure [24k shows the top mass 
template histograms and the fitted template function for six different values of Mtop. The Vecbos 
Monte Carlo program is used to create templates for the background. The background distribution is 
parametrized by a function with fewer parameters and no dependence on Mtop. 

Using the template functions and the observed mass distribution a likelihood function is defined. 
The template parameters and the background fraction are constrained to their central values within 
their uncertainties. The only parameter which is entirely unconstrained is Mtop. The log likelihood 
function is minimized with respect to all parameters in a combined fit to all subsamples. The result 
of the minimization for the combined sample is shown in figure [Mb . The minimum of the likelihood 
is reached for a top quark mass of Mtop — 176.1^5^0 GeV/c^. The given uncertainties are only 
statistical. Systematic uncertainties due to various sources are evaluated. The biggest uncertainty 
originates from the jet energy measurement (4.4GeV/c^). The modelling of initial and final state 
radiation causes an uncertainty of 2.6 GeV/c^. The modelling of the shape of the background spectrum 
induces an uncertainty of 1.3GeV/c^. Smaller sources of uncertainties are: the h tagging efficiency 
(0.4GeV/c^), parton distribution functions (0.3GeV/c^), and Monte Carlo generators (O.lGeV/c^). 
The contributions of all effects are added in quadrature resulting in a total systematic uncertainty of 
5.3GeV/c2. 
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7.2. Matrix element method 

Originally D0 also used a template method to measure the top quark mass. This method used a 
kinematically fitted mass, similar to CDF, and in addition discriminants based on a likelihood ratio 
or a neural network [243j . In 2004 D0 reanalyzed its top mass data sample with a new technique 
based on a matrix element method, which we will describe below 106 . 

The new measurement is based on the same data set of 91 lepton-plus-jets events as the previous 
one. The basic selection cuts are very similar to those used for the cross section measurement [244]. 
The data set corresponds to an integrated luminosity of 125 pb^^. The result of the template analysis 
is Mtop = 173.3 ± 5.6 (stat.) ± 5.5 (syst.) GeV/c^. 

The matrix element method is designed to extract more kinematic information from the events 
than the previous template methods, and thus yield an improved precision. The basic idea of the new 
method is to exploit the fact that the differential cross section for tt production depends sensitively 
on the top quark mass. The differential cross section can thus be used to calculate a probability for 
a certain top mass hypothesis. While the differential cross section consists of a phase space term and 
the matrix element for tt production, the matrix element is the more significant part, which motivates 
the name of the method. Matrix element methods have already been used previously to analyse tt 
dilepton events in CDF |245) and D0 ^246. J47j I226j . Since leading order matrix elements are used 
to calculate the event weights, only events with exactly four jets are used. This jet cut minimizes the 
effect of higher-order corrections and reduces the number of events from 91 to 71. 

The production probability P(a;, Mtop) for a tt event with a measured set of variables a; at a 
certain top quark mass Mtop is given as a convolution of the differential cross section da/dy with the 
parton distribution functions / and the transfer function W{y, x) that maps the measured quantities 
X into the quantities y at parton level: 



The parton distribution functions /((?;) are evaluated for the incoming partons with momentum 
fraction qi. The integral in (|30p is properly normalized by dividing by the total cross section (Jtotai- 
The integration runs over fifteen sharply measured variables, which are directly assigned to the parton 
quantities without invoking a transfer function. These variables are the eight jet angles, the three- 
momentum of the lepton, and four equations of energy-momentum conservation. The jet energies 
are not well measured and a transfer function Wjots(£'part, ^-jot) is needed to map jet energies Ejct 
measured in the detector to parton level energies -Epart- The function Wjots(£^part, ^jct) is a product 
of four functions F{Ep^j.^^, £jct)j one for each jet in the event. The functional form of F is the sum of 
two Gaussians. The parameters of F used for b quarks are different from those for light quark jets. 
After the first integration step there are five integrals left. One integral runs over the energy of one 
of the quarks from the hadronic W decay, the other four are over the masses squared, Mf, of the two 
W bosons and the two reconstructed top quarks in the event. This is an economical choice to safe 
computing time. When computing P{x, Mtop) all possible 24 permutations of jet assignments and the 
neutrino p^ solution are considered and the average is computed. 

Since the candidate sample still contains a considerable amount of background events, it is useful 
to calculate a background probability -Pbkg(2;) based on the W + 4 jets matrix element from Vecbos. 
Figure [^5b shows the Pbkg distribution for the sample of 71 candidate events and an overlay of tt 
and W + 4 jets Monte Carlo events. To reduce the bias in the top mass measurement due to the 
background content in the sample D0 employs a cut requiring log(Pbkg) < —11. After this cut the 
final top mass sample contains 22 events. The final likelihood for the measurement of Mtop is given 




(30) 
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Figure 25. (a) Distribution of the logarithm of the background probabiUty Pbkg for 71 candidate 
events in the D0 top mass sample [243]. To increase the purity of the sample a cut at log(Pbkg) < 
— 11 is performed, as indicated by the vertical line. The data are compared to ft signal Monte Carlo 
(left hatched) and W + A jets background Monte Carlo events (right hatched), (b) Negative of the 
log likelihood as a function of Aftop. (c) The likelihood normalized to its maximum value. The curve 
is fit to a Gaussian. The hatched area corresponds to the 68.27% probability interval. 



by 

N 

- InL(Mtop) = - ^ \n[ciPu{xi,Mtop) + C2Pbkg(a;0] + 

2 = 1 

Nci J A{x)Pu{x,Mt,p)dx + Nc2 J A{x)P^,^^{x) dx . (31) 

The integrals are calculated using Monte Carlo methods. The acceptance function A{x) is 1.0 
or 0.0, depending on whether the event is accepted or rejected by the analysis criteria. The 
sum runs over all N = 22 candidate events. The best value of Mtop and the parameters Ci are 
defined by minimizing — InL(Mtop), which is shown in figure [25b. Figure [25h shows the likelihood 
normalized to its maximum value. The Gaussian fit to the likelihood yields a top quark mass of 
-^top = (179.6 ±3.6) GeV/c^. Monte Carlo studies show that the extracted top quark mass has to be 
corrected by a shift of (SAftop = +0.5 GeV/c^. The total systematic uncertainty of the measurement 
is 3.9 GeV/c^ which is dominated by the uncertainty of the jet energy scale (3.3 GeV/c^). The final 
result is Mtop = 180.1 ± 3.6 (stat.) ± 3.9 (syst.) GeV/c^, which is a substantial improvement over the 
previous template result. The higher precision is mainly due to two differences in the analyses: (a) 
Well measured events contribute more than poorly measured ones, since a probability is assigned to 
each event, (b) all possible permutations of the reconstruction are included. Thus, the correct solution 
always contributes. 



7.3. Top quark mass combination 

The best precision on the top quark mass is obtained if the individual measurements in the different 
tt channels are combined for each experiment. Further improvement is achieved if the results of the 
two Tevatron experiments, CDF and D0, are combined. In the all hadronic channel CDF uses a 
data sample of 136 events with an estimated background of 108 ± 9 events |236| . A similar template 
method as described in section [7TT] gives a result of 186 ± 10 (stat.) ± 12 (syst.) GeV/c^. Since the 
dilepton topology comprises two unobserved neutrinos a straightforward full reconstruction of the 
event is not possible. To get an estimate of Mtop on an event-by-event basis a weighting method 
is used. For each event a weight distribution is calculated as a function of Mtop. The assigned 
weight depends on the agreement of the sum of the assumed transverse neutrino momenta with 
the observed missing transverse energy. The result in the dilepton channel at CDF is Mtop = 
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167.4 ± 10.3 (stat.) ± 4.8 (syst.) GeV/c^ |248j . The statistical uncertainties of the measurements 
in the three ti channels are uncorrelated, since the samples are statistically independent. However, 
the systematic uncertainties are correlated. For simplicity, the correlation is either assumed to be 
100% or zero. The uncertainties concerning the jet energy scale, the signal model (modelling of ISR 
and FSR, PDFs and h tagging), and the Monte Carlo generators are set to 100%. The correlation of 
uncertainties due to Monte Carlo statistics and the background model are assumed to be zero. The 
combination procedure uses a generalized method with full covariance matrix and yields a value 
of Mtop = (176.1 ± 6.6) GeV/c^ [104] . 

In the all hadronic channel D0 has measured Mtop = 178.5±13.7 (stat.)±7.7 (syst.) GeV/c^ [215] . 
In the dilepton channel D0 obtains Mfp = 168.4 ± 12.3 (stat.) ± 3.6 (syst.) GeV/c^ [250]. The 
combination of all topologies is entirely dominated by the lepton-plus-jets result described in detail in 
section [L2] and yields Mtop = (179.0 ±5.1) GeV/c^. Combining all CDF and D0 Run I measurements 
the top quark mass is determined to be 

Mtop = (178.0 ±4.3) GeV/c^ [1]. 

7. ^. Preliminary Run II results 

First preliminary Run II measurements of the top quark mass are now available and prepared for 
publication. CDF has been exploring several different techniques for the measurement of Mtop. The 
best single measurement is obtained from an extended mass template technique in the lepton-plus-jets 
channel. The analysis uses a data sample corresponding to 318 pb^^. The measured invariant mass 
of the hadronic W boson decay is used to reduce the systematic uncertainty on the jet energy scale. 
In a two-dimensional likelihood fit the top quark mass and the jet energy scale (JES) are obtained 
simultaneously. CDF measures Mtop = 173.5l2:6 (stat.) ± 2.5 (JES) ± 1.7 (syst.) GeV/c^. 

D0 has presented a preliminary result which is also based on a mass template method. In b tagged 
lepton-plus-jets events the top quark mass is found to be Aftop = 170.6±4.2 (stat.)±6.0 (syst.) GeV/c^. 
The analysis is based on a data sample of 229 pb~^ and uses 69 candidate events [251] . 

7.5. Future prospects 

The aim of Run II at the Tevatron is to reach a total uncertainty of 2 to 3 GeV/c^ on Mtop. At 
the LHC the top mass will be measured with negligible statistical uncertainty, while the systematic 
uncertainty is predicted to be on the order of 1 GeV/c^ towards the end of the running period [5]. In 
this precision regime the concrete definition of the top quark mass becomes relevant. The experiments 
measure a kinematic top quark mass which is approximately equal to the pole mass that appears in 
the perturbative top quark propagator. The ambiguity between pole and kinematic mass is on the 
order of 1 GeV and due to the fact that at this level the top quark cannot be treated as a free quark 
anymore. At a future linear collider an energy scan in the threshold region for tt production will allow 
a very precise measurement of the MS mass MtopifJ-)- An uncertainty of 20 MeV is envisaged. 

8. Top quark production and decay properties 

While the ti cross section and top mass measurements have established data samples that are 
compatible with SM top quark production, it is crucial to check whether the candidate events are 
in agreement with other predictions made for the SM top quark. In this chapter we discuss the 
investigation of various production and decay properties of top quarks: the helicity of W bosons 
from the top quark decay, the measurement of Rtb, the ratio of branching ratios for t ^ W -\- b and 
t ^ W + q, the search for the tt tau modes, the analysis of spin correlations among the tt pair, the 
measurement of the top quark px spectrum, and the search for electroweak top quark production. 
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8.1. W helicity in top quark decays 

As discussed in chapter 13.31 the V ~ A structure of the electroweak charged current interaction causes 
the W bosons from the top quark decay to be polarized. The fraction JFq of longitudinal W bosons 
(helicity hw = 0) is predicted to be 



^o= /::Z foZr2 -(71±1)%, (32) 



where we have used Aftop — (178.0 ±4.3) GeV/c^. For bosons the remaining 29% are left-handed 
{hw = —1), for bosons 29% are right-handed {hw = +!)• In the SM right-handed W'^ and 
left-handed W~ are strongly suppressed (branching fraction: 0.04%). In the following discussion we 
refer only to the , but imply the CP-conjugate statement for the . 



8.1.1. The lepton pT spectrum We consider further the leptonic decay of the W boson into eVe or 
/xf/i. The V — A structure of the W boson decay causes a strong correlation between the helicity of 
the W boson and the lepton momentum. Qualitatively, this can be understood as follows: The ve 
from the decay is always left-handed, the is right-handed. In the case of a left-handed 
boson angular momentum conservation demands therefore that the is emitted in the direction of 
the spin, that means anti-parallel to the momentum. That is why, charged leptons from 
the decay of left-handed W boson are softer than charged leptons from longitudinal W bosons, which 
are mainly emitted in the direction transverse to the W boson momentum. The spectrum of leptons 
from right-handed W bosons would be even more harder than the one from longitudinal ones, since 
they would be emitted preferentially in the direction of the W momentum. Figure 126b shows the 
PT distributions of leptons for different W helicities. The significant differences for the three helicity 
states are apparent and exploited by a CDF analysis which uses the dilepton and lepton-plus-jets 
data samples |252| . The threshold at 20 GeV/c^ seen in figure [^5k is due to the event selection. Out 
of the standard dilepton sample the helicity analysis uses only those events that feature leptons of 
different flavour, one electron and one muon. This additional requirement removes the major part of 
the Drell-Yan background. Seven dilepton events remain. The lepton-plus-jets sample is divided into 
three subsamples: (1) events with at least one secondary vertex b tag, (2) events with at least one soft 
lepton tag, but no secondary vertex tag, and (3) events with no b tags. In the first two subsamples 
at least three jets with Et > 15 GeV and at least one jet with Et > 8 GeV are required. In the 
third subsample at least four jets with Et > 15 GeV are mandatory. All jets are counted in the 
pseudorapidity region \r]\ < 2.4. 

An unbinned maximum likelihood fit to the observed lepton pT spectrum in each subsample is 
used to determine the fraction of top quark events that decay to longitudinal W bosons. The likelihood 
functions for each subsample are added and simultaneously minimized. In the fit function the lepton 
Pt spectrum for each subsample is parametrized as the product of an exponential and a polynomial. 
Figure 126b shows the observed lepton pT spectrum for the dilepton sample and the sum of all three 
lepton-plus-jets samples. The figure also shows the fit result (signal -I- background) and the background 
component alone. The likelihood fit yields a result of = 0.91±0.37 (stat.)±0.13 (syst.). The largest 
contribution to the systematic uncertainty is due to the uncertainty in the top quark mass, SJ^o — 0.07. 
The second largest source of uncertainty {5J-q = 0.06) is due to the normalization uncertainty of the 
non-W^ background which peaks at low px and thereby mimics the shape of the lepton pT distribution 
coming from negative helicity W bosons. Other sources of systematic uncertainties are the b tagging 
efficiency, Monte Carlo statistics, the modelling of gluon radiation, and uncertainties in the parton 
distribution functions. 

CDF has applied the same method as described above to measure J-q in Run II data. The new 
analysis is based on data samples corresponding to 162 pb~^ or 193 pb~^ for the lepton-plus-jets or the 
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Figure 26. (a) Transverse momentum distributions of Icptons from the W decay for the three 
different W helicities. The sohd dots are from bosons with negative helicity or W~ bosons 

with positive hehcity. The open circles are from longitudinally polarized W bosons. The closed 
squares are for leptons from right-handed bosons or from left-handed W~ bosons. All three 

distributions are normalized to the same area. The figure is taken from reference 12521 . (b) Lepton 
Pt distributions of the CDF Run I W helicity analysis in the lepton-plus-jets and the dilepton 
sample |252| . The data (dots with error bars) are compared with the result of the combined fit (solid 
line) and with the background component (dashed line). 



dilepton sample, respectively. The result for the combined sample is J^o = 0-3lio'23 (stat.)±0.17 (syst.) 
p53. 

8.1.2. Measurement of with the matrix element method D0 has used the matrix element method 
that is employed for the top quark mass measurement, see section [7.21 to determine the longitudinal 
polarization fraction |254j . The data sample consists of the same 22 events that are analyzed for 
the mass measurement. The tt probability density P^i[x,!F{)) includes contributions with longitudinal 
helicity (JFq) and negative helicity (J^-), but only the ratio Tq/T- is allowed to vary. The likelihood 
maximization yields J-q = 0.56 ± 0.31 (stat.) ± 0.07 (syst.). The statistical error also contains a 
contribution from the top quark mass uncertainty which is included in the measurement by integrating 
over the mass from Mtop = 165 GeV/c^ to Mtop = 190 GeV/c^. The result for J^o is in good agreement 
with the SM prediction. 

8.1.3. The helicity angle cosO^w Another possibility to measure the polarization of W bosons from 
top quark decays is to reconstruct the helicity angle OiWi which is defined as the angle between the 
lepton momentum in the W rest frame and the W momentum in the top quark rest frame. Leptons 
from the decay of longitudinally polarized W bosons have a symmetric angular distribution of the 
form 1 — cos^ 9iw ■ Leptons from left-handed W bosons have an asymmetric angular distribution of the 
form {\ — cos9iwYi while right-handed W bosons would have a distribution of the form {\ + cos9iwY ■ 
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The direct reconstruction of cos Og-^j^ requires information on the top quark and W rest frames 
and thus, the reconstruction of the respective momenta. Experimentally, the determination of ptop 
and pw is difficult, since it relies on the relatively poor measurement of jets and the ^t- However, 
there is a very good approximation which relates cos O^w to the invariant mass of the lepton and the 
b quark jet from the top quark decay: 

Ml = \ (Mtop - (1 + cos . (33) 

Using the quantity instead of cosOgw circumvents the difhculties related to direct top quark 
reconstruction. 

CDF has exploited relation to search for an anomalously large right-handed fraction of W 
bosons in Run I data [255] . The large top quark mass has led to speculations that the top quark could 
play an active role in electroweak symmetry breaking, which would lead to anomalous electroweak 
interactions of the top quark [2]. One possibility are left-right symmetric models that lead to a 
significant right-handed fraction of W bosons in top quark decays. An additional V + A component 
would lead to a lower left-handed fraction, but leave the longitudinal fraction unchanged. For 
the analysis CDF uses the same 7 dilepton events as for the lepton pT analysis mentioned earlier in 
section 18.1.1 1 Since the tagging of b quark jets is not employed in the dilepton sample, there are four 

combinations for each dilepton event. In the lepton-plus-jets sample at least one b tagged jet is 
required. There are 15 events with exactly one b tag and 5 events with two tags. The Mf^ distributions 
are fit to a linear combination of template distributions for the three W polarization states. Only the 
shape information is used, the normalization is left floating. The fit maximizes a binned likelihood as a 
function of fv+A, the right-handed fraction of the Wtb vertex. If there are two possible b jets that can 
be matched to the primary lepton, the fit is performed to two-dimensional distributions taking both 
solutions into account. The fit result is fv+A = — 0.211q;24 (stat.) ± 0.21(syst.) which is an unphysical 
value, but gives more preference to the SM V — A interaction rather than a,V + A contribution. The 
dominating systematic uncertainty is due to the uncertainty in the top quark mass, which contributes 
0.19 to the total systematic uncertainty of 0.21. The result can be converted into a one-sided upper 
limit of fv+A < 0.80 at the 95% confidence level (C.L.). If one assumes the longitudinal fraction 
to have the SM value this limit corresponds to < 0.24 at the 95% C.L.. This result can be 
combined with the measurement obtained from the lepton p^ analysis. The correlation of the two 
results is about 40%. The combination yields an improved limit of !F+ < 0.18 at the 95% C.L., which 
is inconsistent with a pure V + A theory at a confidence level equivalent to the probability of a 2.7 cr 
Gaussian statistical fluctuation. 

In Run II CDF has provided a preliminary measurement of J-'q using the method. The 
analyzed data sample consists of 31 lepton-plus-jets events with one secondary vertex b tag. The 
resuhing value is Tq = O.ggtfj j^ (stat.) ± 0.19 (syst.) [253] . 

8. 2. Measurement of Rtt 

In the SM the top quark is predicted to decay to a boson and a b quark with a branching fraction 
of nearly 100%. This is a consequence of the CKM matrix element \Vtb\ being close to unity. Our 
knowledge on | Vif,| comes primarily from measurements of b meson decays using the unitarity condition 
of the CKM matrix (VV^ = VtV = 1). This indirect method yields 0.9990 < \Vtb\ < 0.9992 with 
high precision [10]. However, if a fourth generation of quarks was present, unitarity of the CKM 
matrix could be violated. Therefore, it is desirable to make a direct measurement of the top quark 
branching fraction to Wb. 

In Tevatron Run I CDF has measured the ratio Rtb — BF{t ~* Wb)/BF{t Wq), where q 
can be any down- type quark [256] . The analysis uses 9 events from the dilepton sample and 163 
events from the lepton-plus-jets sample. The ratio Rtb is measured from the data by comparing 
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the observed number of events with secondary vertex b tags or soft lepton h tags to the number of 
expected events based on the kinematic acceptances, tagging efficiencies and background estimates. 
Four event categories are considered: (1) events with no h tags, (2) events with one or more soft 
lepton tags, (3) events with exactly one secondary vertex tag, and (4) events with two secondary 
vertex tags. In the dilepton sample only secondary vertex h tagging is used. A maximum likelihood 
fit yields Rtb = 0.941;q'2i (stat.) Iq 12 (syst.), consistent with the SM prediction. The result can also 
be expressed as a lower limit: Rtb > 0.56 at the 95% C.L.. 

The CKM matrix element |Vtfc| is related to Rtb, although in a model-dependent way. Under the 
assumption that the top quark decays only to Wq final states and using three generation unitarity 
one finds Rtb = \Vtb? ■ As a result, CDF finds \Vtb\ = 0.97lg;i^ or \Vtb\ > 0.75 at the 95% C.L. [256] . 

In Run II CDF has obtained a preliminary result of this analysis using only the secondary vertex 
b tag algorithm and improved the lower limit to Rtb > 0.61 at the 95% C.L. |257j . D0 has performed 
a similar analysis in the lepton-plus-jets sample using the secondary vertex and the impact parameter 
b tag method. The preliminary results are Rtb — 0.701q'24 (stat.) Iq'JJ (syst.) for the secondary vertex 
tagged sample and Rtb = 0.65~^q'^q (stat.) 12 (syst.) for events with an impact parameter tag [258j . 

8. 3. Search for fir and er top decays in ti events 

Up to now the searches for the tt tau decay modes have not produced enough evidence to claim the 
observation of this channel. We discuss in this paragraph the latest search by CDF using Run II data 
corresponding to an integrated luminosity of 194 pb~^ [259j . The search is carried out in the tt tau 
dilepton channel where one W boson decays into ei^e or /x;/^ and the second W boson decays into 
with a subsequent hadronic decay of the tau lepton. The branching fraction of this tt mode is about 
5%, the same size as the classical dilepton channel (ee, e^, /i/i). However, several factors lead to a 
reduced acceptance for the tau mode: the hadronic branching ratio of the tau is 64%, the kinematic 
acceptance is reduced due to the undetected neutrino, and the tau selection is less efficient than the 
electron or muon selection. In total, the event detection efficiency for the tau modes is about five 
times smaller than the one for the standard dilepton analysis. This fact mainly explains why the tau 
mode so far evaded observation. 

The data set for the CDF analysis was selected by a trigger on high-px electrons and muons. 
It is the same data set used for the dilepton cross section measurement discussed in section ISTTl On 
analysis level electrons and muons are identified using standard cuts as mentioned in the description 
of the dilepton analysis. Hadronic tau decays have a distinct signature of narrow isolated jets with low 
charged track multiplicity. A tau candidate is defined by two components: (1) A calorimeter cluster 
containing a seed tower with Et > 6 GeV and a maximum of five adjacent towers with > 1 GeV. 
(2) At least one track with px > 4.5 GeV/c is required to point to the tau calorimeter cluster. Further 
cuts are employed to reduce background. A cone is defined around the seed track using a variable 
radius of ^cono = min{0.17, (5 GeV)/i?ciustor} rad. Within this cone there must be exactly one or 
three tracks with px > IGeV/c including the seed track. In the case of three tracks, the charges 
must not all be positive or negative. Within an isolation annulus extending from the outer edge of 
the tau cone to 30° no tracks or 7r° candidates are allowed. The 7r° candidates are identified in the 
calorimeter by clusters of energy in the shower maximum detector. The transverse momentum of the 
tau is estimated by the sum of the track px plus the i?x /c of tt" candidates identified within the tau 
cone. Additional requirements on the tau candidate are px > 15GeV/c and JTicandidato < 1.8GeV/c^. 
The calorimeter towers within the isolation annulus are required to have Et less than 6% of the tau 
candidate -Ex- Electrons or muons that fake tau candidates are removed by asking that the energy 
in the hadron calorimeter divided by the sum of the tau track momenta is above 0.15, and that the 
Et of the tau calorimeter cluster divided by the tau seed track is above 0.5. Additional requirements 
on candidate events are: ^x > 20 GeV, and at least two jets within jryj < 2.0. The first jet has to 
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have i?T > 25 GeV, the second jet £'t > 15GeV. The scalar sum of the electron or muon pt, the 
tau pt, the : and the i?T of the jets is defined as TJt and must exceed 205 GeV. Including all 
branching ratios, kinematic cuts and efhciencies the event detection efficiency for ti events is found to 
be Cevt = (0.080 ± 0.015)%. 

The number of expected tt events is 1.00 ± 0.17. The estimated background is 1.29 ± 0.25 events. 
The dominant contributions to the background are W + jets events where a jet fakes a tau (58%), 
T~^T~ +jets (19%), WW production (11%), and Z° e"'"e~ +jets where an electron fakes a tau 
(6%). In CDF II data two events are observed, compatible with the SM expectation. The significance 
of the result is not high enough to claim the observation of the tt tau mode. The result can be 
converted into an upper limit on an anomalous enhancement of the t —^ TVrq rate. An enhancement 
of more than a factor of 5.2 is excluded at the 95% confidence level. The Run II analysis has improved 
the expected signal to background ratio over the previous Run I tt tau analysis of CDF [260j from 
0.5 to 0.78. With more data to come in Run II this improvement brings the observation of the tt tau 
mode within reach. In the Run I analysis 4 candidate events are observed over a background of 2 
events and a signal expectation of 1 event. 



8.4^. Spin correlations in tt events 



U Dili U l^UO W . X 

^^^^ = -. • 2., ■ (34) 

'J^A cm 



As discussed in detail in section 13731 top quarks decay as quasi- free quarks due to their short lifetime. 
Thus, the spin information of top quarks is transmitted to the decay products and is accessible 
experimentally. In this way top quarks are a unique laboratory to study spin aspects of heavy 
quark production, contrary to h quarks where the hadronization to h hadrons removes the initial 
spin information. At the Tevatron tt pairs are produced unpolarizcd. However, the spins of the top 
and antitop quark are expected to have a strong correlation and point along the same axis in the tt 
rest frame on an event-by-event basis. 

An optimal spin quantization axis can be constructed using the velocity /3* and the scattering 
angle 9* of the top quark with respect to the centre-of-mass frame of the incoming partons. The 
quantization axis forms the angle ip with respect to the pp beam axis [261 : 

sinr cosr 
1 - /3*2 sin^ 

For the qq ti process the spins of the top and the antitop quark are fully aligned along the same 
direction in this spin basis. If the gluon fusion process gg tt is included the correlation is reduced. 
In the limit where the tt pair is produced at rest (/3* = 0) the spins are pointing along the beam axis. 
The spin correlation of the tt pair is transmitted to the top quark decay products. Experimentally, 
it is particularly advantageous to investigate correlations of the primary leptons in tt dilepton events. 
The observables are the angles 6'_ and 6+ between the momenta of the negatively or positively charged 
leptons in the rest frame of their parent top quark and the spin quantization axis. Using these angles 
the spin correlation can be expressed by the differential cross section [262] 

1 d^cr 1 + KC0s9+ cos6- 

o-totai d(cos6i+)d(cos6l_) 4 

where k is the correlation coefficient which is predicted to be k = 0.88 at the Tevatron (Run I). 
The coefficient k can vary between —1, fully negative correlation, and k — +1 for a fully positive 
correlation. 

The D0 collaboration has searched for evidence of tt spin correlations in Run I data using dilepton 
events [263]. Since the angles 0_ and 9+ are measured in the top quark rest frame, the top quark has 
to be fully reconstructed. In dilepton events this cannot be achieved unambiguously, since the two 
neutrinos are not reconstructed. Up to four different solutions per event have to be considered. A 
weight is assigned to each solution based on how well the sum of the transverse momenta of the two 
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Figure 27. (a) Transverse momentum distribution of hadronically decaying top quarks from 61 
candidate events of reference 12641 . The hatched distribution shows the estimated background, the 
dashed distribution the SM prediction of background plus tt production, (b) Q ■ r) distributions for 
Monte Carlo templates normalized to unit area. The Q ■ rj distribution is used in the CDF Run II 
single top analysis to distinguish between t-channel and s-channel single top events |265| . 



neutrinos in the solution agrees with the measured ^t- For each solution the leptons are boosted into 
the respective top quark rest frame and cos0+ or cos0- are computed. 

To deduce k the two-dimensional phase space of (cos0+, cos0-) is used in a binned likelihood 
analysis. The bins are filled with the event weights obtained from the event fitter. Since an event 
populates each bin with fractional probability, the likelihood fit does not use a simple Poisson 
likelihood. The correlation between different bins is obtained from Monte Carlo events and the bin 
contents of the (cos^+, cos 6-) space are rotated such, that they are uncorrelated. The new variables 
are used to construct the likelihood. D0 uses six dilepton data events, three efi, two ee and one fifi 
event. The fit result from data is used to compute a lower limit of k > —0.25 at the 68% confidence 
level, which is in agreement with the SM expectation of k = 0.88. 

8.5. The top quark spectrum 

The CDF collaboration has further investigated the kinematics of tt events by measuring the top 
quark pT spectrum [264j . This analysis is partially motivated by exotic models that predict alternative 
production mechanisms for top quarks at the Tevatron, in particular an enhancement of top quarks 
with transverse momentum above 200 GeV/c. The Run I lepton-plus-jets data set with secondary 
vertex b tag or soft lepton tag are used in the analysis. The events are subjected to a kinematic fit 
similar to the top quark mass analysis described in section mi The invariant mass reconstructed from 
the top quark decay products is constrained to 175GeV/c^. Events with > 10 from the fit are 
rejected. The final data sample contains 61 candidate events. Figure [27k shows the pt distribution of 
the hadronically decaying top quark. The differences in the spectrum do not constitute a significant 
deviation from the SM prediction. A maximum likelihood technique is used to estimate the fraction 
i?4 of top quarks in the high transverse momentum region of 225 < pt < 425 GeV/c. The measured 
fraction is i?4 = O.Oto'o^'i while the SM expectation is i?4 = 0.025. The data result can be converted 
into an upper limit of R4 < 0.16 at the 95% confidence level. 



Top quark physics in hadron collisions 



71 



8. 6. Electroweak top quark production 

While tt pair production via the strong interaction is the dominant source of top quarks at the 
Tevatron (and also at the LHC), top quarks can also be produced as single quarks via the electroweak 
interaction. The theoretical aspects of single top production are discussed in section 13.21 of this 
review. The relevant modes are i-channel and s-channcl production, see figure 1101 Several analyses 
have searched for electroweak top quark production in Tevatron Run I data |2661 12671 12681 1269j . The 
best upper limits obtained on the production cross sections in Run I are 13 pb at the 95% C.L. for 
the t-channel [268j and 17 pb at the 95% C.L. for the s-channel |267j . 

The CDF collaboration has published the first search for single top production in Run II 265 . In 
the following paragraphs we will briefly review this analysis which uses a data sample corresponding 
to an integrated luminosity of (162 ± 10) pb~^. The event selection is essentially the same as the one 
used for the tt cross section measurement in lepton-plus-jets events with secondary vertex b tagging, 
see section 16.2.21 The jet definition differs in the pseudorapidity range. While the tt cross section 
measurement counts jets up to \r]\ < 2.0, the range is enlarged in the single top analysis to \r]\ < 2.8. 
This change is motivated by an increase in acceptance for the t-channel process in the W^-|-2 jets sample 
by roughly 30%. Exactly two jets are required. At least one of these jets must be identified as a 5 quark 
jet by a secondary vertex tag. To optimize the sensitivity, CDF applies a cut on the invariant mass 
Me„b of the charged lepton, the neutrino and the b tagged jet: 140 GeV/c^ < Mi„b < 210 GeV/c^. The 
transverse momentum of the neutrino is set equal to the missing transverse energy vector Pz (i') is 
obtained up to a two- fold ambiguity from the constraint M^^ = Mw- From the two solutions the one 
with the lower |p2(i^)| is chosen. If the Pzi^) solution has non-zero imaginary part as a consequence 
of resolution effects in measuring jet energies, only the real part of Pz(;^) is used. Two analyses are 
performed on this data sample: (1) a separate search, which measures the rates for the two single top 
processes, t-channel and s-channel, individually, and (2) a combined search where t- plus s-channel 
are treated as one single top signal. For the separate search, the data sample is subdivided into events 
with exactly one b tagged jet or exactly two b tagged jets. For the 1-tag sample, at least one jet is 
required to have Et > 30 GeV. 

The event detection efficiency for the signal is determined from events generated by the matrix 
element event generator MadEvent [2701 1271j . followed by parton showering with Pythia |171j 
and a full CDF II detector simulation. MadEvent features the correct spin polarization of the top 
quark and its decay products. For t-channel single top production two samples are generated, one 
b + q ^ t + q' and one g + q^t + b + q' which are merged together to reproduce the px spectrum 
of the b as expected from NLO differential cross section calculations. This is an improved model 
compared to the Pythia modelling used in the Run I analyses. The event detection efficiency Cevt 
includes the kinematic and fiducial acceptance, branching ratios, lepton and b jet identification as 
well as trigger efficiencies. In the 1-tag sample one finds ecvt = (0.86 ± 0.07)% for the t-channcl and 
e^^j = (0.78 ± 0.06)% for the s-channel. 

Two background components are distinguished: ti and nontop background. The nontop 
background is estimated using the same method as used for the ti cross section analysis, see 
section [^.2.2l The primary source (62%) of the nontop background is the W+hea,vy flavour processes 
qq' Wg with g —f bb oi g ^ cc, and gq Wc. Additional sources are "mistags" (25%), in which a 
light quark jet is erroneously identified as heavy flavour, "non-H^" (10%), e.g. direct bb production, and 
diboson {WW, WZ, ZZ) production (3%). The non-W and mistag fractions are estimated using CDF 
II data. The W+heavy flavour rates are extracted from Alpgen |178j Monte Carlo events normalized 
to data. The diboson rates are estimated from Pythia events normalized to theory predictions [272] . 
After all selection cuts CDF observes 33 events in the 1-tag sample, 6 events in the 2-tag sample, and 
42 events for the combined search. Within the uncertainties, the observations are in good agreement 
with predictions. 
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Figure 28. Search for electroweak top quark production at CDF in Run II [265] . (a) Q-rj distribution 
for CDF data in the 1-tag sample (33 events) versus stacked Monte Carlo templates normalized to 
the SM prediction, (b) /f-p distribution for data (42 events) in the combined single top search, 
compared with smoothed Monte Carlo predictions for signal and background. The inset shows 
the //t distributions for t-channel and s-channel single top production to demonstrate that these 
distributions are very similar and may be combined. 



To extract the signal content in data, a maximum likelihood technique is employed. Separation of 
t- and s-channel events is achieved by using the Q ■ 77 distribution which exhibits a distinct asymmetry 
for t-channel events, see figure [?7b . Q is the charge of the primary lepton and 77 is the pseudorapidity of 
the untagged jet. Figure [28b shows CDF data versus stacked Monte Carlo templates weighted by the 
expected number of events in the 1-tag sample. The separate search defines a joint likelihood function 
for the Q ■ rj distribution in the 1-tag sample and for the number of events in the 2-tag sample. 
The background rates are constrained to their SM prediction by Gaussian priors. The systematic 
uncertainties are taken into account in the likelihood definition. Systematic shifts in the acceptance 
and in the shape of the Q ■ rj template histograms, and their full correlation are considered. The 
largest uncertainties on the acceptance are due to the uncertainty on the b tagging efficiency (7%), the 
luminosity (6%), the top quark mass (4%), and the jet energy scale (4%). The likelihood procedure 
yields a most probable value of O.O^Q Qpb for the t-channel cross section, and (4.6 ± 3.8) pb for the 
s-channel cross section. These results are translated into upper limits on the cross sections which 
exclude an anomalous enhancement of single top quark production. The upper limit for the t-channel 
is found to be 10.1 pb at the 95% C.L., the upper limit for the s-channel is 13.6 pb at the 95% C.L.. 

To measure the combined t-channel plus s-channel signal in data, a kinematic variable is used 
whose distribution is very similar for the two single top processes (see inset of figure [28b ). but is 
different for background processes: Ht, which is the scalar sum of I^t and the transverse energies of 
the lepton and all jets in the event. In figure j^Hb the distribution observed in data is compared 
with the SM prediction. A likelihood function similar to that of the separate search yields a most 
probable value of 7.7^4 gpb for the single top cross section. The resulting upper limit is 17.8 pb at 
the 95% C.L.. The quoted result (combined search) excludes an anomalous enhancement of single 
top production which is more than 6.1 times larger than the SM production. However, separate 
measurements of the t-channel and s-channel are important because the two processes are sensitive 
to different new physics contributions, see for example references |273i 12741 12751 [2761 12771 [278] . 

D0 has presented results on a single top search using neural networks j279j . With a data sample 
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corresponding to an integrated luminosity of 230 pb^^ D0 can set an upper limit of 5.0 pb at the 
95% C.L. on the t-channel cross section and 6.4 pb at the 95% C.L. on the s-channel cross section. 
With more data at the Tevatron to come, a first observation of electroweak top quark production 
is in reach within the next years. Preliminary studies indicate that data corresponding to about 
1.5 fb~^ will allow to establish evidence for electroweak top quark production at a probability level 
excluding a background fluctuation that is equivalent to a 3 ct deviation in a Gaussian distribution. 
Once electroweak top quark production is observed the measurement of the production cross section 
will allow to directly extract the CKM matrix element |Vt6|. 

9. Search for anomalous couplings 

Mainly because of its large mass the top quark has fostered speculations that it offers a unique window 
to search for physics beyond the SM. One possibility to test this hypothesis is to measure top quark 
properties and check whether the observation is in agreement with the SM prediction. This approach 
is based on the measurements that are discussed in previous chapters, involving properties such as the 
tt cross section, the top mass, or the W helicity in top quark decays. A second approach is to search 
directly for new particles coupling to the top quark or for non-SM decays. In this section we present 
those topics which have led to specific analyses at the Tevatron. In section [93] we will additionally 
discuss searches for anomalous single top quark production at HERA and LEP. 

9.1. Decays to a charged Higgs boson 

In the SM a single complex Higgs doublet scalar field is responsible for breaking the electroweak 
symmetry and generating the masses of gauge bosons and fermions. Many extensions of the SM 
include a Higgs sector with two Higgs doublets and are therefore called Two Higgs Doublet Models 
(THDM). In a THDM electroweak symmetry breaking leads to five physical Higgs bosons: two neutral 
CP-even scalars h'^ and H'^, one neutral CP-odd pseudoscalar A", and a pair of charged scalars H^. 
The extended Higgs sector is described by two parameters: the mass of the charged Higgs, Mh+ , 
and tan/? — vi/v2, the ratio of the vacuum expectation values vi and V2 of the two Higgs doublets. 
One distinguishes two types of THDMs. In a type 1 THDM only one of the Higgs doublets couples 
to fermions, in a type II THDM the first Higgs doublet couples to the up-type quarks {u, c, t) and 
neutrinos, while the second doublet couples to down-type quarks {d, s, b) and charged leptons. The 
analyses we discuss in this section are concerned with type II models. A particular example for a type 
11 THDM is the minimal supersymmetric model (MSSM). 

If the charged Higgs boson is lighter than the difference of top quark and b quark mass, 
Mh± < Aftop — Ml,, the decay mode t H'^b is possible and competes with the SM decay t W'^b. 
The branching fraction depends on tan/3 and M^+. The MSSM predicts that the channel t H^b 
dominates the top quark decay for tan (3 1 and tan /3 > 70. In most analyses it is assumed 
that BF{t W+b) + BF{t H+b) = 1. At tree level the does not couple to vector bosons. 
Therefore, the decays only to fermions. In the parameter region tan/3 < 1 the dominant decay 
mode is cs, while for tan/3 > 1 the decay channel t+j^t- is the most important one. For 

tan/3 > 5 the branching fraction to T^Ur is nearly 100%. Thus, in this region of parameter space type 
11 THDM models predict an excess of tt events with tau leptons over the SM expectation. 

First searches for the in top quark events were already performed well before the top quark 
was discovered in 1994/95. The UAl and UA2 experiments at the CERN SppS excluded certain 
regions of the Mtop versus Mh± plane [2801 1281j . First searches of CDF improved these limits using 
events with a dilepton signature [282] or reconstructing tau leptons in their hadronic decay mode |283j . 

At LEP the experiments searched for charged Higgs production in the process e+e~ H'^H~ . 
The analyses encompass the final states cscs, t+i^^t^^v, cst^ i7r, and t^Vt-cs. The combined result of 
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all four experiments excludes charged Higgs masses below 78.6 GeV/c^ |284| . The CLEO collaboration 
has used its measurement of the inclusive 6 — > 57 cross section to set an indirect limit on the charged 
Higgs mass [285] . The dependence of the cross section on Mh± enters via quantum corrections. Since 
the process h ^ occurs only at loop level in the SM, its sensitivity to new particles which might 
be exchanged in the loop is quite high. CLEO sets a limit of Mh± > (244 + 63/(tan/3)^-^) GeV/c^ at 
the 95% C.L., under the assumption that a type II THDM is realized in nature. If the Higgs sector 
has a richer structure, the CLEO limit can be circumvented. Therefore, direct searches for the 
which are less model dependent remain important. 

After the top quark discovery the searches at the Tevatron looked for the decay modes 
tt H^W^bb and tt H^H^bb. CDF published three analyses where the charged Higgs is assumed 
to decay into t~^Vt and the tau lepton subsequently decays semi-hadronically into a tau neutrino plus 
hadrons [286| I287| I288j . The identification of hadronic tau decays in these CDF analyses is similar 
to the one described in detail in section [831 Two analyses (references |286| 1287] ) used an inclusive 
approach, based on final states with missing transverse energy, one identified hadronic tau decay, jets 
or high-pT leptons. One analysis (reference |288| ) requires one high-px lepton (e or /i) in the event, 
which can originate from the decay W —>■ e/n + v or — > t^Vt — > P^. + e+//x+ + v^/^ + Vr- The 
best exclusion limit is obtained from the inclusive analysis |287) : For M//± = 80GeV/c^ the region 
of tan/3 > 25 is excluded. For higher Higgs masses the exclusion becomes less stringent, reaching 
tan/3 > 200 at Mh± = leOGeV/c^. 

The first D0 search for the is based on a disappearance strategy, looking for a deficit of 
events in the SM lepton-plus-jets channel [289]. Since good agreement between the SM prediction and 
the observation is found, certain regions of phase space for charged Higgs production can be excluded. 
This result is superseded by a direct search for hadronic tau decays in top quark events |290j . In the 
following paragraph we will discuss this last Run I analysis on the charged Higgs in more detail. 

The D0 analysis uses data taken with a multijet + trigger corresponding to (62.3±3.1) pb~^ of 
integrated luminosity. Relatively loose preselection cuts are applied, followed by tighter cuts involving 
an artificial neural network. The preselection requires that > 25 GeV, at least four jets with 
E'Y > 20 GeV each, but no more than eight jets with > 8 GeV. The selection is sensitive to 
ti H^W^bb events where the W boson decays into quarks as well as tt H'^H^bb events. 
The network has three input variables: the and two of the three eigenvalues of the normalized 
momentum tensor. The training is done with charged Higgs events generated by Isajet |291| for the 
signal and with 25000 multijet events from data for the background. Figure [29k shows the number of 
remaining events in data versus the cut on the neural network output. Overlaid is the number of events 
expected from the SM as well as one scenario including a charged Higgs boson with M^± = 95 GeV/c^. 
The data show good agreement with the SM expectation, thereby falsifying the chosen scenario. 
Based on a series of Monte Carlo experiments the sensitivity for the search is optimized yielding 
on optimal cutoff at 0.91. 

After the neural network cut hadronic tau decays are identified in the data. A narrow jet with 
AR < 0.25 is required, with 1 to 7 tracks in the jet cone. Events with identified muons or electrons 
are rejected. In addition, a requirement is constructed using calorimeter information from W —> tv 
Monte Carlo events. After all cuts 3.2 ± 1.5 multijet background events are expected, 1.1 ±0.3 from tt 
and 0.9±0.3 from H^+jets. In data D0 observes three events. A Bayesian and a Frequentist statistical 
analysis are used to deduce exclusion regions is the {Mh± , tan /3) phase space shown in figure [29b . 
Under the assumption that Mtop — 175GeV/c^ and a{tt) = 5.5 pb, the region of tan/3 > 32.0 is 
excluded at the 95% C.L. for Mh± = 75GeV/c^. For Higgs masses above 150 GeV/c^ no limits can 
be set. This result can also be interpreted in terms of the branching fraction of i ^ setting a 

lower limit of BF(t -> H+b) > 0.36 at the 95% C.L.. 
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Figure 29. Search for the at D0 in Run I |290) . (a) The number of observed events in data 
versus the cut on the neural network output (dots). Overlaid is the expectation from the SM (hght 
shaded area) and a scenario with additional production in top quark decays for tan f) = 150 and 
Mjj± = 95GeV/c2. (b) Exclusion region at the 95% confidence level in (M^±,tan/3) space. The 
top quark mass is assumed to be 175 GeV/c^ and cr{tt) = 5.5pb. The plot also contains the regions 
excluded by the indirect search method I289| . where D0 looks for a potential disappearance of the 
tt signal. 



9. 2. Search for X'^ —> tt decays 

Several extensions of the SM predict the existence of narrow resonances that decay to tt pairs. One 
such model is, for example, a Z' predicted by top-colour-assisted technicolour |292j . This model 
speculates that the spontaneous breaking of electroweak symmetry is related to the observed fermion 
masses, in particular the large top quark mass, and can be accomplished by dynamical effects |293j . 

CDF and D0 have performed model-independent searches for a narrow resonance X° which 
decays into a tt pair [2941 I295j . Both analyses relate closely to the respective top quark mass 
measurements of the two experiments j224l 1243] and employ fits to the invariant mass spectrum 
of the tt pair. D0 obtains a slightly better mass limit. The D0 analysis is based on data with 
an integrated luminosity of 130 pb~^ and uses events with a lepton-plus-jets signature passing the 
topological and the soft muon b tag selection. The events are subjected to a kinematic fit that 
comprises kinematic constraints on the reconstructed W boson mass and the top quark mass. Out of 
all possible permutations the option with the best is chosen. The is required to be below 10. 
After all cuts 41 events are left in the data sample. Figure[5Dk shows the invariant mass distribution of 
the tt pair reconstructed from data. The data are fit to distributions for the combined QCD multijet 
and W+ jets background, SM tt production, and a signal for X'^ tt. For comparison to the data 
figurelSOb shows the backgrounds and the signal component {Mx — 400 GeV/c^) obtained from the fit. 
There is no significant deviation from the SM prediction. The fit is repeated for several assumptions 
for the X'^ mass, again without providing evidence for a signal. Therefore, upper limits on the X^ 
production cross section times branching ratio to ti are set, ranging from 5.0 pb at Mx = 400 GeV/c^ 
to 1.5 pb at Mx = 850GeV/c^. The confidence level is 95%. These limits are valid as long as the 
width of the resonance, Fx, is small compared to the D0 detector resolution. The limits quoted above 
assume a width of Fx — 0.012 Mx. If one compares the cross section limits with the predictions for 
a leptophobic Z' boson, one can obtain a lower limit on the Z' mass: Mz' > 560 GeV/c^. 




Figure 30. (a) Invariant mass distribution of the reconstructed tt pair of 41 events in the D0 
lepton-plus-jets data sample 12951 . The solid histogram shows the shape of the W+ jets and QCD 
multijets background, the hatched histogram the sum of all SM backgrounds including ti production, 
and the open histogram shows the sum of signal (Mx = 400GeV/c^) and background as obtained 
from the fit to the data, (b) Feynman diagram for the FCNC decay t — > cZ" with Z ^ e+e" . Other 
diagrams involve the d and the s quark within the loop instead of the b quark. All diagrams together 
nearly cancel and cause the branching ratio of this decay to be very small, 0(10"^^), in the SM. 



9.3. FCNC decays 

In the SM flavour changes at leading order (tree level) are induced by charged currents, the exchange 
of W bosons. Flavour changing neutral currents (FCNC) are only possible at higher orders in 
perturbation theory (loop level). FCNC decays of the top quark are strongly suppressed in the SM 
due to the Glashow-Iliopoulos-Maiani (GIM) mechanism j296| . Compared to the top quark mass the 
masses of down-type quarks (d, s, b) occurring in loop diagrams are small and degenerate. Therefore, 
the sum over the respective amplitudes nearly cancels. Fig. 130b shows one of the Feynman diagrams 
that describe the decay t cZ'^. This type of diagram is also called penguin diagram. The SM predicts 
the branching fractions for t — > and i — > 07 to be on the order of 10~^^, while BF(t eg) ~ lO"""^" 
and BF(t cH°) ~ IQ-^^ [M]- 

In extensions of the SM FCNC top quark decays can be considerably enhanced by several orders 
of magnitude if FCNC couphngs at tree level are allowed ' ^55112^511500] . In Two Higgs Doublet Models 
where the neutral scalar h'^ posses flavor changing couplings the decay t — > ch^ can be considerably 
enhanced |301j . A similar enhancement can be reached in supersymmetric models were i?-parity is 
violated |302| . In topcolour-assisted technicolour theories the branching ratio for t cW^W~ can 
reach values up to 10"'^ [303] . Since the SM predictions for FCNC interactions have much smaller 
rates, they are useful probes for new physics beyond the SM. 

The CDF collaboration has searched for the FCNC decays t qZ'^ and t — > 57 in Run I data 
with an integrated luminosity of 110 pb~^ j304j . The Z° is measured in the decay channels e+e" or 
Electron and muon identification is described in section [4.5.11 A Photon is identified as an 
energy cluster in the electromagnetic calorimeter with no track pointing at it. Additionally, a cluster 
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with a single soft track pointing to it is accepted as a photon, if the track carries less than 10% of 
the cluster energy. A typical photon cluster consists of two adjacent towers in the electromagnetic 
calorimeter. The direction of the photon is defined by the line between the event vertex and the 
centroid of the electromagnetic shower as measured in the shower maximum detector. Background 
from hadronic jets is reduced by demanding that the photon cluster is isolated from other energy 
depositions in the electromagnetic or the hadron calorimeter. In addition, the ratio of hadronic to 
electromagnetic energy is constrained to be less than 0.055 + 0.00045 E, where E is the sum of hadronic 
and electromagnetic energy of the cluster measured in the calorimeter. 

The search strategy for the FCNC decays tries to isolate tt events where one top quark decays 
via FCNC and the second top quark decays according to the SM decay t — > W^h. In the t — > 57 
search two signatures are considered, where the W boson decays either leptonically into eVe or ^i/^ 
or hadronically into quarks. In the first subsample a well identified charged lepton (e or with 
PT > 20 GeV/c, $T > 20 GeV, a photon with Et > 20 GcV, and at least two jets with Et > 15 GeV 
are required. The second subsample contains events with at least four jets and one photon with 
Et > 50 GeV. At least one jet must be identified as a 6 quark jet with a secondary vertex tag. In 
both subsamples, there must be a photon-plus-jet combination with an invariant mass in the top 
quark mass window 140 < Mj^ < 210GeV/c^. In the hadronic subsample, the remaining jets must 
have ^ Et > 140 GeV. One event passes all selection criteria. The event features a 72 GeV muon, 
= 24 GeV, a 88 GeV photon, and three jets. The background is expected to be about 0.5 events, 
mainly from production with two additional jets. 

In the t qZ^ search, the W from the second top quark decay is considered to decay 
hadronically, while the is reconstructed in its leptonic decay modes to e+e^ or /i^/i^. This 
yields an event signature of four jets and two leptons with an invariant mass close to the Z^ mass 
(75 < Mu < 105GeV/c^). Since the branching ratio of i'^i" is small, the total detection 

efficiency of the t — > qZ^ mode is much smaller than the one for the t q-f search. Each of the four 
jets must have Et > 20 GeV and |ry| < 2.4. One Z M^/^^ event passes all selection criteria. The 
expected background is 0.6 events. 

Since no statistically significant excess is observed in both searches, the data are used to set upper 
limits on the top quark branching ratios into FCNC decays 

BR(i ^ w/c + 7) < 3.2% BR (t ^ u/c + Z") < 33% (36) 

at the 95% confidence level. These limits are still far from the interesting region. At the LHC 
experiments will reach a sensitivity level of 10~^ [5], and thereby reach a region where some models 
predicting very large enhancements can be excluded |305| . 

9.4- Anomalous single top production 

At the Tevatron the search for anomalously enhanced FCNC in top quark decays is statistically 
restricted by the number of tt pairs produced. More stringent limits on anomalous tq"f of tqZ couplings 
at tree level are set by searching for the production of single top quarks via FCNC at LEP and HERA. 
At LEP all four experiments have searched for the reactions e~^e~ tc/tu and presented upper limits 
on anomalous couplings |306[ 13071 13081 1309j . The limits obtained by the L3 collaboration are shown 
in figure [311 For K-y = the coupling to the tqZ couphng is constrained by kz < 0.41, for = 
the tq'y coupling is confined to K-y < 0.49. Both limits are taken at the 95% confidence level. The 
other LEP collaborations obtain very similar results. Figure |3T] also visualises the limits that result 
from the CDF analysis discussed in section 19.31 The LEP results can be translated into limits on the 
branching ratios: BR(t u/c 7) < 4.2% and BR{t u/c + Z") < 14% [505] . 

At HERA the two experiments HI and ZEUS searched for top quarks produced in the inclusive 
FCNC reaction ep etX |3101 1311] . The SM process for single top production at HERA is the 
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Figure 31. Limits at the 95% confi- 
dence level on anomalous couplings 
Ktii7 and KtuZ (in this plot noted 
as Vfuz) of the top quark to the 7 
or the Z boson from reference 13101 . 
The plot shows the limits obtained 
by the HERA experiments HI and 
ZEUS [3IOII3TT] . by L3 [307] and by 
CDF [304] (see section [931. 




Figure 32. Feynman diagrams 
for single top production via 2 — > 
2 FCNC processes in hadron 
collisions. 



charged current reaction ep vthX [312j with a predicted cross section of less than 1 fb [113] , well 
below the sensitivity of the experiments. An observed excess in this channel can therefore signal 
physics beyond the SM. ZEUS has analyzed a data sample corresponding to 130 pb^^ of integrated 
luminosity. The top quark is assumed to decay according to the SM decay t W^b. The W boson is 
detected in its leptonic and the hadronic decays channels. The observed data are in good agreement 
with the background expectation. Thus, ZEUS is able to compute an upper limit on the anomalous 
twy coupling: Ktu-y < 0.174 at the 95% confidence level, which is the best limit on this quantity to 
date. The ZEUS result is depicted in figure [511 

The HI analysis uses data with an integrated luminosity of 118.3 pb~^ [310] . In the leptonic 
channel five events are observed over a SM background of 1.31 ± 0.22. If this excess is attributed 
to FCNC top quark production a total cross section of a{ep etX) = 0.29^Q'j4pb is found. In 
the hadronic channel there is no excess in the data. Only more data will help to decide whether the 
observation is a real signal or just a statistical fluctuation. 

At the Tevatron there are no searches for anomalous single top quark production yet. Possible 
FCNC processes at tree level could be qiqi tq2, qiqj tq2, qg — > tg, and gg tq. Example 
Feynman diagrams for these processes are shown in figure [32l Another option for FCNC top quark 
production is the 2^1 quark-gluon fusion process g + u/c t. Studies at parton level suggest 
that the anomalous up quark coupling parameter Kt„g/A can be measured down to 0.02 TeV""'^ in 
Run II of the Tevatron [275j . For the anomalous charm quark coupling parameter Ktcg/A a bound of 
0.06 TeV~^ is in reach. Studies with a fast detector simulation for ATLAS indicate a similar reach at 
the LHC ISni- 
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10. Summary and outlook 

Ten years after the discovery of the top quark many of its properties have been thoroughly investigated. 
All measurements are in good agreement with the predictions made by the Standard Model (SM). 
At the Tcvatron and at LHC the production of tt pairs is the dominating source of top quarks. The 
experiments CDF and D0 have established tt signals in the dilepton, the lepton-plus-jets and the all 
hadronic channels and measured the production cross section. The tt channels involving tau leptons 
have not been observed yet, but they arc well within reach in Run II at the Tcvatron. 

The most important parameter of the top quark is its mass. The large value of M^^p distinguishes 
the top quark strongly from the other quarks, since it decays essentially as a free quark. Thus, the 
top quark is an ideal laboratory to study polarization effects in heavy quark production and decay. 
Via higher order perturbative corrections the top quark mass is strongly correlated with the W boson 
mass and the Higgs boson mass. Therefore, electroweak precision tests and predictions of Mh are 
very sensitive to Mtop. The final combined Run I result on the top quark mass yields: Mtop = 
(178.0 ± 4.3) GeV/c^. Run II measurements are underway. First preliminary measurements indicate 
a slightly lower mass value. CDF finds Mtop = 173.5l|^ (stat.) ± 2.5 (JES) ± 1.7 (syst.) GeV/c^ D0 
Mtop = 170.6 ± 4.2 (stat.) ± 6.0 (syst.) GcV/c^. At the end of Run II the top quark mass will be 
measured with an uncertainty of 2 to 3GeV/c^. At the LHC the precision of Mtop will reach about 
1 GeV /(?. A limiting factor will be the theoretical understanding of the relation between the kinematic 
quark mass measured by experiments and the theoretically well defined pole mass or MS mass. The 
final word on Mtop will come from an e^e^ linear collider where a scan of the tt production threshold 
will allow for a precision of about 20MeV/c^. 

After the discovery of the top quark the focus of research by the CDF and D0 collaborations 
shifted to a detailed investigation of its production and decay properties. The hclicity of W bosons 
from the top quark decay was measured and found to agree with the SM V—A theory of the electroweak 
interaction. The ratio of branching ratios Rtb = BF(t Wb)/BF(t Wq) was determined and found 
to agree with a value close to 1 as predicted by the SM. D0 has investigated ti dilepton events for spin 
correlations among the two charged leptons, CDF has measured the top quark pT spectrum. Both 
analyses show no significant deviations from the SM. The search for electroweak single top quark 
production is about to reach a critical phase in Run II. A discovery of this top quark production mode 
is around the corner. 

Its large mass is motivation to view the top quark as an ideal probe for new physics beyond the 
SM. Already before the top quark discovery the search for the decay t H^h started. CDF and D0 
searched also for a heavy resonance which decays into a ti pair. A CDF analysis looked for FCNC 
decays of the top quark. No indication of these phenomena beyond the SM is foiuid. 

The LHC will be a top quark factory where several millions of top quarks will be produced per 
year. This wealth of data will allow to enter a precision era of top quark physics where top quark 
properties will be very accurately investigated. The Tcvatron analyses provide an excellent basis for 
further studies at the LHC. But also new topics where huge amounts of data are needed will become 
accessible at the LHC, for example, the search for like sign top pair production {tt or S), production of 
single top quarks via FCNC, the search for CP violation in production, or the direct measurement of 
the electric charge of the top quark. The abundance of analysis topics and open questions guarantees 
that top quark physics will remain a lively and thrilling field of elementary particle research in the 
coming years which will see the end of Run II at the Tevatron and the turn on of the LHC. 
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